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PREFACE. 



The preparation of this volume was undertaken with the 
view of presenting in as simple and concise a manner as 
possible the fundamental principles of the Method of Least 
Squares. While it is believed that everything essential to 
the solution of all ordinary problems has been included, no 
attempt has been made to develop at length those special 
methods and forms that are so useful and almost necessary in 
case large numbers of observations of certain kinds, such, 
for instance, as those met with in geodetic and astronomical 
measurements, are to be adjusted. 

Frequent references throughout the text, and more particu- 
larly the list of works given on page v of the Appendix, wUl, 
however, enable the student to extend his studies in what- 
ever special direction his profession may require ; it being 
expected that this book will in such cases be looked upon 
merely as an introductory treatise. All of the works men- 
tioned have been freely consulted in the preparation of these 
pages, and the author desires in particular to acknowledge 
his indebtedness for many of the examples. 

Dana p. Baetlbtt. 



Cornell University 
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tine Cornell University Library. 
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THE METHOD OF LEAST SQUARES. 



OHAPTEK I. 

GENERAL PRINCIPLES. 

1. In scientific investigations of all kinds it is frequently 
necessary to determine the values of certain quantities by 
means of actual measurements either with or without the aid 
of instruments. The observations may be made directly upon 
the values of the unknown quantities or upon certain functions 
of the unknowns. In the latter case the values of the required 
quantities must be obtained by computation from the observed 
values of the functions. In order to obtain more accurate 
values of the unknowns than would be given by a single 
measurement, or set of measurements, the observations are 
usually repeated either in the same way and under the same 
conditions or in a .variety of different ways and under vary- 
ing conditions. 

Under these circumstances it will invariably be found that 
the different measurements give discordant results, the amount 
of the discrepancies varying with the character of the observa- 
tions ; and the question that now presents itself is how to 
determine from these discordant observations the true values 
of the required quantities. From the nature of the case, 
however, we can not expect to obtain our values with absolute 
accuracy; all that we can hope for is to obtain those values 
which are rendered most probable after all the observations 
are taken into account, and, further, to determine the degree 
of confidence that can be placed in those values. 
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2. The attainment of the above results constitutes the 
primary object of the Method of Least Squares. The method 
is also employed in comparing the relative worth of different 
measurements of the same quantity, and in determining the 
equation of a curve which shall suitably represent the relation 
between two variables in cases where the exact law connecting 
them is not known. 

Also, before making any observations, we may employ the 
method to determine how precise the component measurements 
of a series must be in order to yield a required degree of pre- 
cision in the final result; or, conversely, to determine what 
the precision of the final result will be, knowing the precision 
attainable in the component measurements. This latter appli- 
cation of the method will be treated at length in the course 
on " The Precision of Measurements." 

3. Errors. The cause of the discrepancies between the 
results of our different observations is that every observation 
that is a measure is subject to error. These errors are of two 
kinds, — Constant or Systematic Errors and Accidental Errors. 

4. Constant Errors are errors which in all measures of the 
same quantity, made with the same care and under the same 
conditions, have the same magnitude, or whose presence and 
magnitude are due to some fixed cause. These constant 
errors may be of several classes, which are designated as 
follows: — 

First. Theoretical Errors, such as those due to the refrac- 
tion or aberration of light, the effect of a definite change in 
temperature or moisture on our standards of measurement, 
etc. As soon as their causes are known the magnitude of 
these errors may be calculated and their effect eliminated 
from the observations. 

Second. Instrumental Errors, such as errors of division of 
graduated scales, defects in micrometer screws, eccentricity of 
circles, etc. These errors will be discovered by an examina- 
tion of the instruments and their effects eliminated from the 
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observations, either by a particular method of using the 
instruments or by subsequent computation. 

Third. Personal Errors. These are due to personal pecu- 
liarities of an observer, who always answers a signal too soon 
or too late, always estimates a quantity smaller than it is, etc. 
The character and magnitude of these errors may be deter- 
mined by a study of the observer, his "Personal Equation" 
may be obtained, and his observations thus corrected for this 
source of error. 

5, Mistakes. Although of a somewhat different character, 
these should be considered in connection with constant errors. 
A mistake is made when a figure 3 is read for a figure 8, or 
when in reading a graduated circle which is numbered in both 
directions the angle is read 43° instead of the complementary 
angle 47°, etc. These mistakes are usually of such a charac- 
ter that they may be detected by an inspection of the observa- 
tions and a proper correction made. 

6. Accidental Errors are errors due to irregular causes, 
whose effect upon the observations is not determined by any 
circumstances peculiar to that particular set of measurements, 
and which cannot therefore be computed and allowed for 
beforehand. Such errors are those due to sudden changes in 
refraction owing to sudden and unobserved changes in tem- 
perature; unequal expansion of different parts of an instru- 
ment with change in temperature ; shaking of an instrument 
in the wind, etc. But most important of all are those errors 
which arise from imperfections in the sight, hearing, and other 
senses of the observer, which render it impossible for him to 
adjust and use his instruments with absolute accuracy. 

After a full investigation of the constant errors, the observer 
should diminish the accidental errors as much as possible, both 
in number and magnitude, by taking every precaution and care 
in the measurements themselves. The problem now remains 
to combine the observations so that the remaining accidental 
errors shall have the least probable effect upon the results, 
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and it is to bring about this combination of observations that 
we employ the Method of Least Squares. 

When no more observations are made than are sufficient to 
determine one value for each of the unknown quantities, we 
must accept these values as the most probable ones. But if 
additional observations are made leading to discordant results, 
we can not take any one of them as the correct value, and in 
fact, as already stated, we shall probably not be able to obtain 
the true values of the unknowns. All that we can do is to find 
values of the unknowns which shall remove the discrepancies 
between the different observations and which shall be those 
values that are rendered most probable by the existence of 
the observations themselves. 

On first thoughts it may seem that these accidental errors, 
being due to so many different and unknown causes, will be 
beyond the scope of mathematical investigation. Neverthe- 
less, the theory of probability requires that these errors shall 
follow in magnitude and frequency a law that is capable of 
exact mathematical expression, and experience confirms the 
correctness of this law. 

For more extended remarks on these subjects see — 

Holman, "Discussion of the Precision of Measurements," pp. 1-14. 
Merriman, " Text-Book of Least Squares," pp. 1-6. 
Chauvenet, " Spherical and Practical Astronomy," pp. 469-473. 
Wright, "Treatise on the Adjustment of Observations," pp. 11-18. 



LAWS OF ERRORS OF OBSERVATION. 

7. The derivation of the general laws of the occurrence of 
errors of observation, and of the processes for determining 
the most probable values of the unknown quantities, will be 
based upon the following 

Axiom. If a series of n direct observations, M^, M^, . . . M„, 
is made upon the value of a quantity M, all the observations 
being made with the same care and under the same circum- 
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stances, the most probable value M^ of that quantity is the 
arithmetical mean of the observations. Or 

Mo- - -^^ (I) 

8. The Real Error (a;) of an observation is the difference 
between the observed value of the measured quantity and the 
real value. 

9. The Residual (v) of an observation is the difference 
between the observed value of the measured quantity and 
the value rendered most probable by the existence of the 
observations. 

10. Exanhple. Eight observations are made upon the 
resistance of a coil of wire, the true resistance being 512. 
Find from these observations the most probable resistance, 
and also the real errors and residuals. 



Observations. 


Real Errors. 


Besiduals. 


M. 


X 


V 


512.4 


+ •4 


+ .30 


512.2 


+ •2 


+ .10 


511.9 


-.1 


- .20 


512.3 


+ .3 


+ .20 


511.8 


- .2 


-.30 


512.3 


+ .3 


+ .20 


511.9 


-.1 


-.20 


512.0 


.0 


-.10 


Mean = 512.10 


%v=. .00 



From the observations, then, we should say that the most 
probable resistance of the coil is 512.10. It wUl also be 
noticed that the sum of the residuals is zero. That this is a 
general result following from the assumption of the arith- 
metical mean as the most probable value may be proved as 
follows : If the observations are Jl^, M^, . . . M„, the 
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arithmetical mean M,,, and thte residuals v^, v^, . . . «„, then 
•we have 

%v = %M— nMo 

'%M 

^%M—%M, since iyi= 

n 

.: 2;v = (2) 

11. Weigfhted Observations. The weight of an observa- 
tion expresses its relative worth compared with other obser- 
vations. Thus, if six observations are made upon the value 
of a quantity, five of which give the same result, while the 
sixth differs, in combining these two different results to 
obtain the most probable value of the unknown, the first 
value ought to have five times the influence upon the final 
result that the second has, since it has taken five times as 
much labor and time to obtain it. Hence in general we may- 
say— 

12. The Weight {p) of an observation may be considered 
as representing the number of times the observation has been 
repeated and the same result obtained. 

The weights assigned to observations may be due to a 
variety of causes, as difference in skill of observers, difference 
in the instruments used or the circumstances under which 
the observations are made, etc. But whatever the cause, the 
effect on the final values of weighting an observation will be 
the same as indicated in the preceding paragraph. 

13. Example. Suppose n observations, i!/[, M.^, . . . M^, 

of weights pi, p Pn, are made upon the value of a 

quantity M. To find the most probable value M^ of the 
quantity. 

From the above interpretation of the meaning of weight, 
we may consider that the whole number of observations is 
Pi~\-Pi-\- • • -Pm or Si>, and that the result Mt has been 
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obtained in p^ observations, M^ in p^ observations, etc., 
Therefore, by (1) 

^»- ^ -~s^ '■^^ 

Ma is called the General Mean. 
If the residuals are Wi, v^, . . . «„, we have 
Vi = i!fi — Jfo, ?)2 = ili^ — Jfo, • • • «»= -31^ — -3fo, 

"Lpv—Q (4) 

Which shows that in the case of direct observations of differ- 
ent weights the sum of the weighted residuals is zero. 

14. If the observations are not made directly upon the 
values of the required quantities, the method of adjusting 
the results so as to obtain the best possible values of the 
unknowns will depend upon the laws which govern the dis- 
tribution of the errors of these observations. It is found in 
practice that the accidental errors of observations follow cer- 
tain well defined laws, and what these are may best be seen 
by taking an actual example. 

15. Example. One thousand shots are fired at a target 
which is divided into a number of horizontal sections by lines 
one foot apart, the centre line of the target being in the 
middle of one of these spaces. The shots were distributed 
as follows : — 



In Space. Shots. 


In Space. 


Shots. 


In Space. 


Shoit. 


-f5ito+4i 1 


+Hto+i 


190 


-2^to-3i 


79 


-f4i«+3i 4 


+ * " - i 


212 


-H " -^ 


16 


+ 3^ " +^ 10 


- i « -H 


204 


_4^ « _5J 


2 


-1-2^ « +4 89 


-1* « -^ 


193 







In this case the errors are evidently the distances of the 
shots from the centre of the target. Further, as far as can 
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— i and — 1^ is 


.204 


-li « -2J « 


.193 


-2^ « -3J " 


.079 


-3i " -4^ « 


.016 


_4i « -bi « 


.002 



be judged from these one thousand shots, if another shot is 

fired the probability that this shot will fall between the 

lines 

+ 5^ and +4^ is .001 

_|_4 " +4 " -004 

4-3^ « +2^ " .010 

-\-^ « -j-li " -089 

+ 1* " + i" -190 

+ i " - i " .212 

The sum of the above probabilities is unity, and, therefore, 
as far as the preceding shots show the 1001st shot will cer- 
tainly hit the target. 

16. Now using as abscissas the distances of the horizontal 
lines from the centre of the target, and as ordinates the num- 
ber of shots falling in the corresponding spaces, we may 
construct the following figure : — 



Figure 1. 



And if the entire area of this figure is taken as unity then 
the area of each rectangle will denote the probability of a 
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9 



shot, if fired, falling within the corresponding space of the 
target. 

The graphical representation of the accidental errors of 
observation will always give a figure similar to the above. 
Hence denoting errors by abscissas, and their frequency by 
ordinates, the law of error of any series of observations may 
be represented by a curve whose general form is determined 
by Figure 1. This curve is called the " Curve of Error," and 
is shown in Figure 2. 

r 



Figure 2. 




In order that this curve may represent exactly the distri- 
bution of the errors in any given series of observations it 
ought to meet the axis of JC at some definite distance to the 
right and left of the origin and coincide with the axis from 
there on, for in all actual observations there is a limit beyond 
which no errors occur. But as the exact point of meeting 
could not be determined for any given case, and as it would 
not be possible to obtain the equation of such a curve, we 
make it asymptotic to the axis of X, taking care that the 
error thus introduced shall in any set of observations be so 
small as to be negligible. 
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17. An inspection of Figures 1 and 2 will now exhibit 
some of the general lawe of errors of observation and the 
corresponding properties of the curve of error. 



Laws of Error derived from an 
inspection of Figure 1. 



Sepresentation of these laws by the 
Curve of Error. 



First. Small errors are more 
frequent than large ones. 

Second. Positive and negative 
errors of the same absolute mag- 
nitude are equally likely to occur. 

Third. The probability of the 
occurrence of very large errors is 
very small. 

Fourth. The frequency of any 
error depends upon the magnitude 
of that error. 



The maximum point of the 
curve is on the axis of T. 

The curve is symmetrical with 
respect to the axis of T. 

The curve is asymptotic to the 
axis of X. 

The equation of the curve will 
be of the form 



y = ^{x) 



(5) 



18. If, now, the total area between the curve and the axis 
of X" be denoted by unity the probability that the error of 
any given observation will fall between the magnitudes x and 
x-\-dx will be represented by the area included between the 
curve, the axis of JE^ and the ordinates of the curve at the 
errors x and x -f- dx, or by 



y dx=<f>(x) dx 



(6) 



And this probability will be known as soon as we find the 
form of the function <^(a!). 

19. The above expressions in (5) and (6) are the ones 
that we should use if we regard the curve of error as repre- 
senting the law of occurrence of errors of observation. If, 
however, we look upon the curve as expressing the law to 
which we must make the residuals conform, in order that the 
values of the unknown quantities obtained from them may be 
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the most probable values, we should replace a; by w and use 
the expresBions 

y=4>{^) (7) 

and ydv=i<^{v)dv (8) 

for (5) and (6), respectively. 

THE EQUATION OF THE CURVE OF ERROR. 

20. Let n observations, all of the same weight, with 
results Ml, M^, . . . J^, be made upon any function or 
functions of a number of unknown quantities Zi, Zj, . . . z^ ; 
and let the residuals of J^, M, . . . M^ be Wi, Wj, . . . w„, 
and the probability of the occurrence of these residuals be 
^(«i) dv, <j}(v2) dv, . . . <l>{v„) dv, respectively. Then the 
probability of the simultaneous occurrence of all these 
residuals will be 

P=.^(«i) <t>(v,) . . . ^(v„) {dvY (9) 

.-. \ogP — \og4,{v{)-\-...\og4>{Vn)+n\ogdv (10) 

Each different method that might be adopted for computing 
the values of the unknowns Zi, Zj, . . . Zg would lead to a dif- 
ferent set of residuals Vi, v^, . . . v^; but obviously that set 
of values of Zi, Zj, . . . Zg should be considered the best which 
corresponds to the particular set of residuals «i, v^, . . . v^, 
the probability of whose occurrence is greater than that of 
any other set. 

Therefore the most probable values of Zi, Sj, . . . z^ are 
those that make J* in (9), or log P in (10), a maximum. 

The values of Zj, Zg, . . . Zg corresponding with this latter 
condition are those that satisfy equations (11). It maybe 
noticed that these equations also express the preliminary 
conditions leading to a minimum value of log P, but the 



12 METHOD OF LEAST 8QUABES. 

nature of the problem is such that a maximum value of P 
evidently exists while a minimum does not, and it is there- 
fore unnecessary to investigate further the mathematical 
conditions for a maximum. Hence we have 

aiogP ^ 1 d^(v,) 1 d <l>(v„) ^ ^ 






(11) 



3iog^^_J_a^0O , 1 a <i>(vn) ^ Q 

SSj ^(«i) dZg "T • ■ • ,^(«„) azg 

and for convenience we may put 

<I)'(V) 

|^ = 'A(«) (13) 

substituting in (11) we have 

-A(Ma^+^(M5^^+-..^(«»)a^"=o 



^(M5z^ + ^(-.)a^^+...^(M5z=0 



(14) 



These equations contain all the unknowns Zj, z^, . . . Zq, and 
there are as many equations as unknowns, hence as soon as 
we find the form of the function i/'(v) we can solve these 
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equations for the most probable values of Zj, %, . . . Zg. Since 
we have considered the general case, and the above results 
are to hold true whatever the number of unknown quantities 
and the form of the functions observed, we may deduce the 
form of ^{v) by solving a special example. 

Example. Let n observations of equal weight be made 
upon the value of a single unknown %, with results M^, M^, 
. . . Mm and let the residuals be v-i, v^, . . . v^. Then the 
most probable value of Zi is given by 

Zi=: Ml — Vi=z Mi — V2 = . . . Mn — Vn 

differentiating with respect to Si, 

1— _— ' — _^^— _^ /\ 

9zi dzi ' ' ' dzi ^ ' 

substituting (a) in (14), changing all the signs, we have 

f(vi)+f(v,)-\-...^(v„)=.0 (b) 

But in this case, as was shown in (2), 

«i + «2 + ••■«» = (c) 

In order that (b) and (c) may both be true the functional 
symbol tj/ must indicate multiplication by a constant. That 
is, in general 

,p(v)z=cv (15) 

Substituting this in (13) and (12), 

dd>(v) dv 



therefore 


1 d<t>(v) 
<l>{v) dz 


dv 

= CV •=;— 

dz 


Integrating, 


log ^(«) : 


= \cv^^d 

— giCT' + C' 
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Since y=<^(w) is the equation of the curve of error, (7), 
we may therefore write it 



y- 



ke^'^ 



But on examination of the curve, y is seen to be a decreas- 
ing function of «, and hence the exponent of e is essentially 
negative. Accordingly we will write our equation in the 
form 

y = he-^'' (16) 

the values of h and A depending upon the character of the 
observations, but in all observations of the same kind and 
weight having the same values. 

This equation represents the law in accordance with which 
the residuals must be distributed in order that the best results 
may be obtained from our observations. But, as before 
mentioned, if we wish our curve to represent the most prob- 
able distribution of the real errors of observation we should 
write the equation in the form 

y=ke-K^' (17) 

Hereafter we shall use without further remark either form 
of the equation according to the aspect in which we are 
considering our curve. 

An inspection of the above equation will show that it 
satisfies all the conditions noted in discussing the form of the 
curve of error in paragraph 17. 

21. It is important to notice that in all discussions in the 
Method of Least Squares the number of observations is sup- 
posed to be large and always greater than the number of 
unknown quantities. As will be illustrated later on, para- 
graph 91, whenever this is the case there is a remarkable 
agreement between the results obtained in practice and those 
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indicated by the theory. And even when the observations 
are few in number the method still affords the best means at 
our command for their adjustment, the results obtained 
merely having a smaller weight than they would have had if 
derived from a greater number of observations. 

THE METHOD OF LEAST SQUARES. 

22. We are now in a position to see whence comes the 
name "Least Squares." 

In paragraph 20 it was pointed out that whenever we 
make a series of observations, each observation of the set 
having the same weight, the most probable system of values 
of the unknown quantities will be that which corresponds 
with the set of residuals the probability of whose occurrence 
is a maximimi. That is, the best set of values of the 
unknowns will be that which gives a maximum value to 

P==^(«i) <^(«,) ...4>{^n) {dvT 

But from equation (16) this reduces to 

P = A™e-'^'(''i' + "2=' + • • • V) (<?„)» (18) 

Since the exponent of e in this expression is negative, 
evidently P will be a maximum when 

Vx-\-Vi-\- ■ . ■ Vn=^v'^ is a minimum. (19) 

Hence the adjustment of observations by the Method of 
Least Squares is based upon the principle that the most prob- 
able system of values of the unknowns is that which renders 
the sum of the squares of the residuals a minimunx. Hence 
the name. 

The conditions for a maximum value of P were expressed 
in equations (14), and since it has been shown that the 
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function ^ means multiplication by a constant, those equa- ' 
tions reduce to the following, called 

NORMAL EQUATIONS. 






(20) 



An inspection of these equations will show that they also 
express the conditions that will make the sum of the squares 
of the residuals a minimum. 

In the adjustment of observations the above are the funda- 
mental equations. In order to obtain the most probable 
values of the unknowns in any set of observations, all that is 
necessary is to form the Normal Equations for that set and 
solve them simultaneously. The examples already solved for 
direct observations are merely special cases of the above 
general solution. 



CHAPTEE II. 

THE ADJUSTMENT OF OBSERVATIONS. 

INDIRECT OBSERVATIONS. 

23. In the determination of the values of quantities by 
means of observations the functions of the unknowns that it 
is necessary to observe may be of any form, but if they are 
not linear the normal equations derived from them are likely 
to be complicated and difficult, if not impossible, to solve. 
Hence if the observations are not upon linear functions of 
the unknowns, the first step will be to reduce them to 
equivalent linear expressions by transformations depending 
on the character of the observed functions ; see paragraphs 
43 and 44. It will be necessary to consider, therefore, the 
method of adjusting observations on linear functions alone, 
and the procedure in cases of this kind may be illustrated by 
the following simple example. 

24. Example. Si, iS^, S^, are three solids whose masses 
are required. Not having standard weights enough to obtain 
all these masses directly, by varying the distribution of the 
solids in the pans of the balance the following observations 
are made : — 

*?i = /Sa -|- 1.7 grams. 
S, — 2.4 " 



"3 == -'•'* 



If, now, the most probable values of Si, S^, S^, are repre- 
sented by Si, Zj, Zs, and the corresponding residuals of the 
observations by Vi, v^, . . . Vi, the values of these residuals in 
terms of Zi, z^, Zj may be found from the above observations 
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by transposing all terms to the first members of the equa- 
tions, and we obtain at once the following, called 

OBSERVATION EQUATIONS. 

gj — Zj — 1.7 = Vi 



Z3 — 2.4 = V2 



(A) 



— Si -f Zj 4- Z3 — 1.0 = Vi 
^2 — ^s — 3.0 = «4 

Applying equations (20) as formulas we have the 

NORMAL EQUATIONS. 
(% - S2 - 1.7) + (- ^1 + Z2 + 2s - 1-0) (- 1) = 

(Si - S2 - 1.7)(- 1) + (- Z1 + Z2+ Za - 1-0) 

+ (z2-Z8-3.0) = (B) 

(Z3 - 2.4) + (_ Z1 + Z2 + Z8 - 1.0) 

+ (Z2-Z3-3.0)(-l) = 

Simplified, these become 







2zi - 2Z2 


— 


88 


-0.7=0 


(a) 






— 2zi + 3z2 






-2.3=0 


(b) 






— % 


+ 


3Z8- 


-0.4=0 


(c) 


3X(a) 




6zi — 6Z2 


— 


3Z3 


-2.1 = 




(c) 




— 2l 


+ 


3z3 


-0.4=0 






5zi — 622 






— 2.5=0 




2X(b) 




-4zi + 6z2 






— 4.6=0 






2i 






-7.1 = 














zi=7.1 


(d) 


substitute 


(d) 


in(b) 






Zj=5.5 


(e) 


(1 


(d) 


in (c) 






Z3=2.5 


(f) 



An inspection of the work in this example will show that 
for the adjustment of observations of equal weight on linear 
functions of the unknowns we may derive the following : — 
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25. Rule. For each observation write an " Observation 
Equation " / then for each unknown form a " Normal 
Equation" by multiplying the first member of each observa- 
tion equation by the coefficient of that unknown in that 
equation, adding the results and placing the sum equal to 
zero. Solve these equations simultaneously for the values 
of the unknowns. 

In solving for the most probable values of the unknowns 
the second members of the observation equations are very 
commonly written zero instead of «i, v^, . . . v^. For this is 
the form in which the equations naturally appear, and if the 
observations were exact the residuals would actually all be 
zero. The method of solution is the same in either case. 

26. Observations of Unequal Weight. If the observations 
are not all of equal weight the same method will apply, 
except that in the formation of the normal equations each 
observation equation will be used the number of times denoted 
by its weight. Thus in the last example if the observations 
have the weights 4, 9, 1, 4, the normal equations will have 
the same form as in (B), page 18, but each part of each 
equation will be multiplied by the weight of the observation 
equation from which it is derived. This will give the 

NORMAL EQUATIONS. 
4(2i-Z,-1.7) + (-Zi + 2, + 23-1.0)(-l)=0 

4(zi - z,- 1.7)(- 1) + (- zi + s, + z, _ 1.0) 

+ 4(^2 -23- 3.0) = 
9(z3-2.4) + (-Zi + z, + Z3-1.0) 

+ 4(s2-S3-3.0)(-l) = 
or reduced 

5zi — 5z2 — Zs — 5.8 ^ (a) 

— 5zi+9z2— 3zs— 6.2=0 (b) 

_ zi _ 332 + 14«3— 10.6=0 (c) 

the solution of which gives 

Zi=7.07 2,= 5.42 z,= 2.42 
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Further it will at once be seen that if Pi, p^, . . . Pm ^^^ ^^^ 
weights of the corresponding observations, equations (20) 
take the general form: — 

WEIGHTED NORMAL EQUATIONS, 



" " — — _. „ 

(21) 



i>i «^i a^ +i>. -y^ a^ + . . ■ i>« -«„ 5^ = 






p^'^^di+^^'^'di + ■ ■ ■^-''-Bi=^ 



Hence for the formation of the normal equations in weighted 
observations on linear functions of the unknowns, we have 
the following: — 

27. Rule. J''or each observation write an " Observation 
Equation " / then for each unknown form a '•'•Normal 
Equation^'' by multiplying the first member of each observa- 
tion equation by the coefficient of that unknown in that 
equation and by the weight of that equation, adding the 
results and placing the sum equal to zero. Solve these 
equations simultaneously. 

28. The same result will be obtained if we begin by 
multiplying each observation equation by the square root of 
its weight and then proceed according to the first rule 
(paragraph 25). 

This result illustrates the important principle that multi- 
plying a set of equations by the square roots of their weights 
reduces them all to equivalent equations of weight unity. 

29. Belation between the Weight of an Observation and 
the Value of Tl, If in paragraph 20 the n observations 
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have weights pi, p^, . . . p„, and the quantity A values 
hi, h^, . . . hn, then equation (18) becomes 

P = Aie-^i'«i' k^e-f'^'"'^^ . . . ;5;„e-'^"'''»' (c?u)" 

= hK ■ ■ ■ *«e-(fti' "i' + "ii' "2^ + • • • W*™'^) (^„)» (22) 

The most probable set of values of the unknowns is that 
which makes P a maximum, and i^ is a maximum when 

hiVi^-{- Jit^Vi^-\- . . . hnVn is tt minimum. (23) 

The conditions for a minimum value of this expression are 
the following, which are then for this ease the 

NORMAL EQUATIONS. 






(24) 






But equations (21) are also the normal equations for 
this case. Hence (21) and (24) must be identical, and 

Pt. : Pi : . ■ ■ Pn= Til : h,^ : ■ ■ ■ Tin (25) 

That is, the square of h is proportional to the weight of 
the observation. Accordingly, since A increases in value as 
the quality of the observations is improved, it is called " The 
Measure of Precision." 
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Further, it follows from (23) and (25) that the most 
probable system of values of the unknowns will be that in 
which 

PxVx^-\-p2Vi-\- . . .pnVn is a minimum. (26) 

And this is the most general form of statement of the 
principle of " Least Squares." The same principle is repre- 
sented in equations (21). 

30. Computation of Corrections. If large numbers occur 
in the observations it is better to compute the most probable 
corrections to apply to the observed values rather than the 
most probable values of the unknowns themselves. In this 
way we can often avoid a large amount of numerical work. 

31. Example. P^, P^, P^, P^, P^, are five points whose 
altitudes above the mean level of the sea are to be determined 
from the following observations of difference of level. 

/>! = 573.08 P4- Pa =170.28 

Pj-Pi= 2.60 ^4—^5 = 425.00 

^2 = 575.27 A = 319.91 

^'3 — ^2 = 167.33 P5 = 319.75 

P^-Ps= 3.80 

An inspection of these observations shows that we may put 

Pi = 573 + Si ^4=745 + 24 

P2 = 575 + Sj Pj = 320-t-S5 (A) 

P3=742 + 2, 

where Zj, z^, Sg, 84, S5, are small corrections whose most 
probable values are to be determined. We now have for 

OBSERVATION EQUATIONS. 

573 4. sj _ 573.08 = or 2i-.08 = 

575 + 02-573 — 21- 2.60=0 or 2j- 21 -.60 = 

575 + 2^-575.27 = or 2^ — .27 = 



or 


Z3-Z2-.33=0 


or 


Z4 — Zs — .80 = 


or 


Z4-Z2-. 28 = 


or 


^4 — ^6 =0 


or 


z,+ .09=0 


or 


Z5+.25=0 
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742 + Zs - 575 - Zj - 167.33 = 
745 -I- z^ — 742 - Zs — 3.80=0 
745 4- ^4 — 575 - zj - 170.28 = 
745 -|- zj — 320 - Z5 - 425.00 = 

320 + Z5- 319.91 = 

320 -I-Z5- 319.75 = 

From these we now form the 

NORMAL EQUATIONS. 

2zi— Zj -|- .52 = 

— Zi -|- 4z2 — Z3 — Z4 — .26 ^ 

— Z2 + 2Z3- Z4 + .47=0 

— za- Zs + 3z4- Zj— 1.08 = 

- Z4 + 3Z54- .34 = 
and solving, 

Zi= — .19; Z2 = .14; Z8=.05; Z4=.43; Z5 = .03 

Substituting these in equations (A) we have for the most 
probable altitudes, 

Pi = 572.81 Ps= 742.05 ^5 = 320.03 
P2 — 575.14 P4 = 745.43 

If the original observation equations had been retained, 
the independent terms in the normal equations would have 
been 

570.48 240.26 163.53 599.08 214.66 

32. Significant Figures. The adjustment by the Method 
of Least Squares of observations which occur in practice, 
although not difficult, is apt to be long and laborious. Hence 
to reduce this labor as much as possible it is of great import- 
ance that careful attention should be given in the solutions 
to the proper use of significant figures. When in doubt, 
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however, as to the proper numher of figures to retain it is 
better to keep too many rather than too few, as the superfluous 
figures can be rejected at the end of the computation ; while 
if too few are retained the results obtained from the compu- 
tations will be worthless. 

For a general discussion of the subject of significant figures 
see Holman's " Precision of Measurements," pages 76 to 84, 
but for the present the following rules will suffice for most 
cases. 

Bule 1. In casting off places of figures increase by 1 the 
last figure retained, when the following figure is 5 or over. 

Rule 2. In the precision measure retain two significant 
figures. 

Rule 3. In any quantity retain enough significant figures 
to include the place in which the second significant figure 
of its precision measure occurs. 

Bule 4. When several quantities are to be added or 
subtracted, apply Mule 3 to the least precise and keep only 
the corresponding figures in the other quantities. 

Rule 5. When several quantities are to be multiplied or 
divided into each other, find the percentage precision of the 
least precise. If this is 

1 per cent or more, use four significant figures. 
.1 " " « " " five " " 

.01 " " " " " six " " 

in all the work. If the final result obtained in this way 
conflicts with Rule 3, apply the latter. 

Rule 6. When logarithms are used, retain as many 
places in the mantissce as there are significant figures 
retained in the data under Rule 6. 

The application of these rules is not always possible in the 
course of the work, since the precision measures may not be 
known until the end of the computation. But as a general 
rule it is sufficient in direct observations to retain one more 
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place of figures than is given by the individual observations, 
and in indirect observations to retain two additional places. 

CONDITIONED OBSERVATIONS. 

33. Conditioned Observations are those in which the 
unknown quantities must be determined not only so as to 
satisfy as closely as possible the observation equations, but 
also so as to satisfy exactly certain other conditions. These 
conditions must be less in number than the unknown quantities, 
otherwise the unknowns could be determined from the con- 
ditions alone. 

The adjustment of observations of this class may be reduced 
to the method already used for unconditioned observations in 
the following way. 

The observations are represented by " Observation Equa- 
tions," and the conditions by certain other equations, called 
" Condition Equations." 

Between these two sets of equations we will eliminate as 
many unknowns as there are conditions. From the resulting 
equations, which will be the same in number as the observa- 
tions, we will form in the usual manner the "Normal 
Equations " for the remaining unknowns. Having solved 
these normal equations and substituted the results in the con- 
dition equations, we shall obtain the values of the unknowns 
first eliminated. 

All conditions of the problem are now fulfilled, for the 
condition equations are satisfied exactly and, moreover, 
according to the principle of Least Squares our results 
are those rendered most probable by the existence of the 
observations. 

As in the example last considered, it is often more 
advantageous to compute corrections to the observed values 
of the unknown quantities rather than the values of the 
quantities themselves. 
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34. Example. Find the most probable values of the angles 
of a quadrilateral from the observations, 



^=101° 13' 22" weight 3 

B= 93 49 17 « 2 

G= 87 5 39 " 2 

D= 77 52 40 " 1 



(A) 



Jl_|_^_|-(7_|_i)=:360° 0' 58" 

The condition to be satisfied is in this problem 

^ + ^ + C+i>=360'' (B) 

Let Si, ^2, 8s) ^4 be the most probable corrections to add 
to the observed values. This gives for 

OBSERVATION EQUATIONS 



Si ^ weight 3 

S2=0 " 2 

S8=0 " 2 

»4=0 " 1 



(C) 



and for the 



CONDITION EQUATION 

Zi + S2+gs + Z4+58 = (D) 

Eliminating »4 between (D) and (C), the equations from 
which the normal equations are to be derived become 

% = weight 3 

z.= « 2 (^^ 

Si + »2 + »3+58=0 " 1 
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Applying the rule in paragraph 27, these give the 

NORMAL EQUATIONS 

4si-)- 2^+ S3 + 58 = 
Si + 3z,+ 23 + 58 = (F) 

2i+ % + 3ss-|-58 = 

Solving, and substituting the results in equation (D), we 
find 

Zi=- 8.29 Zs=- 12.43 ,g,. 

82= -12.43 Z4=- 24.85 ^ ' 

Applying these corrections to the observations (A), the 
most probable values of the angles are 



A = 101° 18' 13".71 

^=93 49 4.57 

(7= 87 5 26.57 

I>= 77 52 15.15 



(H) 



Note. In eliminating unknowns between the observation 
and condition equations care must be taken that the obser- 
vation equations are not combined with each other or multi- 
plied by any quantity. For if this is done the weights of 
the observation equations will be altered. (See §28.) 

35. In the above example it is evident that the corrections 
to be applied to the different observations are inversely as 
their weights. And, in general, when there is but one 
equation of condition, the observations expressing direct 
determinations of the unknowns, the corrections will be pro- 
portional to the coefBcients of the unknowns in the equation 
of condition divided by the weights of the corresponding 
observations. A proof of this is given in paragraph 117. 
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The most common case is that in which these coefficients 
are all unity, as in the example just solved, and we may then 
derive the 

Bule. Find the difference between the theoretical and 
observed results and divide this correction among the 
observations in the inverse ratio of their weights. 

In the last example the sum of the observed angles exceeds 
360° by 58". Therefore the correction to be applied to A is 

1 
- 58 X --^- = - 58 X y = - 8.29 

EMPIRICAL FORMULAS AND CONSTANTS. 

36. In the work so far considered the observations are 
supposed to be made either directly upon the values of the 
unknown quantities or upon some function of the unknowns 
whose form, and the constants entering into it, are definitely 
known. But another sort of problem frequently occurs, 
in which observations are made upon the values of a certain 
variable and the corresponding values of some function of it, 
the exact form of the function not being known. The object 
in this case is the determination of the most probable form of 
the function and the values of the constants involved ; that 
is, the derivation of the algebraic expression best representing 
the law connecting the variable and function. 

This expression may be looked upon as the equation of a 
curve, abscissas denoting values of "the variable and ordinates 
values of the function, and for all values of the variable 
within the range of the observations we may determine from 
it the most probable values of the function corresponding. 
But except in special cases, where the number of observations 
is large, where the law connecting variable and function is 
well defined, and where the equation obtained is an accurate 
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representation of this law, it cannot be assumed to apply 
beyond the range of the observations. And in no case would 
it be safe to make use of the curve very far beyond the limits 
of the observations. 

37. The Method of Least Squares will not assist in deter- 
mining the form of the function. This must be settled upon 
beforehand, either from theoretical considerations or by 
constructing a plot, using values of the variable as abscissas 
and of the function as ordinates, when the smooth curve 
drawn through the points thus obtained wUl indicate the 
form of equation to be used. 

It is to be observed that this is a method of trial and will 
not necessarily give the most probable form of the function ; 
and in fact we may not be able to obtain the form that would 
be absolutely best. Further, several forms of equation may 
be known which would represent well the plotted points. In 
such a case that should be considered the best in which the 
sum of the squares of the residuals is found to be the least. 

38. As soon as the form of the function is decided upon it 
should be reduced to the linear form, and the determination 
of the values of the constants involved is then a simple 
application of the preceding methods. 

As the " Observation Equations " in any given problem will 
all be of the same kind, it is usually advisable to write out 
the general form of the "Normal Equations" and arrange 
the computations in tabular form, while the retention of 
the proper number of significant figures is of particular 
importance in this work. 

39. A case that frequently occurs is that in which the 
quantity y is a constantly increasing function of the vari- 
able X, or where the plotted curve is approximately parabolic 
in form. Here the equation 

y = A + B3c+Cx^+Dx^+... (27) 

may be taken to represent the relation between the variable 
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and the function. The larger the number of terms taken in 
the second member, the more accurately may the equation 
obtained be made to represent the results of the observations ; 
but the labor involved increases rapidly with increase in the 
number of terms, and if the plot shows a very nearly straight 
line the first two terms alone may suffice. 

40. Example. In measuring the velocity of the current 
of a river the following results were obtained : — 



Depths. 


Velocities. 


X 


V 


1 


4.86 


2 


5.14 


3 


5.15 


4 


4.85 


5 


4.24 


6 


3.36 


7 


2.16 


8 


0.67 



The velocity at the surface is 4.250. Find the equation of 
a curve which will express the relation between x and V. 
Plotting the observations we find the curve 

V 
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This is approximately parabolic in form and passes through 
the fixed point (0, 4.25). Therefore the relation between x 
and V vaa,j be expressed by the equation 



F=4.25 + ^a!+ Gx" 



(A) 



and substituting in this the corresponding values of x and V 
as given by the observations, we shall have eight observation 
equations from which the most probable values of B and O 
are to be computed. 

All of the observation equations being of the form (A) we 
have the 

NORMAL EQUATIONS 

C%x''-\-B%x'-\-4:.1b%x^--S,Vx^=(i (a) 

C%x^-{-B S£b''+ 4.25 %x — S Fa; = (b) 

For computing the coefficients in these equations it is most 
convenient to arrange the following table. 



X 


V 


Vx 


x^ 


Vx^ 


x' 


X* 


1 


4.86 


4.86 


1 


4.86 


1 


1 


2 


5.14 


10.28 


4 


20.56 


8 


16 


3 


5.15 


15.45 


9 


46.35 


27 


81 


4 


4.85 


19.40 


16 


77.60 


64 


256 


5 


4.24 


21.20 


25 


106.00 


126 


625 


6 


3.36 


20.16 


36 


120.96 


216 


1296 


7 


2.16 


15.12 


49 


105.84 


343 


2401 


8 


0.67 


5.36 


64 


42.88 


512 


4096 


36 




111.83 


204 


525.05 


1296 


8772 


%x 




%Vx 


■S,x' 


%rx^ 


2iB» 


%x* 
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Substituting these results in (a) and (b) we have 

8772C+1296^ + 341.95 = (c) 

1296 C+ 204i?+ 41.17 = (d) 

and solving, 

C=-.1493 ^=.7465 (e) 

Therefore the required equation is 

V = 4.25 + .7465a! - .1493a!'' (f ) 

Whenever the quantities in the observations are so large 
that the use of logarithms is desirable, these can best be put 
in the same columns directly over the natural numbers. 

41. It is to be remarked that in solutions like the above 
we assume, from our method of forming the observation and 
normal equations, that the observations on the values of the 
function are alone subject to error, the observations on the 
variable being supposed to be exact or to have errors so small 
as to be negligible. 

42. Periodic Phenomena. If as the variable increases the 
function passes through recurring values, that is, if y is a 
periodic function of a!, some form of trigonometric equation 
would be the proper one to select. For instance, a good form 
to use is 

y = A + B8in ^x + C cos ^ac (28) 

where A, £, C are the constants whose most probable values 
are to be found, and m is the number of units of x comprised 
in the entire cycle of values of y. The quantity m is to be 
determined from an inspection of the observations, and if it 
appears that several values of m might be used, all should be 
tried, and that which leads to the smallest sum for the squares 
of the residuals is to be considered the best. If the several 
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cycles are not similar and regular, additional terms involving 
multiples of x will have to be added to equation (28). 

43. Special I Treatment of Exponential Equations. If the 
equation selected is not in the linear form as regards the 
unknown constants the general method of procedure is 
given in paragraph 44, but in some special cases a more 
simple reduction is possible. A quite common case is the 
following: — 

Suppose the relation between x and y is expressed by the 
equation 

y = &«" (29) 

the problem being to obtain the best values of h and m. 
Taking the logarithms of both members of (29), 

log y = log k-\-7n log x 

Denoting log Jc by Jc', this equation becomes 

m log x-\-k' — log y = 

which is in the linear form as regards the unknowns m and 
k', and the normal equations may now be formed in the usual 
manner. 

The most convenient way of determining whether equation 
(29) is a suitable one to select or not is to plot the corre- 
sponding values of x and y on logarithmic cross-section 
paper. If (29) is a proper equation to use the plotted points 
will lie upon a straight line. 

REDUCTION OF OBSERVATION EQUATIONS TO THE 
LINEAR FORM. 

44. When the observation equations are not linear as 
regards the unknowns, the only practicable method of pro- 
cedure, as already mentioned in paragraph 28, is to reduce 
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them to that form. This reduction may be efEected as 
follows : — 

Let the observation equations be 

f,(Z„ Z„... Zg) = M^ 

fi(Zi, Zi, . . . Zg) = M^ 

(A) 



/„(Z. Z„... Zg) = M„ 

in which Z^, Z^, . . . Zg represent the unknown quantities and 

M^, M^, . . . Mn the observations, the functions being of 
known form. 

Let Zi, Z^', . . . Zg be approximate values of Z-^, Z^, . . . Zg, 
found by trial or the solution of a sufficient number of the 
observation equations, and let the most probable values of 
Zi, Z2, . . . Zg be 

Z^'+ZU Z,'-\-Z„...Zg'-\-Zg (B) 

Zj, S2, . . . Zg being small corrections whose values are to be 
determined by the Method of Least Squares. The first 
observation equation in (A) may then be written, 

A{Z,'-\-Z„ Z,'-\-Z„ . . . Zg'-\-Zg) = M, 

Expanding the first member by Taylor's Theorem, denoting 
/x(Zi', Z2, . . . Zg) by ^1, and neglecting terms containing 
powers of Zi, z^, . . . Zg higher than the first, this becomes 

I. ± ^h. _L i^ _L ^^1 _ 31*- 

*^ ■+■ 5Z^ ^' "I" dZ^ 82 -t- • • • Qz^ 2s — Ml 

If now we represent the coefficients of z^, Sj, . . . Zg, by 
fflj, Jj, . . . g-i, also ki — Ml by wii, and treat the other 
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equations in (A) in the same manner, our observation equa- 
tions will take the form 

OiZi + S1Z2 + . . . qxZq -\- mi = Vi 

a^Zi -\- 62% + . . . q2«q + »h = «2 

(C) 

The second members reduce to the residuals since the 
quantities %, Zj, . . . Zg represent merely the most probable 
values of the corrections. Equations (C) can now be solved 
in the usual manner and the most probable values of 
Zi, Z^, . . . Zg found by substituting the results in (B). 

45. Example. From the following observations find the 
most probable values of x and y. 

sin X -\- cos 2y = 1.5 
cos a; -j- 3 sin y = 1.7 (A) 

x'-\-by = 2.1 

By trial it is found that 42° and 18° are approximate 
values of x and y. Then in accordance with paragraph 44 
we put 

a; = 42°+«i y=18°+2, (B) 

and expanding the different functions, we find 

elk die 

k^ = sin 42° + cos 36°, ~ = cos 42°, ^' = - 2 sin 36° 

= 1.48 =.74 =-1.18 

dk dJc 

*2 = cos 42° + 3 sin 18°, g^^ = -sin42°, ^' = 3 cos 18° 

= 1.67 =-.67 =2.85 

*=-\l80/ + S80' dx-\m' Sy-^ 
= 2.11 = 1.47 = 5 

ki-Mi = - .02, k^-M^ = - .03, ks—M, = .01 
= m,, = m,, = ni. 
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Therefore we have for 

OBSERVATION EQUATIONS 

.74% - l.lSgj - .02 = vi 
- .6731 -I- 2.85S2 — .03 = v^ (C) 

1.47zi -|- 5S2 + .01 = Wj 

From these are obtained the 



NORMAL EQUATIONS 

3.16si + 4.5782 + .020 = 
4.5781 4- 34.51S2 - .012 = 

Solving, 

Zy= - .00845 zj = .00147 (E) 



(D) 



These results are in circular measure. Reducing to degree 
measure we have 

Si= -29'.0 S2=5'.l (F) 

Substituting these in (B), 

X = 41° Sl'.O y = 18° 5'.1 (G) 

46. If the observations are conditioned, precisely the same 
method will be followed in reducing all the equations to the 
linear form. The rest of the solution will then be as usual. 

If the values found for »i, z^, . . . Zq should turn out to be 
so large that the terms involving their second and higher 
powers can not be neglected as assumed, the process must be 
repeated using the values of Zi, Z^, . . . Zg first obtained as 
approximations. 

In a few cases the equations may be reduced to the linear 
form by some special artifice of a simple character, as in 
paragraph 43. In such cases the method of expansion by 
Taylor's Theorem should not be used. 



CHAPTER III. 

THE PRECISION OF OBSERVATIONS. 

47. The work so far considered has treated solely of the 
methods by means of which the most probable values of the 
unknown quantities may be determined from a series of 
observations. But in general something more than this is 
desired. We wish to know, if possible, how much reliance 
can be placed upon the results obtained, and how they com- 
pare in precision with other determinations of the same 
quantities. Preliminary to the discussion of the precision of 
our results it will be necessary to consider more fully the 
Curve of Error. 

48. The Constant Je, The law of distribution of errors of 
observation has been shown to be represented by a curve 
whose equation is 

y = k6-^'=^' 

In this, if x — 0,y — k (30) 

Therefore the constant k represents the intercept of the 
Curve of Error on the axis of Y. It is not, however, an 
independent quantity but is determined by the value of A, 
as will now be shown. 

49. To Find k in Terms of Ji. Since, as shown in para- 
graph 18, the total area between the curve and the axis 
of X' is denoted by unity, we have 

le r°°e-'''«' dxz=l or k f'^e-'^'^^dx = i (a) 
This may loe written 






2k 



(^) 
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Let t = hx, .-. dt ^ hdx. Also, when a = oo, < = oo, 
and when a; = 0, J = 0. 

r e-^"" dt = r e-^^ h dx (c) 

Multiplying this equation by 

C e-V'dt = fe-'i'^dh 
we have 

',dh 



dh 
X') Jo 

dx 



[£ e-^ dt T = £ £ e-ft^(l + »') h dx 

Jo 2(1 + 05'') Jo \ JK T J 

_ ^ r°° 

~2Jo l + a;' 

If ,1°° 
2L Jo 

r ,-^at = ^ (31) 

•^ 2 



IT 

4 



Substituting (31) in (c) and (b) 

V^ A A_ 

T ~ 2A *^'' ~ y/;;: 

Therefore equations (16) and (17) become 
y = —e""" 

Vir 



(32) 



(33) 
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50. It has been shown that the quantity A is a "Measure 
of Precision" of the observations, and hence a determination 
of its value in each case would enable us to compare the 
relative reliability of different measurements. In practice, 
however, it is not found convenient to compute the value of 
h directly, and so this quantity is used only in developing 
the theory of the subject, while in the comparison of observa- 
tions certain other quantities now to be derived are used as 
precision measures. These latter quantities are called the 
"Mean of the Errors," the "Mean Error," and the "Probable 
Error," respectively, and may be computed directly from the 
observations. Further, as it will be found that they all bear 
a definite relation to h, the value of this quantity can be 
determined from them if desired. 

In the following discussions no distinction is made between 
positive and negative errors of the same numerical magnitude, 
and unless otherwise stated the observations are all of the 
same weight. 

THE MEAN OF THE ERRORS OR AVERAGE DEVIATION. 

51. The Average SeTiatiou (a.d.) of an observation is 
the arithmetical mean of the errors all taken with the positive 
sign. 

Since from (6) the probability that the error of a single 
observation will fall between x and x-\- dx is 

^(a;) dx, 

if n observations are made, the number of errors falling 
between these limits is 

n <l>(x) dx. 
HencS the sum of all the errors of the observations is 



X <l>(x) dx or Ini X ^{x) 

00 "^fl 



dx. 
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Dividing this last expression through by n we have 

a.d. = 21 x<j>(x) dx = — I e~^'^^x dx 

•^0 Sw'^n 



= ("e-^^'' (-2h^x) dx 



or a.d. = — -= (34) 

feVir 



THE MEAN ERROR. 

52. The Mean Error (/x) of an observation is the square 
root of the arithmetical mean of the squares of the errors. 

The total number of errors being n, the number falling 
between x and x-\- dx is, as just shown, 

n 4>{x) dx, 

and the sum of the squares of these errors is 

n a? ^(x) dx. 

Therefore the sum of the squares of all the errors is 

n I x^ ^(x) dx 

.-. |i^ = ^ r" e-'i^^ x' dx (a) 

But as shown in paragraph 49, 

~ r e-h"^' dx=l or r e-^^ dx = ^ (b) 

Differentiating (b) with respect to h, 

-1h C e-^'^^x^dx = -^ (c) 
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and replacing the integral in (a) by its value as determined 
from equation (c), we have 



THE PROBABLE ERROR. 

53. The Probable Error (r) of an observation is an error 
such that one-half the errors of the series are greater than it 
and the other half less than it. Or it is an error of such a 
magnitude that the probability of making an error greater 
than it in any given observation is just equal to the probability 
of making one less than it, both probabilities being one-half. 

The probability that the error of an observation will fall 
between x and x-\- dx being <f>{x) dx, the probability 
that the error will fall between the limits r and — r is 

JP = f''^(x)dx = — fe-'^^^dx (36) 

If r is the probable error, JP is one-half, or 

-j; e-'^^Ux = \ (37) 

and from this delSnite integral r is to be found. 

Let i=:Ax, .-. dt = hdx. Also when x = r, we have 
t = hr, and when x = 0, < = 0. Substituting these results 
in (37), we have 

2 r*' ,2 , 1 
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Denote hr by p. Then . by interpolation in a table of 
values of this integral, the value of hr in (38) is found to 
be 

p = hr = .47694 (39) 

P .47694 ^... 

h h 

54. When t is small the values of j e~^ dt may be 

found by expanding e~^ into a series and integrating the 
successive terms. Thus, by Maclaurin's Theorem, 

^ ^ [2^ ^ [3 ^ 

... Ce-^dt = r(l-*^+ J1-J1 + . ..)<?. 
^0 -'o ^ |2_ |3_ 

= t- iV — - — + ... 

3 5[2_ 7[3_ 



RELATIONS BETWEEN |1, r, a.d., A, AND p. 

55. From (40) r = ^, and from (34) a.<?. = , 

r ^ p a.d. sfir 

or r = .8453 a.d. (41) 

1 

Also from (35) ft = , 

h\j2 

r = pp.\i2 

or r = .6745 |l (42) 
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The relation between the values of /i, r, a.d., and h may 
be conveniently expressed as follows : — 



(43) 





,x^ = - = 


: — = a.d, t/TT 
P ^ 


or arranging in 


tabular form, 








h 


r 


a.d. 


/* 


1.0000 


1.4826 


1.2533 


r 


0.6745 


1.0000 


0.8453 


a.d. 


0.7979 


1.1829 


1.0000 



From this it will be seen that 

jn > a.d. > r (44) 

56. Further, since in (25) it was shown that p oc W, it 
follows at once from (43) that 

That is, the weights of different determinations of a 
quantity vary inversely as the squares of their Mean Errors, 
their JProbable Errors, or their Average Deviations. 

It is to be observed, however, that the determination of the 
relative weights of quantities from a comparison of their 
precision measures according to (45) applies only when the 
quantities are of the same kind and subject to the same con- 
stant errors, if any of the latter exist. (See § 98.) The 
applications of (45) are numerous and important. 

57. Example A. Suppose n direct observations, all of 
the same weight, be made upon a quantity, and that the 
probable error of a single observation is r. Then since the 
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weight of the arithmetical mean is n, its probable error rj 
will be given by 

^ = 1 or U=^ (46) 

Or in general, suppose 8 is any precision measure of an 
observation of weight p, and suppose p^ is the weight of a 
second similar quantity or observation, then the correspond- 
ing precision measure So of the latter will be 

8o = 8 f^ (47) 

The case of most common occurrence is that in which 
j» = 1, and then we have 

So = 7^ (48) 

Example B. A line is measured five times and the aver- 
age deviation of the mean (A. D.) found to be .016 feet. 
How many additional measurements are necessary in order 
that the A. D. of the mean may be reduced to .004 feet ? 

Let X be the total number of observations required. Then 

X : b = .000256 : .000016 

X = SO 

Consequently the number of additional measurements 
required is 75. 

Example C. In two determinations of the quantity i 
there were obtained 

ii = 427.320 ± 0.040, Z^ = 427.30 ± 0.16 

Find their relative weights, and the most probable value 
of i and its probable error. 
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Note. The above is the method coramonly employed to 
denote that the probable errors of the observations are 0.040 
and 0.16. 

», 162 16 

From (45) — = ' = — 

^ ' iJj 42 1 

From (3), the most probable value of L will be given by 

Z„ = 427 + ^' X -^^ + -^^ 
^ 17 

= 427.319 

and the weight of Zq being 17, by (48), 

.16 
U = = .039 

Therefore we should write the result 

Zo = 427.319 ± .039 



REPRESENTATION OF |Jl,, «.<!., AND r ON THE CURVE OF 

ERROR. 

58. To find the points of inflection of the Curve of Error 
we have 

y = Ice'^^^ 

ax 

For a point of mflection, -^— 5 = 0, [or 

IK'Te e-^'^' (2 A'' a;^ - 1) = 

1 
CK = -— =. /x by (35). 

AV2 
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That is, the Mean Error is represented by the abscissa of 
the point of inflection of the Curve of Error. See OJtf in 
figure 2, page 9. 

Next, for the abscissa of the centre of gravity of the area 
to the r^ht of OY, we have 

I y X dx 
Xo=OD = ^ 

I y ^0" 

2 h /«oo 
= :?= f e-^'^xdx 

Sir Jo 

for, (§ 18), J^"y ^^=\ 
Integrating, 



1 



a.d. by (34). 



Finally, if an ordinate PP' be drawn so as to bisect the 
area to the right of the origin between the curve and the 
axis of X, the Probable Error will be represented by the dis- 
tance of this ordinate from the axis of Y. For the proba- 
bility of the occurrence of an error less than the amount OP 
is then equal to the probability of the occurrence of an error 
greater than OP. This being the case, OP is the probable 
error by definition. 



CHAPTEE IV. 

OOMPUTATION OF THE PRECISION MEASURES. 
DIRECT OBSERVATIONS. 

59. Observations of Equal Weight. Given n direct obser- 
vations all of the same weight on a single quantity M, to 
find the Mean and Probable Errors and Average Deviation 
of a single observation and of the Arithmetical Mean. 

Let the observations be M-^, M^, . . . M^. 

" " arithmetical mean be M„. 

" " real errors be %, x^, . . . x„. 

" " residuals be Vi, v^, . . . v^. 

Denote the mean and probable errors and average devia- 
tion of a single observation by /*, r, and a.d., respectively, 
and the corresponding quantities for the arithmetical mean 
hy fio) ''oj and A.D. Then by definition, 



Wl = Jfl - 


-Mo, 


v^ = Mi — Mo, . . 


. v„ = Mn— Mo 


and 









If Mo represented the true value of M, the residuals 
would be the same as the real errors, and we should have 



=v 



n 



and if n is large this formula is practically exact. But when 
n is small a more accurate expression is necessary. To 
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obtain this let M^ -\- Xg be the true value of M. Therefore 

CBi = Jfi — (Mo 4- Xo) = Wi — a!o 
asj = it^ — (Mo + Xq) = v^ — Xo 

«n = Ml - (M, + a'o) = W« - SBo 

Squaring, adding, and dividing by n 

= ^' + Xo' for by (2) S« = 

The value of Xg is not known and can not be found 
exactly, but it is approximately equal to the mean error of 

Mo, that is, by (48), to /*„ = -^■ 

Substituting this value in the above we have 

.^ ^ , ^ 
^ n ' n 

(n - 1)11.'' = 2w' (49) 






(50) 






by (42) ». = .6745 SJ ^^ (52) 



and n = .6745 V ^, j'_ ,, (53) 



n(n— 1) 
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■ 60. In order to avoid the use of the squares of the 
residuals we may proceed as follows : From (49) 



%v' = (n - 1) 



n 



On the average, the values of the residuals will then be 



j n - 1 
ln—1 



«»= V 



n - 1 

— T. — "^ 



Adding and dividing by n, neglecting the signs of the 
residuals, 



S« /w— ISa; j n — 1 

n \ n n \ n 



a.d. 



a.d. = — (54) 

^ n (n — 1) 

by (48) A.D. = — (55) 

n \/n — 1 

v /^ix .8453 Si? ,„, 

by (41) r = - (56) 

\/n {n — 1) 

, .8453 Sv 

and To = — ^zr (57) 

n^n — 1 

The mean errors may also be computed from the above by 
using the table of equivalents in paragraph 55, but this is 
not customary, formulas (50) and (51) being used for this 
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purpose. Results derived from (50) are to be regarded as 
more accurate than those obtained from (54), the latter 
being a second approximation. 

61, Example. From the following measurements on the 
length of a base line find the most probable length and the 
values of the various precision measures. 



M 


V 


v^ 


455.35 


.02 


.0004 


.35 


.02 


4 


.20 


- .13 


169 


.05 


- .28 


784 


.75 


.42 


1764 


.40 


.07 


49 


.10 


-.23 


529 


.30 


- .03 


9 


.50 


.17 


289 


.30 


- .03 


9 


3.30 


+ .70 • 


.3610 


455.330 


-.70 




M, 


Sw = 


%v' 



By (50) and (51) 
^ / .3610 

r = .6745 ^ 

By (54) and (55) 

1.40 
a.d. = 



V'90 
r — .8453 a.d. 



.20 
.13 

.15 
.13 



/^ 



V'lO 
ro = .6745 Ha 



A.D. 



1.40 
"30" 

.8453 A.D. 



.063 
.042 

.047 
.040 
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and we should write for the most probable length of the 
base line 

i!fo = 455.330 ± .042 

62. Observations of Unequal Weight. Using the same 
notation as above, with slight modifications, we will 

Let Ma represent the General Mean. 
" i'u Pii • • • Pn l^e the weights. 
" a.d.i, a.d.^, . . . «.<?.„ be the average deviations. 
" Hi, ix^, ... fin be the mean errors. 
" r-i, r^, . . . r„ be the probable errors. 
" a.d., ju., r, and v refer to observations of weight unity. 

Then by (48) 

/*o = , 1^1. = —;=., etc. 

\f^ \/Pk 

If the " Observation Equations " are formed for this case 
they will be 

Ml — Mo = Vi, Mi — Mo ^= v„ . . . Mn — Mo = Vn. 

And, as was shown in paragraph 28, if these equations are 
each multiplied by the square root of the weight of the corre- 
sponding observation, they will all be redftced to equivalent 
equations of weight unity. On performing this operation 
it will be seen that the residuals of the new equations 
become 

■»1 sTpx, Vi nTJb^, ... Vn sTPn 

And evidently to these reduced observations the formulas 
of paragraph 59 apply. 
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Therefore 



.=V" 



^pv^ 



n —\ 



r = .6745 |i. (58) 



^^ = ^ pjr-x) *•* = '^''' ^^ (^»> 



Also, by a method similar to that used in paragraph 60, it 
may be shown that 

a.d. = - — ^ r = MSt^a.d. (61) 

\l n{n — 1) 



a.tJ.j, = ^ n = .8453 a.d.,, (62) 

S Pf,n(n — 1) 



JL.JD. = ^^ -^ — n = .8453 4.X>. (63) 

)J'Spn{n — 1) 



Formulas (61), (62) and (63) are not of much value in 
practice unless all the weights are perfect squares, for other- 
wise no labor is saved in the computations, since the square 
root of each weight will have to be determined, and the results 
obtained can not be considered as reliable as those given by 
formulas (58), (59) and (60). 
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63. JExample. Given a series of observations on M, the 
difference in longitude between two stations. To find M^, 
IX, /xo, a.d., A.J)., r, r^, r^. 



M 


P 


pM 


V 


v\[^ 


v" 


pv'' 


V 4' 30" 


4 


120 


1.6 


3.2 


2.6 


10.4 


41 


1 


41 


9.4 


9.4 


88.4 


88.4 


43 


1 


43 


11.4 


11.4 


130.0 


130.0 


37 


9 


333 


5.4 


16.2 


29.1 


261.9 


48 


4 


192 


16.4 


32.8 


269.0 


1076.0 


34 


16 


544 


2.4 


9.6 


5.7 


91.2 


25 


9 


225 


6.6 


19.8 


43.6 


392.4 


46 


1 


46 


14.4 


14.4 


207.4 


207.4 


28 


25 


700 


3.6 


18.0 


13.0 


325.0 


24 


4 


96 


7.6 


15.2 


57.7 


230.8 


1» 4' 31".6 


74 


2340 




150.0 




2813.5 


M, 


^P 


'Sipm 




%v\[p 




^pv" 



2 vM 2340" 

By (3) i[f„ = ^ = lo 4' + 



/ 



9 



%p ' 74 

= 1° 4' 31".6 



2813.5 



= 17.7 



/^ = -A- = 2.1 



\/74 



.6745 u. = 11.9 



n = - = 4.0 



^90 



A.I).^^= 1. 

V/74 



r, = .6745 jiio = 1-4 = .8453 A.D. = 1.5 



Mo 



r 4' 31". 6 ± 1".4 
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FUNCTIONS OF OBSERVED QUANTITIES. 

64. Theorem. Given any number of quantities and their 
Mean and Probable Errors and Average Deviations, to find 
the Mean and Probable Errors and Average Deviation of any 
function of the quantities. 

Let the quantities be Jif^, M^, . . . Mg. 
" mean errors be ju.i, fi^, . . . fiq. 
« function be M=f{M^,M^, . .. M^). 

" mean error of M be E. 
" probable error of M be It. 
" average deviation of M be D. 

The derivation of the general formula will be simplified if 
we consider first a few special forms of functions. 

65. Case I. Suppose M = Mx ± M^- 

The number of observations from which M^ and M^, and 
hence fi-^ and /tj, have been determined is not necessarily 
known, but we may assume that for each quantity it is any 
large number n, and that the real errors of the observations 
are 

for Ml, Xi', Xi", iBj'", . . . 

Then the real errors of M, computed from the separate 
observations on M^ and M^ will be 

ail' ± x^, Xi' ± sca", asi'" ± a^'" . . . 



E^ = 



Saji^i 2SaHa;2 + Sa;/ 



or W = |Jii^ + 11^^ (64) 

since in the most probable case the term Sa^ajj will dis- 
appear, as there most likely will be as many positive as 
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negative products of the same absolute magnitude of the 
form aji x^. 

By successive applications of the above the same principle 
may be extended to cover the algebraic sum of any number 
of quantities. So that if 

then :E' = iii^ + |i,2^ + . . . iij" = S |JL» (65) 

Since probable errors and average deviations difier from 

mean errors merely by a constant factor, we shall likewise 

have 

_g2 ^ ^^2 _|_ ^^2 _|_ _ _ ^^2 ^ jr'' (66) 

D^ = a.d.^^ + a.d.^^ + . . . = -Za.d.^ (6T) 

66. Example. Given the telegraphic longitude results, 

ft. m. sec. sec. 

(A) Cambridge west of Greenwich, 4 44 30.99 ± 0.23 

(B) Omaha " " Cambridge, 1 39 15.04 ± 0.06 

(C) Springfield east " Omaha, 25 8.69 ± 0.11 

Find X, the longitude of Springfield, and its probable error. 

i = ^ + ^ — C = 5* 58*" 37'.34 ± 0'.26 
for by (66) R = ^(.28" + .Oe'' + .11^) = .26 

67. Case II. Suppose M = ttiMi 

Using the same notation as in paragraph 64, the real 



errors of M will be 




aiKi', aiKi", aiiBi'", 




n 


<hW 


or ^ = ailJii 





(68) 
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68. Case III. Suppose M =^ a-iMx ± a^Mi ■ • . ttqMq 
By combining (65) and (68) 

also M^ = ^aW^ (69) 

and J>^ = 'Za^a.d.' 

69. Example. The length of a bar at 20° Centigrade is 
found to be 75".0041 ± 0".0037, and its expansion per 
centigrade degree 0".0036 ± 0".0018. What is the length 
of the bar at 56° Fahrenheit ? 

20° C. = - X 20 + 32 = 68° F. 
5 

5 5 

The expansion for 1° F. will be - X -0036 ± - X -0018 

y y 

or .0020 ± .0010 

Z = 75.0041 - 12 X .0020 
= 74.9801 



B = ^.0087"+ (12 X .001)" 

= .013 
Z = 74".980 ± 0".013 

70. Case IV. Suppose M = f {Mx, M^, . . . Mg). 

Let Mi = ai-{- mi, M^ = a^ -}- m^, . . . M^ = a^ -\- mg, 
where a^, a^, . . . ag are arbitrarily assumed quantities very 
nearly equal to Mx, M^, . . . Mg, so that mj, m^, . . . m^ are 
so small that their second and higher powers may be neg- 
lected. Then the errors of Mx, M^, . . . Mg may be con- 
sidered to be in mx, m^, . . . mg, and hence /ni, /j^, . . , /xg 
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may be regarded as the mean errors of the quantities 
m,, m^, . . . nig. We now have 

M := f («! + mi, aj + jWjj, . . . a, + Wg) 

Expanding this expression by Taylor's Theorem, and 
denoting f{ai, a^, . . . ag) by M', we have 

M = M' 

BM' , SM' , dM' 

' ' 5% ' " da^ ' « dag ^ ' 

-\- negligible terms in the second and 
higher powers of wij, Wa, . . . nig. 

Then the mean error of M will be the same as the m^an 
error of terms (A), and by (69) this is given by 

or what is practically the same thing, by 

- -— -[sj +■•■•-.• [:fj<"> 

71. Example A. Two sides of a right triangle are meas- 
ured with results 

a = 49.53 ± 0.59 
* = 50.38 ± 0.93 
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Find the length of the hypotenuse c and its probable 
error. In this example 



M — V Mi" + M^^ = sja"- + b"- = 70.65 
BM _ a BM _ b 

By (70) R^ = r,= ^, + ./ ^, 



r49.53n'' , „ „ r5o.38n2 
.59^ — + .93^ — 

L7O.65J ^ [_70.65j 



It = .78 
and c = 70.65 ± 0.78 

JSxample JB. If the probable error of x is r, what is 
the probable error of the common logarithm of a; ? In this 
case 

M = logio Mi = logio X 

BM B logio Jg logio e 

BMx dx X 

it = — r 



Example C. If the weight of x is p, what is the weight 

Po of sin i" ? 

Denoting the mean errors of x and sin a; by ju, and E, 
respectively, we have, by (45) and (70), 

„ rs sin x~Y 

E^ = y? \ = i>? C08= X 
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and Pi = !^ 

p W 

P 

Piy = = p sec' X 



72. Equation (70), which expresses the law of propaga- 
tion of error in functions of observed quantities, is one of the 
most important in the whole theory of the Method of Least 
Squares. Upon it in particular is based the discussion of the 
" Precision of Measurements." This subject treats, in the 
first place, of the methods of finding the precision of a 
quantity obtained by computation from a series of measured 
quantities ; and, in the second place, it investigates the pre- 
cision with which the component measurements of a series 
must be made in order to obtain a required degree of pre- 
cision in the final result. The following simple example will 
illustrate the character of the solutions : — 

73. Example. In the determination of a current by a 
tangent galvanometer we have 

I = 10 — tan <A 
G ^ 

where T is the current in amperes, £r the horizontal com- 
ponent of the earth's magnetic force, G the galvanometer 
constant, and <^ the angle of deflection. Given the errors 
§1, §2, Sg, in S, G and tan ^, to find the error A in Z 

By (70) 

100 „, , , IddH^ 100^2 

= -^ ^' tan^ ^ + —^^ S^ tan^ ,^ + -^^ 83= (a) 
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100^2 
Dividing this equation hj I^ = tan^ <^ we have 

That is, the square of the percentage error in I is equal 
to the sum of the squares of the percentage errors in S, G, 
and tan ^. Hence if 

S is determined within .4 per cent, 
tan <^ " " " .1 " " 

then — = V -16 + .04 + .01 = .46 per cent. (c) 

Next, suppose the value of I is required to within .1 per 
cent. To find the necessary accuracy in the determinations 
of S, O, and tan ^ when the error in each of these quan- 
tities is to have the same influence upon the total error. 

From (b) we shall now have 

A = -^ =. .000577 

8i = .00058 J2" 

Sj = .00058 O (e) 

S3 = .00058 tan <^ 

It is comparatively easy to obtain the necessary accuracy 
in the measurements of G and tan <^, but difficult in the 
case of M. 

For additional work of this kind see Holman's " Precision 
of Measurements." 
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74. Combination of Functions of the Same Variables. 
It is to be noticed tbat equation (70) applies only when M 
is a function of independent quantities. If M^, M^, . . . M^ 
are merely different functions of the same quantities we 
must proceed as follows : — 

Let Mx — 4, (zi, «2, . . . sj,) 

Mi — xj/ (%, Zj, . . . zj) 
M = f ( Jfx, M,) 

If any single observations of Zj, Zj, . . . Zjc are subject to 
errors ajj, x^, . . . ajj, the corresponding errors in M^ and 
Mi will be 

for J^, Xi = ffliaji + a^x^ + • ■ ■ «*"'* (a) 

« M^, Xi = tti'lKi + fifa'^'a + • • • «*''"* (15 ) 

Where %, «2, . . . a^ are the differential coefficients of 
Mi, and aj', a^ . . . a^ the differential coefficients of i!^ 
with respect to Zj, Zj, . . . Zj. The corresponding error in 
M will then be 

X = AXi + A'Xi (c) 

Where A and ^' are the differential coefficients of M 
with respect to Mi and Jlf^. Substituting in (c) from (a) 
and (b), 

X = (Aui -\- A'a\) aji + {Aa^ -\- Aa\) x^ . . . 
= aaii -|- pXi + - . . Xxjc 

Then if the number of observations or values of X be 
denoted by n, 

2 -T' 2 Ki' %Xi^ 



n 



+ 2a/3 ^ + 



n 
E^ = a.W + PW + ■■■ ^Vk' (d) 
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since in the most probable case the product terms will 
cancel out. 

Expanding (d) we have 

E^ = {Aa^ + A'a^'f ^^ + {Aa^ + A'a^'Y i>^^ + . . . 

= A\a,W + a^W ■■•) + ^'\<W + «2'V • • •) 
+ 2^J.'(aiffliV+«AV+...) (71) 

75. Example. As a very simple problem take 

Ml = 2zi, M, = 3%, jai = 0.1, 
and M = Mi-\- M^. 

Then A = \, A' = 1, a-^ = 2, ai' = 3. 
By (71) ^2 = 4 X .01 + 9 X .01 + 2 X 2 X 3 X .01 = .25 
or E = 0.5 

In this particular example the result may be found directly 
from (68) by substituting at first in M the values of Mi 
and M^. 

Thus M = 2si + 3zi = 5si 

E = 5^ = 0.5 

If J!^ and M^ had been independent quantities, by (64) 
or (69) we should have had 

E = V(2 X 0.1)^4- (3 X 0.1)^ 
= 0.36 

INDIRECT OBSERVATIONS. 

76. The determination of the precision measures of the 
unknown quantities in case the observations are indirect 
involves a knowledge of the weights of the unknowns, and 
consequently the method of computing these weights must 
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first be demonstrated. It will he assumed at present that all 
the observations are of weight unity. 

FIRST METHOD OF COMPUTING THE WEIGHTS. 

Let the observations be 

ai% + *iZ2 + • ■ • ^iZa + *i = -^ 



«2«i + *2S2 + • • • q^eg + h 



M, 



atfil + *A + • ■ • qrflq + *m = -3^ 

in which M^, M^, . . . M^ denote the actual observations, 
and Si, %, ... Zq the most probable values of the unknown 
quantities. Let 

*i — i!fi = Wi, k^— M^ =1 m^, . . . kn— Mn = m^ 
Then the above equations give rise to the 

OBSERVATION EQUATIONS 

«iSi + *i% + • • • q\Zq + m, = ■yj 



«2% + h^i + 



««2i + irfii + . 



q^Zq -\- m^ = Vi 



qrflq + m„ = Vn 



(A) 



By the rule, paragraph 25, we now form the 
NORMAL EQUATIONS 

Zi S aS + Ka 2 i^ + . . . Zy S 5? + 2 Swi = 



(B) 



Zi S ay + Zj 2 J? + . . . Sg 2 S'^ + 2 ?»» =: 
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Multiply the first of (B) by Qi, the second of (B) by 
Q2, • ■ ■ and add the results. Then 



^iCSiSa^ + Q.^ctb + . 


.. Qg^aq) 


+ s, (^iSaS + e^SJ' + . 


■■ Qa^H) 


+ 0g(CiSa?+ Q.^bq + . 


■ ■Qa^sn 


4- §iS am + Q^%bm -{- . . 
= 


. Qqtqm 



(C) 



Let Qi, Qi, ■ ■ ■ 'Qg be determined so that the coefficient 
of 8i in (C) shall be unity, and the coefficients of z^, 
Sj, ... Zg each equal to zero. That is, let 

Qi%a^ + Q^^ab + . . . Q^^aq -1 = 

^iS abJ^Qab^ + • . . Qq^ bq = 

(D) 

Qr%aq-\- Q,%bq -{-... Qg^q'' = 

Then (C) becomes 

3v + 61S ««* + 62S Sm + . . . §g2 g-m = (E) 

Equations (D) may be derived from the normal equations 
(B), if in them we replace Zj, Zj, . . . z^ by Q^, Q^, . . . Qg, 
2 am by — 1, and 2 5m, 5 cm, . . . % qm. by zero. Hence 
the solution of the normal equations with these changes will 
give the values of Qi, Q^, ... Qg. 

TO SHOW THAT Qi IS THE RECIPROCAL OF THE WEIGHT OF «i. 
Expanding the coefficients of (E) we have 



COMPUTATION OF THE PBEOISION MEASURES. 65 



^1 + Cl(«l»»l + «2»W2 + • • • ««»*«) 

+ Qiihmi + hm^ + • • • *n»^«) 

+ QqiSifni + ^jms + . . . S'nmn) 
= 



or collecting the coefficients of m^, m^, . . . m„ 

zi + mi( Q^ai + Q^bi + . . . Q^qi) 
+ m^(Qia2 + §jSj + ■ • . Q,qi) 

+ m^iQian + QJ>n + • • ■ Qqq^) 

— 



(E') 



For convenience in writing we will let 

ai =: giQSi + QJ)^ + • ■ • 6g?l 

"2 = 6l«2 + ^2*2 + ■ ■ • Qq^i 



On,. = Ql^n + §2*71 + • • • QqQn 



(F) 



Zi + aim, + 02^2 + • ■ • "n"*" = *^ (Gr) 

Multiply the first of (F) by a,, the second by a^, ... 
and add the results. Then 

■S,aa = Qita' -^ Q^%ab -{- . . . §,2 aq 

— 1 by first of (D) (H) 

Multiply the first of (F) by 5,, the second by b^, ... 
and add the results. Then 
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%ba. = Q{S,ab + Q^%¥ -\- ... Qgtbq 
= by second of (D) 

Likewise S ca = (H') 



S^a = 

Multiply the first of (F) by ai, the second by ouj, , . . 
and add the results. This gives 

Sa^ = §iS Oa + §2^ 5a + . . . QgS, qa. 

= Q^ by (H) and (H') (I) 

Let /* be the mean error of an observation of weight unity. 
Let /Ag, be the mean error of z^. 
Let Pz^ be the weight of z^. 

The mean errors of m^i, wij, . . . m^ are the same as the 
mean errors of J^, M^, . . . M^ and each is accordingly 
equal to ju. Therefore from (G), by (69) 

2 2 2 1 2 2 1 2 2 

= Q^^ by (I) (J) 

But by (48) ;a,.= = ^ 

Comparing this with (J) we see at once that 

«^ = K <^) 

Therefore, for the First Method of computing the weights 
we have the following : — 
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77. Rule I. In the normal equation for z^ write — 1 
for the absolute term % am, and in the other equations 
zero for each of the absolute terms % bm,, 2 cm, . . . S qm. 
The value of Zi found from, these equations, is the recipro- 
cal of the weight of the value of Zj obtained by the solution 
of the normal equations. 

To Jind the weights of z^, z^, . . . Zg, proceed in a similar 
way, forming a corresponding set of equations for each 
unknown. 

SECOND METHOD OF COMPUTING THE WEIGHTS. 

78. Write equations (B) of paragraph 76 in the form 

Si S a^ + ^2 2 aft + ... Zg'Siaq -\- % am = A 
Zi%ab-\- z^'Sib^ -\- . . . Zg'S,bq -\- ■Sibm = B 

(L) 

Zi S a? + Za 2 *§- + . . . 2g 2 g''^ -\- 'Xqm = Q 

Then in the solution by the preceding method, equation 
(E) becomes 

2i + QA ain + Qi% bm -\- . . . Qg% qm 

= §1^ + ^2-5 + ... QgQ (M) 

in which, as was proved in (K), Q-^ is the reciprocal of the 
weight of Zi. Whatever method of elimination is employed 
in the solution of the normal equations, the coefficient of A 
in the value of Zj must necessarily be always the same. 
Hence we have 

79. Rule II. Write A, B, . . . Q instead of zero in the 
second mem,bers of the norm,al equations and carry out 
their solution in any convenient way. Then the most 
probable values of Zj, z^, . . . z are given by those terms 
in the results which are independent of A, B. . . . Q, 
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The weight of z^ is the reciprocal of the coefficient of A 
in the value of Zj. The weight of z^ is the reciprocal of 
the coefficient of B in the value of z^, etc., etc. 

THIRD METHOD OF COMPUTING THE WEIGHTS. 

80. From the second, third, . . . equations of (L) find 
the values of z^, z^, ... Zj in terms of Sj and substitute 
in the first of (L) without reduction. Then the first of 
(L) becomes 

Bzi — T -\- A -\- terms in B, C, . . . Q 

Where T is the sum of all the numerical quantities result- 
ing from the substitutions. Dividing through by B, 

Zi = 4- _ + terms m B, C, . . . Q (N) 
si Jti 

T 

in which — is the most probable value of Sj, and, as was 
B 

shown in deriving the second method. 



(O) 



81. Eule III. Substitute in the normal equation for Sj 
the values of z^, ^a, • • • Zg in terms of z^ as found from, 
the remaining equations. Then before freeing of fractions 
or introducing any reduction factor, the coefficient of Sj 



1 

B ~ 


«i 


B = 


P^ 


From this follows at once 
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in this equation is the weight of the value of z^ obtained 
in the solution. 

To find the weights of Sj, e^, . . . Zq, proceed in a similar 
way with the normal equations for each of these unknowns. 

For the solution of an example by the three different 
methods see paragraph 84. 

THE MEAN ERROR OF AN OBSERVATION. 

82. The next step will be to derive /*, the mean error 
of an observation of weight unity. In the following demon- 
stration the equations referred to by letters are those . in 
paragraphs 76 to 81. 

Let the real values of Zj, z^, ... Zq be 

«1 + «h., Z2 + Ka, . . . Zq -\- Xg 

and substituting in (A) we have 

«i(% + aa) + *i(z2 + x^) + ■ ■ ■ gi(Zg -f a;,) -]- m^ = A^ 
a^izi + aji) + bi{z^ J^ x^) -\- . . . q^{Zq -f Xq) + m^ = A^ 

(P) 

a„(gi -f Xi) 4- 5„(Z2 -\-x^) -\- . . . q^{Zq -f- Xq) 4- m„ = A„ 

where Aj, Aj, . . . A„ are the real errors of M^, M^, . . . M^, 
or of Wi, m^, . . . m^. 

Multiply the first of (P) by a^, the second by flj, ... 
and add the results. This gives 

Zi 2 a'' + z^%ab -\- ... Zq^aq -\- 2 am 
-\- Xj^ta'^ -\- x^'tab -\- . . . Xq'taq = S a A 

But by the first of (B) the first line in this equation is 
equal to zero, and therefore 
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x^'S, a^ -\- Xi'Si ab -\- . . . Xq'S, aq — % a^ = 
Also cBi S a5 -I- jca S S^ + . . . a;„ S Sg- — 2 5A = 

(Q) 

Xi'%aq -\- Xi'Sibq -\- . . . Xq%^ — S 2'A = 

These being of the same form as the normal equations 
(B), the value of Xy resulting from their solution will be 
of the same form as that of Si in the solution of those equa- 
tions, with only the substitution of — A for m. From 
(G) we shall therefore have 

Xi — ajAi — ttjAa — ... o,iA„ = (R) 

Multiply the first of (P) by v^, the second of (P) by 
«2, . . ■ and add the results. Then 

■ (^1 + Xi)%av -\- (Z2 + Xi)-S,bv + . . , (Zg -f Xq)%qv 
-\- ^mv = SAz) 

and multiplying the first of (A) by fflj, the second by 
a^, . . . and adding the results 

S au = Zi'S,a^ -\- z^'Siab -]-... Zq^aq -\- % am 
— by the first of (B). 

Also SSy = 

S?w = 
Substituting these values in the above, we find that 

Simv = SAW (S) 

Now multiply the first of (A) by v^, the second by 
v., . . . and add the results. This shows that 
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Zi'Zav -|- Zi%hv -\- ... Zj 2 5"w + 2 mv = %v^ 
and as above, 

2 aw = %bv . . . = 2 S'w == 
Combining this result with (S) we have 

%mv = 2^)^ ~ 2 Aw (T) 

Next multiply the first of (A) by A^, the second by 
Aj, . . . , add the results and compare with (T). This gives 

Zi 2 aA -|- Zj 2 5a + • ■ • 2g 2 2'A -|- 2 wiA 

= 2 Aw — 2w^ (U) 

And finally, multiplying the first of (P) by Aj, the second 
by Aj, ... and adding the results, we have 

Zi 2 «A -\- Za 2 5a -\- ... Zj 2 g-A -|- 2 mA. 
-\- Xi-S,aA -\- x^-S,bA -If- . . . Xg^S,q^ = 2 A" 

Therefore from (U) 

2 w^ + Ki 2 aA + aij 2 5A + . . . Kj 2 2'A = %A'' (V) 
, ^ 2A^ ^ 2w= aii2aA + . . . a;^ 2 gA 

We must now find the mean values of the terms x-^ 2 aA, 
CBj 2 JA, ... a;, 2 jA. Expanding 2 aA, 

2 aA =: aiAi -|- a^A^ + ■ • • a„A„ 
from (R) ail = aiA, -\- a^A^ + • • • 'hAri 

Multiplying, 

ail 2 aA = aiaiAj^ -|- a^a^A.^ + • • • a„a„A„' 
-|- terms in AjA^, AiA,, . . . 
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In the most probable case these product terms vanish, and 
substituting for Aj", A^", . . . A„^ the mean value /*", we 
have 

Similarly x^libA = ft!^ 



Kg S g'A = /*" 
From (W) then, q being the number of unknowns, 

n n 






n — q 



(72) 



1^. = J^ = \l ^""^ . (73) 



^f^ y p^{n- q) 



r, = .6745 Ji. = .6745 V — ^^^ 

y Pz{n — 



Q) 



(74) 



83. By a method similar to that used in paragraph 60, we 
may derive 

a.^. = (75) 

\l n{n — q) 

Si? 
a.d.z = - (76) 

>J Pzn{n - q) 

84. Example. In illustration of the above processes 
take the example in paragraph 24, where we found for 
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OBSERVATION EQUATIONS 



and for 



Sl — ^2 — 1-7 = ■Ml 

Zg — 2.4 = ■Wj 

— Si + % + Zs — 1-0 = ■"s 

% — ^8 — 3.0 = Vi 

NORMAL EQUATIONS 



(A) 



Z3 — 0.7 = (a) 

- 2 Zi + 3 Z2 _ 2.3 = (b) 

— Si + 8 Zs — 0.4 = (c) 

SOLUTION FOR THE WEIGHTS BY THE FIRST METHOD. 

For finding the weight of Zj the above normal equations 
would be written 





2 -^1 - 2 z, - Z3 - 1 

- 2 zi 4- 3 Z2 

— 2l + 3 Zj 


— 







(a') 

(b') 
(c') 


Solve for 


2l 








3 X (a') 


6 Zi — 6 Z2 — 3 Z3 — 3 

- zi + 3 Z3 











2 X (h') 


5 Zi — 6 Z2 — 3 
- 4 zi + 6 z. 


=: 


~0 







Zi — 3 


= 









Zi = 3 .-. ^^, = 


1 




(d') 



For finding the weight of z^ the equations are 



2 Zi — 2 Zj — Z3 


= 


(a") 


- 2 Zj + 3 z. 


-1 = 


(b") 


- Zi + 3 Zs 


= 


(c") 
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Solve for z^ 

3 X (a") 6 Zi - 6 32 - 3 Zs = 

(c") - zi + 3 ^3 = 

5 Zi — 6 Z2 =0 
_ 6 

Zi ~ Zo 



5 3 

Substitute in (b") z^ = - . . p^^ = - (d") 

o 

For finding the weight of Zg the equations are 

2 Zi - 2 22 — Z4 =0 (a'") 

- 2 Zi + 3 S2 =0 (b'") 

- Si 4. 3 z, - 1 = (c'") 

Solve for z, 

4 
from (b'") 2 Z2 = - z^ 

O 

substitute in (a'") 2 Zi — 3 Zj :^ 

2 X (c'") — 2gi-|-6z3-2 = 

8 Z3 - 2 = 

2 3 

Z3 = - ... ^,3 = - (d'") 

SOLUTION BY THE SECOND METHOD. 

The normal equations will now be modifi^ed so as to appear 
in the following form : — 

2 zi - 2 Zj - Zs - 0.7 = ^ (a) 

- 2 zi + 3 Z2 - 2.3 = ^ (b) 

- Zi + 3 Zs - 0.4 = (7 (c) 
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Solving, 

3 X (a) 6 Si - 6 2j — 3 Ss - 2.1 = 3 4 

(c) - % + 3 z, - 0.4 z= C 

5 Zi - 6 S2 — 2.5 = Z A -\- C 

- 4 gi + 6 2 3 - 4.6 = 2 -g 

Zi -7.1 = 3^ + 25+C (8) 

2i = 7.1 and p^ = ~ (d) 

Substituting (8) in (b) 

3 z, = 16.5 + 6^ + 55 + 2C 

3 

2s = 5.5 and p^ =z - (e) 

a 

Substituting (8) in (c) 

3zs = 7.5 + 3^ + 2^ + 2(7 

3 

2, = 2.5 and /j^, = - (f) 

SOLUTION BY THE THIRD METHOD. 
The normal equations are now taken in their original form. 

2 2i - 2 2j - 2, - 0.7 = (a) 

- 2 Si + 3 Zj - 2.3 = (b) 

- 2i + 3 S3 _ 0.4 = (c) 
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To obtain Si and its weight we proceed as follows : 

from (c) z, ■= — -\- — 

2 Si 2.3 

from (b) g^ = __ _^ __ 

Substitute in (a) 

4 4.6 zi .4 „„ 

2zi_-zi- — - — 0.7 = 

3 3 3 3 

Collecting terms, — — -^ = 

o o 

zi = 7.1 and Pz^ = - (d) 

For Zj 

3 X (a) + (c) 5 Zi - 6 zj - 2.5 = 

6 

Zi = 7 ^2 + 0.5 

12 

Substitute in (b) Za — 1.0 + 8 Zj — 2.3 = 



3 
Collecting terms, — Zj — 3.3 ^ 
5 

3 
Zj = 5.5 and p^ ^ ~ (e) 

o 

For Zs 
3 X (a) + 2 X (b) 2 Si - 3 Z3 - 6.7 = 

3 , 6.7 

Substitute in (c) Sg — -— + 3 Sg — 0.4 := 

^ 2 
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^ „ . 3 7.5 „ 

Collecting terms, — e^ — — = 

3 

23 = 2.5 and p^ = - (f) 

It is evident that the three methods give identically the 
same results and that the work is about the same in each case. 
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Substituting the values found for z^, z^ and z, in the 
observations equations (A), we have 

7.1 - 5.5 - 1.7 = «! = - .1 .01 — V 

2.5 — 2A — v^ — -\- .1 .01 = V 

- 7.1 + 5.5 + 2.5 - 1.0 = ■Wj = - .1 .01 — V 

5.5 - 2.5 - 3.0 = Vi — .0 .00 = < 



.03 = 21)' 



In this 


example w = 4, S* = 


3. 




By (72) 




^ = \/f = 


.17 




By (74) 




r = .6745/i = 


.12 




By (73) 


/^0i 


= J^ = .30 


'•zi = 


.20 




/*Z2 


= ^J^ = .22 


r^ = 


.15 




M23 


= ^J = ■'' 


^z, = 


.10 


By (75) 




ad - -^ — 
2 


.15 





r = .8453 a.d. = .13 
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85. Observations of XTnequal Weights. If the observations 
are not all of the same weight the formulas and operations 
are merely modified in the usual manner and equations (72) 
and (75) take the more general form 



_ JlP^ (77) 

a.d. = — (78) 

V n{n — q) 



CONDITIONED OBSERVATIONS. 

86. Suppose there are given n observations, n' condi- 
tions and q unknown quantities. Then by paragraph 33, 
the method of solution is to eliminate n' unknowns between 
the " Observation " and " Condition " equations, leaving 
q — n' independent unknowns in the " Normal " equations. 
Consequently formula (77) now applies to this case and it 
would be written 



V^;,^^^ <«) 



n — q -\- n' 

also n, = y^ (80) 

The weights of the q — n' unknown quantities can be 
found by any one of the three methods already given and 
then the mean error of each unknown may be computed by 
using formula (80). If the mean errors of the n' quantities 
that were first eliminated is wanted, their weights must be 
determined by eliminating a different set of n' quantities 
from the original observation equations and solving the neces- 
sary sets of equations for these weights. The first method of 



%H 


h z^ - 3.0 = 


weight 1 


«i 


- 23 + 1.5 = 


« 4 




Z2 - 2.2 = 


« 3 


Sl 


+ Z3 - 8.4 = 


« 2 


and the 
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solution for the weights would perhaps be hest here as the 
actual values of the unknowns have already been found. 
87. Example. Given the 

OBSERVATION EQUATIONS 



(A) 



CONDITION EQUATION 

Zs - Za - 0.5 = (B) 

To find the most probable values of Zi, Zj, and Zg, and 
also their mean errors. 

Eliminating Zj between (A) and (B) there remains 

Zi -|- % — 3-0 = weight 1 

3j _ 3^ + 1.0 = "4 

Z2 — 2.2 = "3 

Zj _)_ Zj - 2.9 = "2 

From these we have the 



NORMAL EQUATIONS 

7 zi - Za — 4.8 = 
- Zi 4- 10 Za - 19.4 = 



Solving, 








zi = 0.98 


p,, = 6.9 




Za = 2.04 


Pz, = 9.9 


om (B) 


Za = 2.54 





(C) 



(D) 
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Now eliminating Zg between (A) and (B), we find 



Zi _|_ Zj _ 3.5 = 


weight 1 


z^ - z^ -\- 1.4 = 


4 


z, - 2.7 = 


3 


z^ -^ Zg — 3.4 = 


« 2 


From these we derive the new set of 





(E) 



(F) 



NORMAL EQUATIONS 

7 Zi — Zs — 4.8 = 
- Zi 4- 10 Z3 — 24.4 = 

and solving for Zg 

Z3 = 2.54 p^ = 9.9 (G) 

Substituting the values of Zj, Zj and Zj in equations (A), 
we have the residuals 



V v^ 


pv^ 






.02 .0004 


.0004 






.06 .0036 


.0144 






.16 .0256 


.0768 






.12 .0144 


.0288 








.1204 


= 


Si^w" 


Since in this example w = 4, 


n' = 
: .25 


1, 


2- = 3, 


, (79) . - v/-'f - 






: .09 







^6.9 
;.,, = -^ = .08 

H-^ = -^ = .08 
^9.9 
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The first significant figure in the mean errors is so large 
that it is not worth while to retain the second place as usual. 

If desired, the probable errors and average deviations can 
now be computed by the usual formulas. 

88. In case the observations are made directly upon the 
values of several quantities subject to certain conditions, we 
have n = q, and equation (79) reduces to 



^ = \/^ (81) 



n 



DBSERVATIO 


N EQUATIONS 


zi = 


weight 3 


Z2 = 


2 


S3 = 


2 


z. = 


" 1 



from which the mean error of any observation may at once 
be computed from its weight. 

89. jExample. Taking the example in paragraph 34 on 
the measurement of the angles of a quadrilateral, we had for 



(A) 



for the 

CONDITION EQUATION 

% + 2. + sa + Z4 + 58 = (B) 

and for the 

NORMAL EQUATIONS 

4 Sj _|_ zj _|_ Z3 + 58 = 
zi + 3 Z2 + zs + 58 = (C) 

z, + Z2 + 3 Z3 + 58 = 

Solving these equations for the values of the unknown 
quantities and also for their weights, we have 



Zi = 


— 8.3 


Pz, = 3.5 


S2 = 


- 12.4 


Pz^ = 2.5 


2s = 


- 12.4 


Pz, = 2.5 
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and forming a new set of normals containing S4, we find, on 
solving for that quantity, 

z, = — 24.9 p^^ = 1.75 

These values of Sj, Sg, 83, z^ are also in this case the 
residuals of the ohservations, and therefore to compute the 
precision measures we have 

8.3 68.9 206.7 

12.4 153.8 807.6 

12.4 153.8 307.6 

24.9 620.0 620.0 



= 1 



1441.9 = 2jJ«* 



By (81) ^ = y/llil:^ = 88 



/*«! = 



/ta, = /*a = 






= 20 r,, = 13 

= 24 r^ = r^ = 16 



/*.. = -== = 29 r^ = 20 

Vl.75 

We shall accordingly write for the most probable values of 
the angles of the quadrilateral 



A 


= 101° 


13' 


14" 


± 


13" 


B 


= 93 


49 


5 


± 


16 


G 


= 87 


5 


27 


± 


16 


D 


= 77 


52 


15 


± 


20 



In the original solution in paragraph 34 it is evident that 
the results were carried out to a greater number of places of 
significant figures than the character of the observations 
warranted. 



CHAPTEK V. 

MISOEI.LANEOUS THEOREMS. 

THE OrSTRIBUTION OF ERRORS. 

90. Having developed the processes for the adjustment of 
observations according to the Method of Least Squares, it 
will now be interesting to show how closely the distribution 
of errors found in actual practice corresponds to the theo- 
retical distribution upon which our methods of solution are 
based. 

By formula (36), the probability that the error of a single 
observation will be numerically less than a is 

P = M r\-h^^a^ (82) 

\ ir Jo 

Let < = hx, .•. dt ^^ h dx. Also when a; = 0, 

< = ; and when a; = a, t ■= ha = p—. Sub- 

r 
stituting in (82), 



JP == -^ rre'^dt (83; 

■n J 



Values of P for values of the argument — are given in 

Table I. Also for any series of observations this quantity 
P will represent the fraction of the entire number which 
should have errors less than the amount a. Hence if P is 
multiplied by the whole number of observations the result 
will be the number of errors which should be less than the 
limit a. 
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91. Example. Forty measurements on the diameter of 
Saturn's ring were made by Bessel, with the following 
results : — 



M V 


M 


V 


M 


V 


M V 


38".91 - .40 


39".35 


+ .04 


89".41 


+ .10 


39".02 -.29 


39 .32 +.01 


39 .25 


— .06 


39 .40 


+ .09 


39 .01 - .30 


38 .93 - .38 


39 .14 


— .17 


39 .36 


+ .05 


38 .86 - .45 


39 .31 .00 


39 .47 


+ .16 


39 .20 


-.11 


39 .51 + .20 


39 .17 — .14 


39 .29 


-.02 


39 .42 


+ .11 


89 .21 - .10 


39 .04 — .27 


39 .32 


+ .01 


39 .30 


-.01 


39 .17 - .14 


39 .57 +.26 


39 .40 


+ .09 


39 .41 


+ .10 


39 .60 +.29 


39 .46 +.15 


39 .33 


+ .02 


39 .43 


+ .12 


39 .54 +.23 


39 .30 - .01 


39 .28 


-.03 


89 .43 


+ .12 


39 .45 + -14 


39 .03 - .28 


39 .62 


+ .31 


39 .36 


+ .05 


39 .72 +.41 



From these the most probable value of the diameter is 
found to be 



B 



89".808 ± 0".022 



the probable error of a single observation being r = 0".136. 
Compare the theoretical and actual distribution of errors 



between 


0".00 


and 


0".05 


u 


.05 


(( 


.10 


<( 


.10 


(( 


.20 


it 


.20 


a 


.30 


(C 


.30 


u 


.40 


over 


.40 







In the following table the first column gives the successive 
values of the limiting error a, the second column the values 

of -, and the third column the corresponding values of P. 
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The fourth column contains the differences between the suc- 
cessive values of P, and by multiplying each of these 
differences by 40, the number of observations, we have the 
quantities in column five, which are the numbers of errors 
that according to the theory should fall within the corre- 
sponding limits. Column six shows the actual number of 
residuals occurring between these limits. 



a 


a 
r 


P 


d 


n 


n' 


0.00 


0.000 


0.000 














0.196 


8 


9 


0.05 


0.368 


0.196 














0.184 


7 


6 


0.10 


0.735 


0.380 














0.299 


12 


12 


0.20 


1.471 


0.679 














0.184 


7 


8 


0.30 


2.206 


0.863 














0.090 


4 


3 


0.40 


2.942 


0.953 














0.047 


2 


2 


00 


oo 


1.000 









This is a close agreement considering that the number of 
observations is not very large. Also the number of errors 
greater than the probable error should be equal to the num- 
ber less than it. On counting the residuals we find twenty-one 
less than 0".136 and nineteen greater. 



THE REJECTION OF OBSERVATIONS. 

92. After a series of measurements have been made, it is 
frequently found that one or two of the observations differ 
widely from the others, and hence it becomes a matter of 
great importance to establish, if possible, some criterion by 
which we may determine whether such discordant observa- 
tions should be rejected or not. We are not concerned here 
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with the question of the detection of a mistake or constant 
error, which a consideration of the circumstances of the 
observations or of the instruments might reveal, but it is 
assumed that there is nothing whatever to guide us except 
the mere fact of the unusual size of the residuals of the 
observations under discussion. To reject an observation 
merely because it differs considerably from the others is 
entirely unjustifiable, while to retain it without any investi- 
gation is a neglect of the evidence furnished by the observa- 
tions themselves. 

The adoption of any rigid criterion based upon the magni- 
tude of the residuals is perhaps more satisfactory from a 
mathematical standpoint than from that of a practical observer, 
and some of the latter are of the opinion that no observation 
should be rejected entirely, even the most widely discordant 
ones being given a certain weight. In this latter case, 
however, the Theory of Probability will furnish a guide as to 
the proper weights to assign to the different observations. 

Of the various criteria that have been proposed, that 
developed by Pierce (see Chauvenet, page 558) is perhaps 
the most complete. The derivation and application of this 
criterion is, however, somewhat long and complicated, and 
for all ordinary cases the following simple methods will give 
practically as good results. 

93. Criterion for the Rejection of a Single Doubtful 
Observation. It was shown in (83) that in a series of n 
observations the number of errors numerically less than a 
should be nP, and therefore the number of errors greater 
than a should be 

n - nP = n(l - P) (84) 

If the value of the expression in (84) is less than one-half, 
the occurrence of an error of magnitude a will have a 
greater probability against it than for it, and hence the obser- 
vation corresponding may be rejected. 
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Accordingly the limit of rejection, a, of a single doubtful 
observation is obtained from the equation 



n{l - P) = 



or 



P = 



2n - \ 
2n 



(85) 



94. Example. Fifteen observations on the value of an 
angle are made. Ought any of the observations to be 
rejected ? 



M 


V 


v^ 


v' 


«'2 


V 


" 


J," 2 


2° 23'.90 





.80 


.090 


_ 


.41 


.168 





.33 


.109 


23.76 


— 


.44 


.194 


— 


.55 


.303 


— 


.47 


.221 


25.21 


+ 


1.01 


1.020 


+ 


.90 


.810 








24.68 


+ 


.48 


.230 


+ 


.37 


.137 


+ 


.45 


.203 


23.96 


— 


.24 


.058 




.35 


.123 




.27 


.073 


24.26 


+ 


.06 


.004 


— 


.05 


.003 


+ 


.03 


.001 


24.83 


+ 


.63 


.397 


+ 


.52 


.270 


+ 


.60 


.360 


24.07 




.13 


.017 




.24 


.058 




.16 


.026 


23.98 


— 


.22 


.048 


— 


.33 


.109 





.25 


.063 


24.14 


— 


.06 


.004 


— 


.17 


.029 


— 


.09 


.008 


24.40 


+ 


.20 


.040 


+ 


.09 


.008 


+ 


.17 


.029 


24.38 


+ 


.18 


.032 


+ 


.07 


.005 


+ 


.15 


.023 


24.59 


+ 


.39 


.152 




.28 


.078 


+ 


.36 


.130 


24.10 




.10 


.010 


+ 


.21 


.044 




.13 


.017 


22.80 


— 


1.40 


1.960 














2° 24 .20 






4.256 






2.145 






1.263 


M, 




%v' 






%v'^ 






S„"2 
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Using all the observations we find 



Mo = 2° 24'.20 r = ,6745 i/i:?^ = .37 



14 



By (85), P = — = .967 

30 



By Table I, - = 3.17 .-. a = 1.17 

As the residual — 1.40 is larger than a, we reject the 
last observation. 

From the remaining observations we now compute a new 
mean value and a new set of residuals. And we find 



./ 2.145 



M'o = 2° 24'.31 r' = .6745 1/^—- = ,27 

^ 13 

By (85), P = ^ = ,964 

By Table I, - = 3.11 ,-. a = .84 

The third observation may accordingly be rejected. 
From the thirteen observations that remain we find 



M"o = 2° 24'.2.3 r" = .6745 Jh^ = .22 



P = — = .962 
26 



- = 3.08 ,-. a = 
r 



Therefore no more observations are to be rejected. 
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95. The Huge Error. In cases where the number of 

observations is not unusually large, a simple and safe criterion 
for the rejection of a doubtful observation is found in the use 
of the " Huge Error." 

This is an error of such a magnitude that 999 out of every 
1000 errors are less than it and only 1 as large as or greater 
than it. 

Therefore the probability that the error of any given 
observation will be less than the " Huge Error " is .999, and 
from Table I, 

when P = .999, — = 4.9 

T 

a := Huge Error = 4.9 r 

= 3.3 [»• (86) 

:= 4.1 a.d. 

Then in any limited series of observations, if an error 
greater than the huge error is found, we should reject the 
observation corresponding. 

See also, Jlolman, page 30 ; Wright, page 131. 



CONSTANT ERRORS. 

96. Throughout our discussion of the methods of adjusting 
observations so as to obtain from them the most probable 
values of the unknown quantities, all constant errors are 
supposed to have been eliminated before the Method of 
Least Squares is applied in deducing the results. 

If this is not done, and each observation is subject to the 
same constant error, the final result will be affected by an 
equal amount, and in short, the Method of Least Squares is 
not capable of removing or reducing the effect of errors of 
this kind. All that is accomplished by the use of the method 
is to reduce to a minimum the effect of the Accidental Errors. 
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Hence it will be seen that although by increasing the 
number of observations of a given kind we may increase the 
precision, that is, reduce the probable error, of our final 
result as much as we choose, yet we do not in this way 
necessarily increase the accuracy of the determination. 

But if the unknowns can be determined in several ways, or 
under a variety of different circumstances, with various 
instruments, or by different observers, then it is most proba- 
ble that the constant errors of the different sets of measure- 
ments will be grouped about the true values of the unknowns 
according to the exponential law of error. Accordingly a 
combination of such observations will enable us to increase 
not only the precision, but also the accuracy of the final 
result, the constant errors of the different sets tending to 
cancel each other in the same way that the accidental errors 
of a single set do. 

It is for this reason that determinations of a quantity from 
observations made in a variety of ways are more valuable 
than those obtained merely from different sets of measure- 
ments of the same kind. 

97. The probability of the existence of a constant error 
may often be expressed in the following manner. 

Example. A standard 100 ohm coil is compared with 
a Wheatstone's bridge and the mean result found to be 
100.90 ± 0.20. To find the probability that thfere is an error 
in the bridge between -|- 0.30 and -\- 1.50 ohms. 

Suppose the result 100.90 ± 0.20 is treated as a single 
observation, and we find by an application of (83) the prob- 
ability that the error of this observation is numerically less 
than 0.60 ohms. 

Here a _ .60 ^ ^^^ ^ _ ^^^ 

»• .20 

Hence, as far as is shown by the observations, the proba- 
bility that 100.90 ohms is within 0.60 ohms of the true 
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value is .957. But since it is known that the true resistance 
is 100 ohms, it follows that there is the same probability 
that there is a constant error in the bridge between -|- 0.30 
and -\- 1.50 ohms. 

98. Combination of Determinations having Different Con- 
stant Errors. In case two or more determinations of a 
quantity, together with their probable errors, are obtained, 
the method of combining them so as to secure the best final 
result was considered in paragraph 56, and in Example C, 
paragraph 57. But it was there assumed that all the results 
were subject to the same constant errors, while if this is not 
true the probable errors of the separate determinations bear 
no relation to their weights, and accordingly in such cases 
another process must be adopted. 

To determine whether the different measurements may 
fairly be considered to have the same constant errors we may 
proceed as follows : — 

Let the determinations of the quantity M be 

JTi ± ri (a) 

-3^2 ± r^ (b) 

and let the difference between these results be 

d = Ml — Mi (c) 

Then the probable error of d is by (66), 



R = Vi' + '•»' (d) 

If d is of such a magnitude that an accidental error as 
great as it may reasonably be expected, we may assume that 
the constant errors of (a) and (b) are the same, and pro- 
ceed as in Example C, paragraph 57. 

But if the probability of making two determinations which 
differ by the amount d is very small, we had best consider 
Ml and M^ to have the same weight, provided there are no 
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special reasons for regarding one better than the other. The 
final value of Jf will then be the arithmetical mean of Mi 
and M2, and its probable error will be found by (53). 

99. Example A. An angle is measured by a theodolite 
and by a transit with results 

By Theodolite, 24° 13' 36".0 ± 3".l 
By Transit, 24° 13' 24" ± 14" 

What is the most probable value of the angle and its prob- 
able error? 

Referring to the preceding paragraph, 

ri = 3.1, r^ =14, d z= 12, 

and the probable error of d is 



R = v'3.1^ + 14^ = 14 

Then from Table I the probability that the accidental 
error of a determination will be at least as large as 12 is 
found from 

« = 1^ = .86 ... 1 - P = .57 
r 14 

That is, there is more than an even chance that two such 
determinations of the angle will differ by as much as 12. 

Hence it is fair to assume that the two determinations are 
not affected by constant errors of different magnitudes, and 
they would be combined as in Example C, paragraph 57. 

Example S. Suppose the zenith distance, i!^ of a star, 
observed at two different culminations, is found to be 

My_ = 14° 53' 12".10 ± 0".30 
M^ = 14° 53' 14".30 ± 0".50 

What is the best final value ? 



Here d = 2.2, B = V -09 + .25 = .58 

a 2.2 

and for - = — = 3.8, 1 _ P = .01 
»• .58 
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Therefore the chance that the difference in the two deter- 
minations, due to accidental errors, will be as large as 2.2 is 
only one in a hundred. It is to be concluded then that the 
constant errors of observations at the two culminations differ 
by about 2.2, and as there is nothing to show that one meas- 
urement is more accurate than the other we will give them 
both the same weight and take the mean. Then the best 
value for the zenith distance is 

jr„ — 14° 53' 13".20 ± 0".74 

12".10 + 14".30 
For JJf„ = 14° 53' + ^ 

= 14° 53' 13".20 



/ 2.42 
By (53) r„ = -6745 y^-^^ = .74 

For a more extended treatment of this subject see Johnson, 
" The Theory of Errors and Method of Least Squares," 
chap. vii. 

THE WEIGHTING OF OBSERVATIONS. 

100. In case the relative worth of observations is not 
settled by methods already discussed, the proper weight to 
assign to each quantity in the final adjustment can only be 
determined from a full knowledge of all the circumstances of 
the measurements. Even then considerable experience in 
the particular work in hand is required before the best values 
for these weights can be assigned. The weight given to a 
quantity should never be considered final, but always subject 
to revision whenever new information with regard to the 
quantity is obtained. Thus an observation which at first is 
supposed to deserve a high degree of confidence is often 
found on later investigation to possess very little value, and 
vice versa. 

See also Wright, page 118. 
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OTHER LAWS OF ERROR. 

101. Although in the great majority of cases the distribu- 
tion of errors follows the exponential law thus far considered, 
there are a few special cases in which some of the suppositions 
made in deriving that law do not hold, and hence for the 
adjustment of such observations the corresponding special 
laws of error must be determined. 

For instance, in applying the exponential law we assume a 
large number of observations, that each observation is sub- 
ject to the same law of error, that small errors are more 
likely to occur than large ones, and that positive and negative 
errors are equally probable. Now it is easy to conceive of 
cases where only positive errors can occur, or where the 
probability of the occurrence of a small error may not be 
greater than that of a larger one, etc. If we can determine 
the different sources of error in any case and the relative 
effect of each upon the quantity sought, we shall arrive at 
the law of error for that particular set of observations. The 
case of most common occurrence is the following. 

102. Suppose all errors between the limits a and — a 
are equally probable, and that there are no errors beyond 
these limits. Then li y = <l>(x) is the equation of the 
Curve of Error, and its area is represented as in para- 
graph 18, we have 

X<t>{x) dx = 1 (a) 

a 

or 2 <^(a;) C"dx = 1 (b) 

since by the supposition made <j){x) must be a constant. 
Integrating and solving for </>(»;), we have 

y = ^(oc) = ~ (-87) 

2 <5t 
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To find the Mean Error we have by definition as in para- 
graph 52 

y? = f x^tj, (x) dx 

= 1 ("x'dx 
a Jo 

_ a^ 

3 



V- = 



a 



The Probable Error is derived from the equation 

Xr 1 

4>(x) ax = — 
»■ 2 



_ I dx = — 

a Jo 2 



r = 



a 



(88) 



(89) 



Finally, for the Average Deviation we have 
a.d. = ( X tj} (a) dx 

J ~a 

= _ j X dx 

a'^0 

a.a. = — 
2 

And the Curve of Error has the form 
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That the average deviation and probable error are in this 
case equal to one half of a may also be seen from the defini- 
tion of these quantities. 

Example. In taking a logarithm from a four place ' table, 
what is the probable error of the mantissa ? 

In this case the maximum error is .00005, and all 
errors between — .00005 and .00005 are equally probable. 
Therefore 

r = .000025 

103. The only other special case of common occurrence is 
that in which the error of a quantity is due to two sources, 
each of which can with the same probability assume all 
values between a and — a. Here it may be shown that 
the curve of error consists of two straight lines whose equa- 
tions are 

V = and y = ■ (91) 

2 

Also ^^ = -a\ r = (2-\r2)a (92) 

For a more extended discussion of special laws of error, see 
Wright, paragraphs 31 to 39. See also Example 151. 

104. Contradictory Observations. Suppose three observa- 
tions of a quantity give results 55, 56, 91. It is obvious that 
to take the arithmetical mean as the most probable value in 
this case would be contrary to the evidence furnished by the 
measurements, while in so small a number of observations it 
would not be allowable to reject the value 91 entirely. With 
such a series of observations no satisfactory solution can be 
obtained, but probably the best thing to do would be to take 
the value which has as many observations less than it as it 
has greater, or 56. 



CHAPTER VI. 

GAUSS'S METHOD OP SUBSTITUTION. 

105. The most laborious part of the application of the 
Method of Least Squares to the adjustment of observations 
consists in the formation and solution of the "Normal Equa- 
tions," and this labor increases enormously with increase in 
the number of observations, of unknowns and of conditions. 
It is not at all unusual to find that the adjustment of a single 
set of observations takes several weeks, even with all the aid 
that can be obtained from tables of logarithms, of squares, of 
products, and of reciprocals, and also from the use of calcu- 
lating machines. In general the computations can be per- 
formed more rapidly and with less fatigue by using a machine 
or a table of products than by using logarithms, but a combi- 
nation of methods is often desirable. 

106. In any case however it is of the utmost importance 
that the formation and solution of the normal equations 
should be effected in a systematic way, and that as far as 
possible checks on the numerical work be carried along in 
the computations. By a slight amount of additional work 
checks upon the results at successive stages of the solution 
may be obtained by means indicated in the following dem- 
onstrations, while the most satisfactory form for the solution 
of the normal equations is given by the "Method of Substitu- 
tion" proposed by Gauss. 

This method will now be explained; and it is to be observed 
that it is customary in this subject to enclose a quantity in 
brackets when the sum of a number of quantities of the same 
kind is to be denoted. Thus [a5] means the same as Sa6. 

For the sake of simplicity in demonstration it will be 
assumed that the observations are all reduced to weight unity. 
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107. Checks on the Formation of the Normal Equations. 
If, as in paragraph 76, we take for 



OBSERVATION EQUATIONS 



ajSi + JjSj, + 



• • fi'a^ff + mj = 






(A) 



«»«l + *«22 + • ■ • 9'»^« + n^n = Vn 

we shall have for 

NORMAL EQUATIONS 

[aa] Si + [aZ>] s^ -[- . . . [ag-] «g -f" [«»*] = 



Let 



[«2] ^i + [^2] 22 + 

«i + *i + • 

Oz + ^2 + • 



, . Ibq] z^ + [Jm] = 

( 

■ • C??] 2a + [?»»] = 

• qi -{- mi = Si 
. g'2 -|- ma = Sg 

• 3'n + »»» = «« 



(B) 



(C) 



(D) 



Multiplying the first of (C) by m^, the second by jWj, ... 
and adding, there results 

[a«i] + [6»*] + • • • [gwi] + [*»**»*] = ism] (93) 

Next multiplying each of equations (C^ by its a and 
adding, and then each by its b and adding, etc., we have 
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[««] + [a6] + . . . \_aq] + Vain] = [as] 
[aft] + [66] + . . . [6gJ + [6m] = [6s] 

(94) 

[ag] + [60-] + . . . \_qq-\ + [^m] = [<?s] 

Equation (93) will be satisfied if the absolute terms in 
the normal equations are correct, and equations (94) when 
the coeflBcients of the unknown quantities are correct. These 
check the formation of the normal equations. 

108. The Reduced ITormal Equations and the Elimination 
Equations. 

The value of e-i in terms of the remaining unknowns, 
derived from the first of equations (B), is 

, _ [«*] , [^ „ [«"^] .T^N 

^1 — ~ laa] ^"^ ~ \_aa-\ ^^ " • • • [aa] y^) 

Substituting this in the remaining n — \ equations, they 
become 

\w - [«*] [Sj} ^^ + • • ■ jc*'«3 - c«*J i3| = •^ 

\V>0-] - [_ac] f^l % + . . . {[cm] - [ac] M| = o 

jc*?] - [«?] gy[ -^ + • • • j[^H - [«?] ^1 = 

And letting 

[,m, 1] = [,^] - L^^ 



'f 



100 METHOD OF LEAST SQUABES. 

the above equations take the following form, which, being 
the same as that of the original normal equations, they are 
called the 

FIRST REDUCED NORMAL EQUATIONS. 

[5c, 1] Za + [cc, 1] Zs + . . . [eg-, 1] s, + [cm, 1] = 

(G) 

Ibq, 1] «a + [cq, 1] z, + . . . [qq, 1] Zg + [^m, 1] = 

An inspection of equations (F) will render it easy to form 
a rule for writing out any one of them. 

Now by means of the first of equations (Gr), eliminating 
S2 from each of the others in the same way that Zi was 
eliminated from the normal equations, there results the 

SECOND REDUCED NORMAL EQUATIONS 
[cc, 2] s, + . . . [eg-, 2] Zg + [cm, 2] = 



(H) 



[cq, 2] z, + . . . [qq, 2] z^ + [qm, 2] = 
In which 

fee 21 — rcc n [^c, 13 IK 1] 

[cc,2] _ [cc,l] ^^^-^ 



(I) 



[qm,22 = [,m,l] - ^-^^^^aSP 

Continuing this process we shall finally arrive at the single 
equation 

[y?. q-1] »« + [?»». q-1] = (J) 



from which the value of Zg is determined. 






/ > * ' .A 



' ,/ 
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The value of Zq_i will then be obtained by substituting the 
numerical value of Zg in the first of the preceding set of 
equations, and so on, until finally % is obtained from the 
first of the original normal equations. The equations from 
which the unknowns are actually determined are then the 
following, called the 

ELIMINATION EQUATIONS. 

[aa] »i + [a&] Sz + • • • [<*Q'3 «? + [am'] = 
[6&, 1] «2 + . . . \_bq, 1] «s + [6m, 1] = 

(95) 

Iqq, q-1] «, + Iqm, q-1] = 

It may be seen from the rule in paragraph 81 that 
Iqq, q— 1] is the weight of Zg, and the weight of any 
unknown might be found at the game time as its value by 
making it the last in the order of elimination, but except in 
special cases the weights had best be obtained by the general 
process of paragraph 115. 

109. Check on the Solution of the Normal Equations. 
Multiplying the first of the Observation Equations (A) by 
nil, the second by m^, . . . and adding the results, we have 

Imv'] = [am'] z^ -f- [Sm] z^ -\- ... [qm] Zg -|- [mm] 

But in equation (T), paragraph 82, it was shown that 
\m,v\ = [■««]. Therefore 

\yv\ = [am] Sj -|- [5m] z^ -\- . . . \_qm] Zg -\- [mm'] 

Substituting in this the value of z^ from the first of (95), 
we get the result 

[vv] = [hm, 1] Za + [cm, 1] z, -)- . . . [qm, 1] Zg -f- [mm, 1] 
in which 
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[aS] [^am] 



\bm, 1] = [Jm] — 



[mm, 1] = [wm] — 



[«a] 



[aa] 



(K) 



being similar in form to equations (F). 

Next eliminating z^ in a like manner, we get 

[yv'] = [cm, 2] Zs -|- . . . [qm, 2] e^ -j- [mm, 2] 

and continuing this process it finally appears that 

[u v] = [mm, q] [96] 

110. Arrangement of the Computations. In computing 
the coefficients that appear in the "Auxiliary" or "Reduced 
Normal Equations" it is most convenient to arrange the work 
in tabular form. The arrangement of the solution will be 
illustrated for an example containing four unknowns but it 
will be evident that the process can be extended to cover any 
case. 

Let the Observation and Normal Equations be represented 
by equations (A) and (B), there being only four unknowns 
Zx, Zj, Sj, Zi, and arrange a table as on the next page. In this 
scheme of solution the upper lines of the rows in the first 
compartment contain all the quantities that appear in the 
Normal Equations, together with [»n»i] and the quantities 
in the column headed s which are used in checking the 
results in accordance with equations (93) and (94). The 
other compartments contain the corresponding quantities 
for the Reduced Normal Equations, and the first line in 
each compartment gives the coefficients in the Elimination 
Equations. 



GAUSS'S METHOD OF SUBSTITUTION. 



103 



Scheme A. — Solution of the Normal Equations. 





a 


h 


e 


d 


m 


« 


a 


Iog[oo] 


[ai] 
log[ai] 


log [oc] 


[ad] 
log[ai] 


[am] 
log [am] 


Iog[oj] 


b 


log A, 


[JS] 

Ailab] 

log AM 


^»M 
log^jM 


[6rf] 

A,[ad] 

log A,[ad] 


[Jm] 

At [am] - 

log ^j [am] 


[fcj' 

A[<«] 
logA['"] 


c 


log J, 




M 

AM 

log AM 


[.d] 

A, [ad] 

logAEod] 


A, [am]. 
log ^g [am] 


A[»-] 
logA[<".l 


d 


logAa 




[dd] 

Ai[ad] 

logAC"-'] 


[rfm] 
- Aft lam] 
log A^ [am] 


[A] 

A[«»] 

iogA[''»] 


m 


log^„ 




[mm] 

-4ffl[am] 

log An lam] 


;[m.] 
. A['"] 
logA[m] 


h 




Iog[M.l] 


[6«,1] 
log[J«,l] 


log[M,l] 


[fc.,1] 

log[Sm, 1] 


[fe, 1] 
log [4s, 1] 


e 


log-B. 




[™, 1] 

logB.[J<:,l] 


[cd.l] 
log B, [bd, I] 


[cm. 1] 

-Ba [*"■■!]■ 

log j5„ [bm, 1] 


[". 1] 

B»[S»,1] 

Iog£,[4«, 1] 


d 


logS^ 




[dd, 1] 

Ba [hd, 1] 

logS^[W,l] 


,[rf».,l] , 

B,,[bm,l] 

logfij[4».,l] 


[A,l] 

-Bi[J», 1] 

log Si [4s, 1] 


m 


\og.B^ 




' [»m,l] 
£» [fo», 1] 
logSmifon, 1] 


[ms, 1] 

A„[SM] 

log B„ [fo, 1] 


e 




log[«,2] 


[0*2] 
log[«rf,2] 


[c;r!,2] 
log [cm, 2] 


[«.2] 
log [cs, 2] 


d 


logCi 




[<«rf,2] 

C^[ed,i] 

logCi[crf,2] 


[dm,2] 

C^ [cm. 2] 

log Ci [cm, 2] 


[*, 2] 
• Oa [cs, 2] . 
log Ci [cs, 2] 


m 


logCi. 




[mm, 2] 

C„ [cm, 2] 

log On [cm, 2] 


[mi, 2] 
C„[cs,2] 
log C„ [cs, 2] 


d 




[dd,3] 
\oe[dd,3] 


[dm,3] 
log [rfm, 3] 


[*, 3] 
log [A, 8] 


m 


log/?. = log-«. 




[mm, Z] 

D„ [dm, 3] 

logZ)„[ifm,3] 


■ [ms, 3] 
i>™ [A, 3] 
iogi>„[&,3] 


m 


M = 


[mm, 4] 


[ms, 4] 
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The logarithms of the quantities in the first row of each 
compartment are also written in, and from these by proper 
subtractions are obtained the logarithms in the margin, where 

A — l^ A — l2:S} A — l^^"'] 

_ [he, 1] _ [bm, 1] 

" ~ Ibh, ij ••• -^^ — [-55, 1] 

(L) 

Now in each compartment adding the logarithms at the 
margin to each of the logarithms in the first row of that com- 
partment we obtain the corresponding logarithms written in 
the other rows. The numbers represented by these logarithms 
are next written above them, and if each of these quantities 
is then subtracted from the one above it the result will be 
the corresponding quantity in the compartment below. Some 
of the squares in each compartment are left vacant as the 
quantities belonging to them have already appeared above. 

111. Application of Checks. Also, by (93) and (94), in 
the first compartment the quantities in the first lines of the 
last column should be equal to the sum of all the quantities 
in the first lines of the corresponding rows plus the quantities 
similarly situated above the first terms of the rows. Similar 
checks will apply in each compartment; for if from the 
second of equations (94) we subtract the product of the first 
equation multiplied by Ab we have 

[bb, 1] + \bc, 1] H \bm, 1] = lbs, 1] (97) 

In the same manner we may show that a corresponding 
check holds throughout, so that finally we shall have 

[mm, 4] = [ms, 4] (98) 
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The last compartment of the table is added to give this 
final check and the value of [vv'] in accordance with equa- 
tion (96). 

If the multiplications and divisions are simple or if a table 
of squares or products or a computing machine is used the 
logarithms will of course be omitted from the scheme of 
solution. 

112. Example. In order to illustrate the systematic for- 
mation of the coefficients that appear in the Normal Equations 
as well as the solution of the latter by the above method we 
will take the 



OBSERVATION EQUATIONS 



— 2l + «2 + «3 + ^4 



Zj -)- Zj — Zj 

Zi — 2z, 



; + Z4 = 0.1 



0.1 


Zl 


— Zj — 2Z4 = 


0.3 


0.6 




Z2 - 2a -f 2* = 


— 0.1 


0.1 


Sl 


- Zs = 


0.4 



First form a table containing the coefficients in these equa- 
tions and also the sums s. As a first check the sum of the 
quantities in column s should be equal to the sum of all the 
quantities in all the other columns. 



I. 


COBPFICIBNTS IN 


THE Observation Equations. 


Ifo. 


a 


b 


c 


d 


m 


s 


1 


- 1 


1 


1 


1 


- .1 


1.9 


2 


1 


1 


— 1 


- 1 


- .6 


— .6 


3 


1 


— 2 


- 1 


1 


— .1 


- 1.1 


4 


1 


— 1 





- 2 


- .3 


- 2.3 


5 





1 


- 1 


1 


.1 


1.1 


6 


1 





- 1 





- .4 


- .4 


Sum. 


3 





- 3 





- 1.4 


- 1.4 
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From these we now compute the coefficients in the Normal 
Equations (B), and also the necessary quantities for the 
check equations (93) and (94). 

II. COEITICIENTS IN THE NoEMAI, EQUATIONS. 



ISTo. 


aa 


ah 


ac 


ad 


am 


as 


bb 


be 


bd 


bm 


1 


1 


-1 


-1 


-1 


.1 


-1.9 


1 


1 


1 


-.1 


2 


1 


1 


-1 


-1 


-.6 


- .6 


1 


-1 


-1 


-.6 


3 


1 


-2 


-1 


1 


-.1 


-1.1 


4 


2 


-2 


.2 


4 


1 


-1 





-2 


-.3 


-2.3 


1 





2 


.3 


5 














.0 


.0 


1 


-1 


1 


.1 


6 


1 





-1 





-.4 


- A 











.0 


Sum 


5 


-3 


-4 


-3 


-1.3 


-6.3 


8 


1 


1 


-.1 



Ifo. 


bs 


cc 


cd 


cm. 


cs 


dd 


dm 


ds 


m,m 


ms 


1 


1.9 


1 


, 1 


-.1 


1.9 


1 


-.1 


1.9 


.01 


-.19 


2 


- .6 


1 


1 


.6 


.6 


1 


.6 


.6 


.36 


.36 


3 


2.2 


1 


-1 


.1 


1.1 


1 


-.1 


-1.1 


.01 


.11 


4 


2.3 








.0 


.0 


4 


.6 


4.6 


.09 


.69 


5 


1.1 


. 1 


-1 


-.1 


-1.1 


1 


.1 


1.1 


.01 


.11 


6 


.0 


1 





.4 


.4 





.0 


.0 


.16 


.16 


Sum 


6.9 


5 





.9 


2.9 


8 


1.1 


7.1 


.64 


1.24 



If the coefficients in the Observation Equations are large, 
so that it becomes convenient to use logarithms, other tables 
corresponding to I and II would be formed to contain these 
logarithms. 
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Substituting these quantities now in the general tabular 
scheme of paragraph 110 we have the results on the following 
page. The work of the first compartment is performed 
without the use of logarithms, as the numbers are simple. 
The quantities in the last column should all be zero accord- 
ing to the check equations, what small differences there are 
being due to the rejection of figures beyond the third place 
in the decimals. The decimal points in logarithms to which 
correspond negative numbers have been replaced by the 
letter n. 

The demonstrations that have been made now enable us to 
see at once from an inspection of the results in this table 
that 

Z4 = 0.238 

P., = 1.6 
[w] = .007 

Therefore substituting in equations (73) and (74) we 
have 

1^^ = .047 r^ z^ .032 

If the two values of [w] obtained in the solution had 
differed at all we should have taken the mean of the two. 

113. If the Elimination Equations (95) are divided by 
[_aa'], [66,1], [cc, 2], [<?<?, 3], respectively, they become 

«i -\- Ab^i -{- Acffa + Adffi ■+■ Am = 

«2 + -Bc«s + J5d«4 + Sm = .aqv 

«s + C^«4 + Cm = ^ ^ 

«4 + 1>™ = 

And the solution for the unknowns can be effected most 
conveniently by arranging the computations in the manner 
illustrated on page 109. 
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Scheme A. Solution of the Normal Equations. 
a b c d m s 8 



a 


5 


-3 


-4 


-3 


-1.3 


-6.3 





b 

A, 


-.6 


8 
1.8 


1 
2.4 


1 

1.8 


- .1 

.78 


6.9 

3.78 





c 
Ac 


-.8 




5 
3.2 



2.4 


.9 
1.04 


2.9 
5.04 





d 
As 


-.6 




8 
1.8 


1.1 

.78 


7.1 

3.78 





m 
Am 


-.26 




.64 
.34 


1.24 
1.64 





b 




6.2 
0.7924 


-1.4 
0„1461 


- .8 
9„9031 


- .88 
9„9445 


3.12 

0.4924 





c 
-Be 


9„3537 




1.8 
.316 
9.4998 


-2.4 

.181 
9.2568 


- .14 
.199 

9.2982 


-2.14 

- .705 
9„8479 





d 


9„1107 




6.2 
.103 
9.0138 


.32 
.114 
9.0552 


3.32 
- .403 
9„6049 





m 

Bra 


9„1521 




.30 
.125 
9.0966 


- .40 

- .443 
9„6463 





c 




1.484 
0.1715 


-2.581 
0„4118 


- .339 
9„5302 


- 1.435 
0„1569 


1 


d 


0„2403 




6.097 

4.489 

0.6521 


.206 

.589 

9.7705 


3.373 

2.496 

0.3972 


1 


m 

Cm 


9„3587 




.175 
.077 

8.8889 


.043 

.828 

9.5156 


1 


d 




1.608 
0.2063 


- .383 
9„5832 


1.227 
0.0888 


2 


ni 

Bm 


9„3769 = log-z* 




.098 

.091 

8.9601 


- .285 

- .292 
9„4657 





m 


0.238 = Zt [w] = .007 


.007 


.007 






GAUSS'S METHOD OF SUBSTITUTION. 109 

Scheme B. — Solution of the Elimination Equations. 



--z>™ 


-Caz, 


- B, Za 


- A^ 

— AaZi 

— Ac 23 

- Ab Zj 


Zl 

logzi 


Zs 

logZj 


Za 
logZ2 


Zl 


logC,^ 
log-Bd 
log^d 


log -Bo 
log-4o 


log ^6 




log A^ Zj 


log-BoSa 
log ^C 28 


log dz^ 
log-Bd Z4 
log AaZi 







114. Filling out this table for the example just solved, we 
have 



.238 


.228 


.142 


.260 




.414 


.031 


.143 






.145 


.514 
.191 


.238 


.642 


.318 


1.108 


9.3769 


9.8075 


9.5024 




0„2403 


9„3537 


9„7782 




9J107 
9„7782 


9„9031 






9„2806 


9„1612 
9„7106 


9„6172 

8„4876 












9„1551 
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115. The Weights of the Unknowns. In order to deter- 
mine the precision measures of ai, Zj, and 8g, it would 
next be necessary to compute the weights of the latter quan- 
tities. The demonstration of the processes by which these 
weights may be found will not be taken up here, as the best 
method to adopt varies a good deal with the character of the 
example, but a statement of the results in the general form 
of solution will be given. 

By treating the Elimination Equations in a way similar to 
that used in deriving equation (E) of paragraph 76 from 
equations (B) of the same paragraph, we may show that 

2l + ^m + -Sm ai + ^^ a.^ -\- -Z>m aj =: 

Z, 4- ^„ + C^p, + JD^p, = (100) 

s. + c^ + JD^y, = a 

where the a's, ^'s, y's, are determined from the equations 

Aa + Baa^ + C<ja2 + a, = ^d + C<j A + /Sa = 
^<, + ^„ ai + aa = Bo +P^ = ^ (101) 

^6 + ai = Q + 78 = 

Then by an application of the principles of Rule 1, para- 
graph (77), it may be shown that 



1 



+ r/./._ 91 + 



i>.. - \_aa-] "^ Ihh, 1] -^ \_cc, 2] "^ Idd, 3] 
1 1 , _Jl__ , iSa" 



i>.. — [&&, 1] "T \_cc, 2] ^ Idd, 3] 

1 1 v ' 

^ = ice, 2] + [t?<?, 3] 

J 1 
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These equations can of course be extended to cover any 
number of unknown quantities, and tabular schemes for the 
computations of the a's, j8's, y's, . . . and of the weights 
Fv -^V -^^8' • • ■ <"*ii readily be arranged. 

For a general demonstration of these results, and also for 
a discussion of special methods of solution, consult — 

Johnson, "The Theory of Errors and Method of Least Squares," 

chap. ix. 
Wright, " Treatise on the Adjustment of Observations," chap. iv. 
Chauvenet, " Spherical and Practical Astronomy," pp. 530-649. 



THE METHOD OF CORRELATIVES. 

116. The method of adjusting " Conditioned Observa- 
tions " explained in paragraph 33 is perfectly general, but 
where there are many conditions to be satisfied the solution 
is apt to be very laborious. For the case that occurs most fre- 
quently in practice, in which the observations are direct and 
equal in number to the number of unknown quantities, the 
process of solution devised by Gauss and called the " Method 
of Correlatives " is the most convenient. This method is 
derived as follows : — 

Let q observations, Jf^ M^, . . . Mg, of the respective 
weights Pi, jOj, . . . Pq, be made directly upon the values 
of q unknown quantities, and let the most probable values 
of the unknowns be 

Zl = -S^l + "ll 23 = -3^2 + ^'2) . . .Zg = Mg-^ Vg. 

Where v^, v^, , . . Vg, are the most probable corrections to 
apply to the observed values as well as in this case the 
residuals of the observations. 

If the n' condition equations are not linear they may be 
reduced to that form by the method of paragraph 44, so that 
we may assume for our 



(A) 
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CONDITION EQUATIONS 

«i«i + aa'Oi + • • • «g«'g -\- m-L = 

biVi -|- &2W2 + • • • ^q'Oq -\- m^ = 
hVl + 4^2 -f . . . IqVg -\- m^- = 

In which the quantities mj, m^, . . . m„', would all be zero 
if the observations were exact. It is to be observed that the 
coeflScients a, b, ... I are not arranged in the same order in 
these equations as they are in the observation equations of 
paragraph 107. 

The values of Vi, v^, ... Vg, must be determined so as to 
satisfy the above equations and also by the principle of Least 
Squares, so as to make 

PiVi^ -\- P'iO^ -(-... PqOq = a minimum. 

Corresponding with a minimum value of this last we have 

PiVidv^ + PiV^dv^ + . . . PgVgdVg = (B) 

for all possible simultaneous values of dvi, dv^, . . . dvg ; that 
is, for all values which satisfy the equations. 



aidvi -\- a^dVi -|- . . . Ogdvg = 
bidvi -\- b^dv^ + • • • bgdvg = 



(C) 



lidvi -\- l^dVi + . . . IgdVg = 

obtained by differentiating equations (A). 

Therefore, denoting the first member of (B) by R and the 
first members of (C) by Si, /S^, ... S^', it will be nec- 
essary that 

Ji _ kiSj, - k^S^- ... K,8n- = (D) 

where ^1, k^, ... k^', are undetermined coefficients. 
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This last equation will be satisfied if the coefficient of each 
differential in it is made equal to zero, that is, if 

PlVi = ^iffli + lcj)i + ... K<h 



(103) 



PqVg — \a^ + k^\ 4- ... hn\ 



All that remains therefore is to find values of v^, v^, ... Vg 
and ^1, ^2, ... Jc^,, which will satisfy simultaneously equa- 
tions (A) and (103), and that this may be done is easily 
seen from the fact that we have the same number of equations 
as unknowns. 

Substituting the values of v^, v^, ... Vg from equations 
(103) inequations (A) we have the following : 



fci S — 7- + fe S ^^ + . . . &„' S — + mi = 



2> ' i> ' p 

ab , , -bb , , ^bl 

4- fcj S — + . . . A!„' S — 



, „ ab , , -bb , , ^bl , 

fel S — - + fcj S — - + . . . fe' S — - + Wta = 



(104) 



, „ al , , „ bl , , „ll , 

Aji S — +^e^^^ + . . . fc„H - + w„' = 



The solution of these equations gives at once the values of 
Tci, &2, . . . A„', which are called the "Correlatives" of the 
Condition Equations. The values of v-^, v^, . . . Vg are then 
found by substituting the values of the k's in equation (103). 

117. As equations (104) are of the same general form as 
a set of Normal Equations, Gausses Method of Substitution 
can be advantageously employed in the solution. 
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When there is but a single equation of condition the second 
members of equations (103^ reduce to their first terms, and 
equations (104) reduce to the single equation 



(105) 



A;iS 


aa 


+ wi 


= 


)f «i, v^, 




. Vq in 


(103) be( 


ai 






a^ 


^«- 


t>2 = 


i>2 
_ aa 


^ P 




aq 


^ P 


Vg 


zrr 


Pi 

aa 
^ P 


mi 



Vi = ==j-— nil, Vi = — mi 



(106) 



It is from these results that the rules in paragraph 35 are 
derived. 

118. Example. Suppose we have given the observations 



Ml = 


2.02, 


weight 3 


M, = 


4.13, 


« 2 


M, = 


2.52, 


" 5 


M, = 


2.67, 


" 7 


M, = 


2.84, 


" 4 



(a) 



and let the most probable values of the unknowns be repre- 
sented by 

Zi = ilfi + vi, z^ — M,-\-v^,...z,= M^-\- vs (b) 

Also suppose that the unknowns are subject to the conditions 

% + 22 + ^8 + ^4 + ^5 = 14.0 



= 1.5 



(c) 
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Then expressing these conditions in terms of the corrections 
by means of (a) and (b), we have the 



CONDITION EQUATIONS 

«i + ■"2 + ■"s + W4 + •ws + -18 = 

■yj - Vi - .04 = ^^^ 

Referring to paragraph 116, we see that in this example 

n' ■ = 2, jwi = 0.18, wij = — 0.04 

For the purpose of computing the coefficients in equations 
(104) we next arrange the following table. 



p 


a 


h 


aa 
~P 


ah 


bb 
~P 


3 
2 
5 
7 
4 


1 
1 
1 

1 
1 



1 

- 1 



1 

T 
1 

1 

1 

7 
1 
4 



1 



1 
7 





1 

¥ 


1 

7 










599 

420 


5 
14 


9 

14 
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Substituting these results in equations (104), we have 



599 5 

Ai + — Aj + .18 = 

420 ^14 ^ 



— hi A Aj - .04 = 

14 ~ 14 

Solving, Ai = — .1647 

ife, = .1537 

Then from equations (103) we get at once 



(f) 



vy, = — .0549, W4 = — .0455, 

Vi = - .0055, Vi = — .0412. (g) 

w, = _ .0329, 

And by substituting these results in (b) we can obtain the 
values of %, z^, s,, s^, Sj. 

The above is the solution of Example 70, page 127. 



EXAMPLES. 



1. An urn contains five black balls, three red balls and 
two white balls. If three balls are drawn from the urn what 
different combinations may result, and what is the probability 
of each ? 

2. In a single throw with a pair of dice what is the 

probability that neither ace nor doublets will appear ? 

3. Four cards are drawn from a pack. What is the 
probability of getting four aces ? Of getting one of each 
suit? 

4. From a lottery of thirty tickets, marked 1, 2, . . . 30, 

four tickets are drawn. What is the probability that num- 

2 

bers 1 and 15 will be among them ? 

^ 145 

5. Find the odds against the appearance of 7 or 11 in a 
single throw with a pair of dice. 7 : 2 

6. I toss up n coins. What is my chance of getting just 
one head ? 

7. In a single throw what are the relative chances of 
throwing 9 with two dice and with three dice ? 24 : 25 

8. From 2 n counters marked with consecutive numbers 
two are drawn. What are the odds against having an even 
sum ? n : n — 1 

9. In two trials with a single die what is the probability 
of throwing (a) an ace the first time only? (S) at least one 
ace? 

10. Find the probability of throwing doublets one or 

91 
more times in three trials with a pair of dice. 

216 
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11. Find the probability of throwing exactly three aces 

125 

in five trials with a single die. 

3888 

12. A certain stake is to be won by the first person who 
throws 5 with a die of twelve faces. What is the chance of 
the sixth person ? 

13. A and B play chess. A wins on the average two 

games out of three. What is A's chance of winning just 

80 
four games out of the first six ? 

^ 243 

14. A and B shoot alternately at a mark. A hits once in 
n times and B once in n — 1 times. Find their chances of 
first hit, and the odds in favor of B if A misses on his first 
shot. Even, n : n — 2 

15. In how many trials will it be a wager of 4 to 3 that 
double five will be thrown with a pair of dice ? 30 

16. Find the probability of throwing one and only one 

5 
ace in two trials with a single die. — 

^ 18 

17. If I have three tickets in a lottery of four prizes and 

41 
eight blanks, what is my chance of drawing a prize ? — 

55 

18. Find the probability of throwing at least four aces in 

six trials with a single die. 

23328 

19. On an average seven ships out of eight return to port. 

Find the chance that out of five ships expected at least three 

.,, ^ 16121 
will return. 

16384 

20. In a lottery containing a large number of tickets, 
where the prizes are to the blanks as 1 : 9, find the chance 

of drawing at least two prizes in five trials. 

100000 
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21. In a purse are ten coins, all nickels except one which 
is a five-dollar gold piece ; in another are ten coins, all 
nickels. Nine coins are taken from the first purse and 
placed in the second, and then nine coins are taken from the 
latter and placed in the former. If you now had your choice 
which purse would you take ? 

22. A and B engage in a game in which A's skill is to 
B's as 2:3. What is A's chance of winning at least two 
games out of five ? 

23. If A's skill at a certain game is double that of B, 
what are the odds against A's winning four games before B 
wins two ? 131 : 112 

24. A party of twenty-five take seats at a round table. 
What are the odds against any two specified persons sitting 
next to each other ? 

25. A has three shares in a lottery in which there are three 
prizes and six blanks. B has one share in another where 
there is but one prize and two blanks. What are their 
relative chances of getting a prize ? A : B = 16 : 7 

26. Expand through the terms involving h' and k", the 
expression 

^-T + (y + *)' 

When a: is 1 and y is -, does ~ -\- y' increase or 
. 2 X 

diminish when x and y begin to increase at the same rate ? 

27. Given f {x, y) = x^(a -\- y)', expand the expres- 
sion (x -\- hy {a -{- y -\- ky. 

28. Find the value of logio a -\- cos b, when a = 1001 
and b =: 0.1°. Give the result first to five places and then 
to seven places of significant figures, in each case without the 
aid of tables. 4.000433 

29. Transform to the new origin, (2, — 3, 1), the equa- 
tion, iB* + y= + a* — 4 a: + 6 2/ - 2 z — 11 = 0, the 
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axes remaining parallel to the original ones. The equations 
of transformation are x = 2 -\- x' , y = — 3-|-y', 
s = 1 + «'. a;'' + 2/2 _|_ g2 _ 25 

30. Find the minimum value of x^ -\- xy -\- y^ — ax — by. 

[ab - a^ ~ b^) 

3 

31. Find the values of SB, y and a that render a maximum 
or a minimum the function x^ -\- y^ -\- z^ -\- x — 2 z — xy. 

32. Find the values of x and y that render a maximum 
or a minimum the expression sin x -\- sin y -\- cos {x -\- y). 

88. Find the co-ordinates of a point the sum of the squares 

of whose distances from three given points, {x^, y^), (Xj, yj), 

, . . . (x, -4- x„-i- Xx) 

(*8j ye)) 18 a minimum. x = ^ ' = — ! ~ 

o 

34. Given the volume, a\ of a rectangular parallelopiped, 
find its shape when its surface is a minimum. 

35. Find the volume of the greatest rectangular parallelo- 
piped that can be inscribed in the ellipsoid 

x^ y^ z^ 1 ^"'^'^ 

-. + ;^ + -^ - -3-7? 

86. In the Physical Laboratory apparatus for illustrating 
the Estimation of Tenths, the reading of the micrometer head 
for a certain setting is computed to be 2.3088, which may be 
taken as exact. Setting the apparatus by the eye the follow- 
ing readings are obtained : — 

2.314 2.324 2.310 2.519 2.326 
2.820 2.302 2.318 2.305 

What are the accidental and real errors and the residuals ? 
Are there any constant errors or mistakes ? If there are 
constant errors are they of the first, second, or third class ? 
How do you tell ? 
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37. Are the following observations such as to call for an 
application of the Method of Least Squares in their adjust- 
ment? 

K-J-y— a-)-M=: 5 3a! — y — z — u =z 1 

a;— 2y-)-2z-(-M= — 1 a; — y-j-4z-(-6M = 9 

38. In the case of direct observations, what other quanti- 
ties besides the arithmetical mean might reasonably be 
assumed to give plausible values of the unknowns ? Why is 
the arithmetical mean preferred to these ? 

39. Find the most probable value of a quantity Jf from 
the observations 

216.58 



216.27 


216.16 


216.04 


216.19 


216.44 


.29 


.43 


215.99 


.39 


.51 


.33 


.09 


216.23 


.14 


215.94 



Also test the result by finding the sum of the residuals. 

M = 216.251 

40. In the determination of a certain wave length, Row- 
land made the following observations. Find the most prob- 
able wave length. 

4.524 4.515 4.518 4.507 4.501 4.485 4.517 4.493 4.505 
.500 .508 .511 .497 .502 .519 .504 .492 

4.5055 ± 0.0017 

41. Ten measurements of the density of a body gave 
results as follows. Find the most probable density. 

9.662 9.664 9.677 9.663 9.645 
.673 .659 .662 .680 .654 

9.6639 ± 0.0022 

42. In a triangulation of the U. S. Coast Survey an angle 
was measured twenty-four times with results 
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116' 



43' 44".45 


48".90 


47".40 


47".85 


51".75 


50 .55 


49 .20 


47 .75 


50 .60 


49 .00 


50 .95 


48 .85 


51 .05 


48 .45 


52 .35 



51".30 49".05 46".75 
51 .05 50 .55 49 .25 
51 .70 49 .25 53 .40 

Find the best value for the angle. 116° 43' 49".64 ± 0".28 
43. The following observations were made with a sextant 
in order to determine the latitude of a place : — 



43° 4' 46" 


7" 


59" 


52" 


47" 


36" 


24 


28 


39 


52 


15 


40 



What is the most probable latitude? 43° 4' 27" ± 5" 

44. Determine the quantity M from the observations 



M = 81.55 .41 .68 


.25 .27 


.77 .13 .86 


.10 .03 


Wts. 


25 25 16 


9 9 


9 4 1 


1 1 


45. An 


angle M is measured with the following results : 


M p 


M p 


M p 


M p 


M p 


45".00 5 


42".50 5 


27".50 3 


36".25 2 


45".00 2 


31 .25 4 


37 .50 3 


43 .33 3 


42 .50 3 


40 .83 3 


45 .00 3 


38 .33 3 


40 .63 4 


39 .17 3 





Find the most probable value of the angle and test the 
result by computing %pv. 39".78 ± 0".94 

46. In one hundred measurements of angles made in the 
primasry triangulation in Massachusetts by the Coast Survey 
there were found between 



6".0 


and 


5".0 


1 error. 


1".0 


and 0".0 


26 errors. 


5 .0 


t( 


4 .0 


2 errors. 


.0 


" - 1 .0 


26 " 


4 .0 


cc 


3 .0 


2 " 


-1 .0 


" - 2 .0 


17 " 


3 .0 


(C 


2 .0 


3 " 


-2 .0 


« - 3 .0 


8 » 


2 .0 


<( 


1 .0 


13 « 


-3 .0 


" — 4 .0 


2 « 



Plot these results as in Figure 1, page 8. 
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47. In sixty-six determinations of the velocity of light 
made at Washington the percentage of errors of different 
magnitudes was found to be as follows : — 



Over — 7.0 


.8 


Equal to — 7.0 


.8 


-6.0 


1.6 


-5.0 


2.0 


-4.0 


3.9 


— 3.0 


7.8 


-2.0 


12.8 


-1.0 


17.1 



Over 


7.0 


0.0% 


Equal to 


7.0 


1.6 




6.0 


3.5 




5.0 


4.2 




4.0 


5.8 




3.0 


8.6 




2.0 


12.4 




1.0 


17.1 



Plot the results and draw a smooth curve. 
48. Given the observations 



X — y + 22 = 3 


4a; + y + 4s = 21 


3a; + 2y — 5s = 5 


— a; -j- 3y -(- 33 = 14 



find the most probable values of a;, y, and s. 

a; = 2.470 ± .038 
49. O being the level of the sea and Pj, Pj, Pj, three 
points whose altitudes are to be determined, the following 
observations are made : — 



Pi above = 10 ft. 


Pj above P, = 9 ft. 


P, « Pi = 7 


Pi « P^ = i 


P» " = 18 





Find the most probable altitudes. P, ::= 8.50 ± .29 

50. The altitudes of A above 0, B above A, and S 
above are found by measurements to be respectively 
12.3, 14.1, and 27.0 feet. What is the most probable value 
of each of these differences in level ? A = 12.50 ± .17 

51. Measurements of the ordinates of points on a straight 
line corresponding to abscissas 4, 6, 8, 9, are made with 
results 5, 8, 10, 12. What is the most probable equation of 
the line in the form y = mx -\- bl b = — 0.29 
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52. Find the altitudes in Example 49 if the observations 
have weights 5, 3, 6, 2, 4, respectively. 

53. Solve the example in paragraph 31, giving the obser- 
vations the weights 25, 25, 4, 4, 4, 4, 4, 4, 1, respectively. 

Elevation of Pj = 320.25 

54. Find the most probable values of Zi, Sj, and e^ from 
the observations 

Zi = 552.10 wt. 16 »i — S2 = .75 wt. 1 

— S2 -|- gj = .15 " 9 Zi 4- S2 — Sj = 552.05 " 1 
Zs = 551.23 "4 zi - Zb = .70 " 1 

Sa = 551.30 " 4 

Sa = 551.2345 

55. In the triangulation of Lake Superior there were 
measured at station the angles 

F P = 62° 59' 40".33 wt. 5 



F B = 64 


11 


34 


.92 


(C 


7 


FOB = 100 


20 


29 


.12 


(C 


4 


P B = 37 


20 


49 


.55 


(C 


7 


B B = 86 


8 


55 


.86 


C( 


4 



Required the adjusted values of the angles. 

F O P = 40".28 ± 0".34 
56. In the TJ. S. Lake Survey the following angles were 
measured at station North Base : — 

(1) Crebassa — Middle 

(2) Middle — Quaquaming 

(3) Crebassa — Quaquaming 

(4) Quaquaming — South Base 

(5) Middle — South Base 

Find the adjusted values of the angles. 

(1) = 58".965; r = 0".28 



55° 


57' 


58" 


.68 


wt. 


3 


48 


49 


13 


.64 


(£ 


19 


.04 


47 


12 


.66 


CC 


17 


54 


38 


15 


.53 


ct 


13 


.03 


27 


28 


.99 


(( 


6 
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57. Adjust the following observations of differences in 
level : — 



Altitude of A 


401.3 


wt. 16 


C above £ 


72.5 


wt. 9 


A above B 


220.8 


" 16 


A '^ £ 


222.0 


" 1 


A " C 


150.2 


« 4 


Altitude of JB 


180.7 


" 1 



58. In "Conditioned Observations" can the number of 
observations required be less than the number of unknown 
quantities? Why must the number of conditions be less 
than the number of unknowns ? 

59. From the following measurements of the angles formed 
at the centre of a disk by four radial lines, find the most 
probable values of the angles. 



A = 104° 25' 13" 


C = 86° 33' 20" 


B = 98 13 47 


i> = 70 48 23 




A =■ 104° 25 



2".25 

Also solve giving the observations the weights 5, 2, 1, 4, 
respectively. 

60. Four observations on the angle ^ of a triangle gave 
a mean of 36° 25' 47", two observations on JB gave a 
mean of 90° 36' 28", and three on G gave 52° 57' 57". 
Adjust the triangle. A = 36° 25' 44".2 ; r = 7".7 

61. Five angles at a station are measured, and also their 
sum. The observed sum differs from the sum of the five 
observed parts by the amount d. What are the adjusted 
values of the angles ? 

62. The three angles of a spherical triangle are measured 
with results 

A = 46° 17' 38".32 B = 73° 35' 16".15 
= 60° 7' 5".16. 

Adjust the triangle, knowing that the spherical excess is 
2".475. A = 39".3; ^ = 1".6 
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63. At the station Pine Mountain the following angles 
were observed between surrounding stations : — 

Jocelyne — Deepwater 65° 11' 52".500 wt. 3 

Deep water — Deakyne 66 24 15 .553 " 3 

Deakyne — Burden 87 2 24 .703 " 3 

Burden— Jocelyne 141 21 21 .757 « 1 

Find the most probable values of the angles. 

64. Solve Examples 55 and 56 by the method of " Con- 
ditioned Observations." 

65. .^ is a station whose altitude is known to be 5240.1 
feet. JB and O are floats on a lake, and Z> is a signal point. 
From the following observations determine the most prob- 
able altitudes of £, O and J). 

C below A 720.1 wt. 3 JB below A 719.7 wt. 3 
Z> " ^ 200.3 "5 B <■'■ D 520.9 « 2 
C " D 520.4 " 2 

66. Given the following observations, subject to the con- 
dition Zi -f- ^2 = 8s, find the most probable values of 

Si, 32. and Zj. 

2zi — Za + z, = 3.0 2S2 - z, = 1.0 

2zi — 3s2 = — 4.5 Zi + 2z2 = 5.1 

Zi + z, = 3.8 

67. The chemical composition of a specimen was found 
by several observers to be as follows : — 

Pb = .52 Other substances := .09 Au and Ag = .39 
Ag = .27 -P6 and Ag = .78 Impurities = .10 

Au = .11 P6 and impurities = .62 Au = .12 

From these observations find the most probable composition 
of the specimen. 



23'" 


41».041 


wt. 


30 


42 


14 .875 


(C 


7 


47 


27 .713 


C( 


8 


23 


46 .816 


(C 


7 


5 


12 .929 


C£ 


5 



EXAMPLES. 127 

68. From the following observations what are the best 
values of the unknowns, supposing that y and z must be 
equal ? 

a; -|- 2/ = 5.2 wt. 4 y -|- z = 4.2 wt. I 
a; = 3.0 " 9 z = 2.0 " 4 

a: — z = 1.1 "1 

69. In determining the difference in longitude between 
various cities the results obtained were 

( 1 ) Cambridge — Washington 

(2) Cambridge — Cleveland 
(8) Cambridge — Columbus 

(4) Washington — Columbus 

(5) Cleveland — Columbus 

Adjust these observations. 

70. The capacity of a condenser is known to be 14.0 m. f. 
It is divided into five sections, a, b, c, d, e, and it is known 
that the difference between b and d is 1.5 m. f. Find the 
most probable capacities of the sections from the observa- 
tions 

a = 2.02 wt. 3 d = 2.67 wt. 7 

b = 4.13 "2 e = 2.84 " 4 

" == 2.52 "5 ^ ^ 19651 

71. If the unknowns in the following observations are 
subject to the condition a; -(- 2y -)- 3z =: 36, what are 
their adjusted values ? 

X = 4.3 wt. 1, y = 5.7 wt. 4, z = 7.3 wt. 9 

X = 3.77 

72. A cannon is discharged horizontally from the top of 
a bluff. Observations on the time, and distance of fall of the 
ball gave the results 

t = 0.5 1.0 1.5 2.0 seconds 
s = 1.2 4.0 9.1 15.0 metres 
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What curve, passing through the point of departure of the 
ball, will represent the above observations ? 

73. An Argand burner shows the following efficiencies 
with varying rates of gas consumption : 

g = 2.0 2.8 2.8 3.3 4.0 4.5 5.0 feet 
E = 2.1 2.4 2.5 3.0 3.2 3.8 4.1 

Find the equation of the straight line which best rep- 
resents the relation between g and JE. The measurements 
on g are without appreciable error. 

74. Observations are made upon the expansion of Amyl 
alcohol with change in temperature as follows : — 



F = 1.04 


1.12 


1.19 


1.24 


1.27 cu. cm. 


t — 13.9 


43.0 


67.8 


89.0 


99.2 C. degrees 



If V = \ -\- JB t -\- O t^ expresses the law connect- 
ing the volume and temperature, find the most probable 
values of S and G. 

75. In a Hooke's joint where the angle between the axes 
is 45°, X being the angular rotation of the driver, and y 
that of the follower, from the following measurements find 
the equation of a curve that will represent the relation 
between x and y — x. 



X 


y — X 


X 


y — X 


X 


y—i 


0" 


0°.0 


80° 


- 5°.8 


140° 


8°.8 


20° 


- 5°.7 


90° 


- 2°.0 


160° 


5°.3 


40° 


— 9°.9 


100° 


2°.3 


180° 


o°.o 


60° 


- 10°.4 


120° 


8°.0 







y - X = — 0.85 - 9.82 sin 2a! + 0.92 cos 2a! 

76. A series of observations extending over a period of 
thirty years was made by Quetelet to determine the daily 
variation in temperature at Brussels. The mean results of 



EXAMPLES. 129 

the measurements are given below. From them derive an 
equation to express the temperature at any time of the year. 



Jan. 


4°.66 


May 9°.83 


Sept. 8M6 


Feb. 


5°.42 


June 10°.09 


Oct. 6°.56 


Mar. 


6°.77 


July 9°.71 


Nov. 5°.10 


Apr. 


8°.59 


Aug. 9M4 


Dec. 4°.41 




y = 


7.369 + 0.9854 sin 


30x — 2.7084 cos 30a; 






+ 0.0100 sin 


60a: - 0.1950 cos 60a; 






— 0.0133 sin 


90a; +0.1783 cos 90a; 



In this answer the values of x begin at the 16th of Janu- 
ary, and represent the time in months. 

77. The law connecting the time of vibration of a pendu- 
lum with its length is assumed to be of the form, T =^ m i". 
From the following observations find the most probable values 
of m and n. 

T — 12.9 11.6 10.4 9.7 5.3 4.6 
i = 164.4 132.9 107.6 98.5 28.4 20.6 

L is in centimetres, T in tenths seconds. n = 0.5000 

m = 1.0044 

78. Determine the equation of a curve which will repre- 
sent the following observations : 



X 


y 


X 


y 


X 


y 


0.0 


0.00 


1.5 


1.09 


3.0 


8.65 


0.5 


0.04 


2.0 


2.56 


3.5 


13.72 


1.0 


0.31 


2.5 


4.99 


4.0 


20.47 



79. Determine the equation of a curve which will repre- 
sent the relation between x and y in the observations, 



X 


y 


X 


y 


X 


y * 


X 


y 


0.0 


4.51 


0.3 


4.09 


0.6 


3.03 


1.2 


0.92 


0.1 


4.44 


0.4 


3.76 


0.8 


2.24 


1.5 


0.38 


0.2 


4.31 


0.5 


3.42 


1.0 


1.49 


2.0 


0.05 
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80. At a station P the angles between a straight line 
passing through P parallel to the axis of JT and the direc- 
tions from P of four points Pi, P^, P3, P^, are measured. 
Having given the coordinates, (a, b), of the four points, 
find the coordinates of P. 



"aint. 


Coordinates. 


Angle. 




a 


b 




A 


4.21 


3.24 


39° 18' 


A 


1.21 


2.10 


147° 54' 


^. 


- 0.51 


0.22 


205° 24' 


P. 


2.50 


- 1.10 


277° 15' 



If the coordinates of the point P are (», y), and the 
angle is denoted by A, we have 

y — * 

tan A = 

X — a 

81. If aeinbx = M, and values of M are observed 
for known values of a and J, determine the most probable 
value of X. If xl is an approximate value of x found by 
trial, and m = a sin bx' — M, we shall have 

^ 'Sitt b m cos bx' 
S (« 6 cos 5a:')'' 

82. If in one series of observations the value of h is 
twice what it is in another, what is the relative probability of 
the occurrence of an error of given magnitude a in the two 
series ? Show what the curves of error will be in the two 
cases. What error has the same probability for its occurrence 
in each series ? What is the relative probability of the occur- 
rence of an error not greater than a in the first case and not 
greater than 2a in the second case ? 

83. From 64 observations the latitude of a station was 
found to be 49° 10' 9".110 ± 0".051. What was the prob- 
able error of a single observation ? 0".41 
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84. If twenty measurements of an angle give a result with 
an A.D. of 0".38, and it is required to find the angle so 
that the A.D. shall he only 0".25, how many more observa- 
tions must be made ? 27 

85. From the following determinations of the area of a 
field find the most probable area and its probable error. 

5674 ± 12, 5680 ± 4, 5685 ± 3, 5682 ± 1, 5678 ± 2 

A — 5681.41 ± 0.84 

86. From the following measurements by Fizeau and 
others, find the most probable value for the velocity of light 
together with its probable error. Measurements are in kilo- 
meters. 

298000 ± 1000 299990 ± 200 299930 ± 100 
298500 ± 1000 300100 ± 1000 

V = 299917 ± 88 

87. Two different instruments give for the value of an 
angle, 

84° 55' 33".0 ± 4".l, 34° 55' 36".0 ± 6".3 

What is the best value to take for the angle ? 

34° 55' 33".9 ± 3".4 

88. Determinations of the difference in longitude between 
Washington and Key West made on seven different days 
gave the results 

19"* 1'.42 ± 0'.044 19*" 1*.60 ± 0'.046 

1 .37 ± .037 1 .55 ± .045 

1 .38 ± .036 1 .57 ± .047 

1 .45 ± .036 

What is the best value and its probable error ? 

1'.460 ± o'.oie 
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89. In the triangulation of Lake Ontario two different 
instruments gave for an angle, 74° 25' 5".429 ± 0".29 from 
sixteen readings, and 74° 25' 4".611 ± 0".22 from twenty- 
four readings. Find the most probable value of the angle 
and its probable error. 

90. In each of Examples 39-45 find the mean and prob- 
able errors and average deviation of each observation and of 
the most probable value, using formulas from (50) to (63) 
according as they apply. 

91. In Example 42 divide the observations in their order 
into six groups of four observations each and compute the 
mean of each group. Then determine the probable error of 
the first of these means : ( 1 ) considered as a single measure 
of four times the weight of those in Example 42 ; (2) directly 
as one of six observations of equal weight; (3) as a deter- 
mination from its four constituents. 0".67 ; 0".72 ; 1".00 

92. The following twenty-nine measurements on the den- 
sity of the earth, made by Cavendish, give as a mean result 
5.48. What is the probable error of an observation ? Solve 
by the usual method and also by taking the residual that 
occupies the middle position. 0.14 



5.50 


5.55 


5.57 


5.34 


5.42 


5.30 


.61 


.36 


.53 


.79 


.47 


.75 


.88 


.29 


.62 


.10 


.63 


.68 


.07 


.58 


.29 


.27 


.34 


.85 


.26 


.65 


.44 


.39 


.46 





93. What is the probable error of the mean of two obser- 
vations which differ by the amount a ? 

94. A base-line is measured five times with a steel tape 
reading to hundredths of a foot, and five times with a chain 
reading to tenths of a foot, with results 

By tape, 741.17 741.09 741.22 741.12 741.10 
By chain, 741.2 741.4 741.0 741.3 741.1 



8.78 


-8.78 


-8.68 


.76 


.51 


.63 


.85 


.64 


.58 
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Find the probable errors and weights for a single observa- 
tion in each case, and also the adjusted length of the line 
and its probable error. 741.146 ± 0.015 

95. Twenty-one determinations of a chronometer correc- 
tion gave results 

8.80 - 8.96 - 8.83 - 8.79 
.75 .64 .70 .90 

.78 .65 .64 .93 

Find the probable error of the mean by using both formulas 
(53) and (57), and also determine the probable error of a 
single observation by taking the middle residual. 

0.017; 0.018; 0.09 

96. In the following observations show that Mo := 49.64, 
p. = 1.95, r = 1.31, /*„ = 0.40, »•„ = 0.27, /.ts = 0.87, 
»•, = 0.59. 

M = 48.81 48.76 49.53 51.56 50.38 49.84 
p = b 45 3 2 5 

97. Observations on the time of ending of a transit of 
Mercury are made by different observers with a variety of 
instruments and under more or less favorable circumstances. 
If the weights assigned by the computer are as indicated, find 
the best value for the time and its probable error. 



1* 38" 23' 


wt. 1 


38™ 26' 


wt. 3 38"" 19' wt. 8 


37 55 


" 


38 21 


" 2 38 21 « 2 


88 10 


« 1 


38 18 


" 2 38 15 « 2 
<„ = 5*38'"19'.9 



98. An angle is measured five times with a theodolite, and 
seven times with a transit, giving results 

Theodolite, 31".7, 39".8, 40".7, 28".6, 32".3 

Transit, 32 .8, 36 .7, 38 .2, 29 .3, 41 .6 85".3, 36".2 
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If the relative values of readings by the two instruments 
are as 3 to 2, what is the most probable value of the angle ? 
What is the mean error of the result ? 

99. Given M^ = 65.58 ± .59, M^ = 35.15 ± .93, 
M^ = 49.64 ± .27, find the probable errors of 4 ilfi - 

M W M 

3 Jf, + 2 JIfs and of ^ -\- :!^ - t^. 3.69 ; 0.43 

^ o 4 

100. The three angles of a triangle are measured, and the 
probable error of each observation is r . What is the prob- 
able error of the triangle error ? ^ i/^ 

101. The zenith distance of a star on the meridian is 
observed to be z = 21° 17' 20".3 ± 2".3. The declina- 
tion of the star is given as <^ = 19° 30' 14".8 ± 0".8. 
What is the latitude of the place and its probable error ? 

i = s + (? = 40° 47' 35".l ± 2"4. 

102. The zenith distance s of a star at upper culmina- 
tion is observed n times, and its zenith distance z' at 
lower culmination n' times. If the latitude is given by 
X, = 90° — i (2 + ^')> ^^^ t^^ probable error of an 
observation is r, what is the probable error of the latitude ? 

103. The horizontal force necessary to start a 100-pound 
weight sliding along a table is observed to be 15.5 ± 0.2 
pounds. Find the probable error of the coefficient of friction. 

104. If a line is measured by the continued application of 
a unit of measure, and r is the probable error of the placing 
and reading of this measure, what is the probable error of 
the length I ? ^ J^ 

105. If the average deviations of Sj, e^, Zg, are a, b, o, 
respectively, what is the average deviation of z,^ -\- z^ -j- z^ ? 

106. If the radius of a circle is measured with result 
1000.0 ± 2.0, how should the circumference and area be 
expressed ? 

107. Two sides, a and J, and the included angle G of 
a triangle are measured with results a = 252.52 ± .06 
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feet, b = 300.01 ± .06 feet, C — 42° 13' 00" ± 30". 
What is the area and its probable error ? 2b4b2 ± 9 

108. Measurements of adjacent sides of a rectangle gave 
a ± j-i, and b ± r^. What is the probable error of the 
area, and for what kind of a rectangle will this probable 
error be the least ? 

109. If the measured sides of a rectangle have the same 
a.d., what is the a.d. of the diagonal determined from 
them ? Same 

110. If the sides of a rectangle are measured in the 
manner indicated in Example 104 and found to be a and b, 
what is the probable error of the area ? 

111. The correction to be applied to a chronometer is 
found to be + 12"" 13'.2 ± 0*.3. Ten days later the cor- 
rection is again determined and found to be 12™ 21'.4 ± O'.S. 
What is the mean daily rate and its probable error ? 

0'.820 ± 0'.042 

112. Measurements of the compression of the earth's 
meridian have resulted in 

± .000046 

294 

What is the probable error of the denominator 294 ? 3.98 

113. The current flowing in a circuit is due to two 
sources whose electromotive forces are determined to be 
«! = 200 ±2, ej = 400 ± 3. The resistance of the 
circuit is 30 ± 1. Find the current and its probable error. 

20 ± 0.68 

114. The side b and angles B and C of a triangle are 
measured with results b = 106 ± .06 metres, £ z= 
29° 39' ±1', C = 120° 7' ± 2'. What is the most 
probable value of the angle A and of the side c ? 

A = 30° 14' ± 2'.2; c = 185.35 ± 0.15 

115. The distance between two divisions on a graduated 
scale is measured by a micrometer. Show that the average 
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deviation of the mean of two results is the same as the aver- 
age deviation of a single reading. 

116. If the weights of the determinations of three angles 
A, JB, C, are 3, 3, 1, respectively, what is the weight of 
the sum of the three angles ? 0.6 

117. If the weight of x is p, what is the weight of 
log„ X ? 

118. If a; = — and the weight of y is p, what is the 
weight of a; ? c^p 

119. In Example 107, how closely must the parts be 
measured in order to obtain the area within 0.5 per cent ? 

120. From observations on I and t the value of g is 
to be computed by the pendulum formula 



/ 



What changes in g will be produced by changes in I and 
t of 8i and Sj units, respectively, and what are the allow- 
able errors in I and i if ^ is to be determined within 
1 per cent ? 

121. The moment of inertia of a cylindrical bar is to be 
obtained from measurements on its mass m, its length h, 
and its diameter <?. The error in the determination of m 
is negligible, the precision of the determination of d is four 
times that of h. If the measurements give m = 48, 
h = 8.000, d = 1.200 ± 0.10, and 



I = ml I 

\l2 ^ 16/ 



what is the probable error of T, and what should be the 
ratio of 0? to A to determine T most accurately ? 

d: h = 256 : 9 
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122. If observations give for a certain quantity x the 
value 303, with a mean error of 2, what is the mean error 
of the expression 3 a; -|- logio 2 a;? 

123. The probable error of the determination of the 
angle A is 20". What is the maximum probable error of 
sin A -\- cos A ? 

124. If the probable error of an observation on an angle 
is 10", is there any difference between the probable error 
of the function sin A -(- cos A -\- sin C and of the func- 
tion sin A -\- cos Ji -|- sin C, supposing A and -B are 
of the same magnitude ? 

125. Given the observations, 

Zi - 2z2 + Ss - 3 = Bzi + Z2 + 2s3 - 17 = 
2zi + 3z2 — 4z8 — 2 = - Zi -f 4z2 + 3sj - 10 == 

Find the most probable values of Zi, z,, Zj, and also their 
weights and precision measures. 

Zi = 3.541 ; p„ = 29 ; r,, = .024 

126. Find the weights and precision measures of the 
unknowns in Examples 48 to 57. 

127. Determine the probable errors of the constants in 
Examples 72 to 79, inclusive. 

128. The length of a pendulum which beats seconds is 
given by 

I = r J^ (— q — s\l' sin^Z 

where l' is the length at the equator, q the ratio of 

289 

the centrifugal force at the equator to the weight, and s the 

compression of the meridian regarded as unknown. Putting 

I' = 991"- + X, /^q - s\l' = y, 

observations in different latitudes gave in millimetres the 
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following equations, from which we are to determine I and 
s together with their probable errors : — 

X + 0.969y = 5.13 x + 0.152y = 0.77 

X -\- 0.7492/ = 3.97 x + 0.327y — 1.70 

X -\- 0.426y = 2.24 x -f QMby = 8.62 

X -\- 0.0952/ = 0.56 x -j- 0.7932/ = 4.23 

X = 0.19 

r = 991.069 ± .026; s = — ± 0.00046 

294 

129. Find the weights and precision measures of the 
unknowns in Examples 64 to 70, and also in Examples 59 to 
63, and in 71. 

130. In Example 95 the probable error of a single obser- 
vation is 0.08 seconds. Find the number of errors which 
should fall between 0.00 seconds and 0.10 seconds, between 
0.10 seconds and 0.20 seconds, and also the number that 
should be over 0.20 seconds. Compare the results with the 
number actually found. 

131. In 470 determinations of the right ascensions of 
Sirius and Altair made by Bradley, the probable error of a 
single observation was 0".2637. The number of errors falling 
between specified limits was as shown below. Compare this 
result with the distribution of errors called for by the theory. 

Limits. Errors. Limits. Errors. 



0".0 to 0".l 


94 


0".6 to 0" 


'.7 


26 


.1 to .2 


88 


.7 to 


.8 


14 


.2 to .3 


78 


.8 to 


.9 


10 


.3 to .4 


58 


.9 to 1 


.0 


7 


.4 to .5 


51 


Over 1 


.0 


8 


.5 to .6 


36 









132. What is the probability that the error of a single 
observation will be as large as twice the probable error? 
As large as five times the probable error ? 
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133. On the average how many observations must be made 
before an error as large as three times the mean error will 
occur ? 

134. In Example 46, assuming that all errors between any 
two limits fall half way between those limits, compute the 
average deviation and mean error of an observation and com- 
pare their ratio with the theoretical value given in the table 
in paragraph 55. 

135. A line is measured 500 times and the probable error 
of each observation is 0.6 cm. How many errors should 
occur between 0.4 cm. and 0.8 cm. ? 

136. Show how the value of t could be determined 
experimentally from observations such as those in Example 
131. 

137. In a system of observations all equally good, r being 
the probable error of a single observation, if two observations 
are taken at random, what quantity is their difference as 
likely as not to exceed, and what is the probability that the 
difference will be less than r'i ^sT^', 0.367 

138. In the following measurements of an angle, ought 
any -of the observations to be rejected ? 

12' 51".75 47".85 47".40 48".90 44".45 
48 .45 51 .05 48 .85 50 .95 
50 .60 47 .75 49 .20 50 .55 

139. Determine whether any of the observations in 
Example 44 should be rejected. 

140. A quantity M is measured with the results given 
below. Ought all the observations to be retained ? 

M = 236, 251, 249, 252, 248, 254, 246, 257, 243, 274 

141. A certain angle has been laid out with such accu- 
racy that its true value may be taken as exactly 90°. Twenty- 
five observations are made upon it with a transit that it is 
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desired to test, and the result obtained is 89° 59' 57" ± 0".8. 
What are the odds in favor of a constant error in the instru- 
ment between — 1" and — 5"? Between 0" and — 6"? 

908:92; 86:1 

142. Repeated measurements of a standard metre bar 
with a decimetre scale gave a result 10.032 ± 0.010. What 
are the odds in favor of a constant error in the scale between 

0.010 , 0.054 . ,^ ^ 

and ? 43 : 7 

vTo v^lO 

143. Two determinations of the length of a line gave 
683.4 ± 0.3 and 684.9 ± 0.3, respectively. Show that the 
best value for the length is 684.15 ± 0.51, and that the 
probable systematic error of each determination is 0.65. 

144. Two men A and B observe an angle repeatedly 
with the same instrument with results 

A. B. 

47° 23' 40" 23' 35" 47° 23' 30" 24' 00" 

23 45 23 40 23 40 23 20 

23 30 23 50 

Is there any relative personal error, and what is the best 
final value ? 47° 23' 38".2 ± 1".6 

145. Three independent determinations of the capacity 
of a condenser made with three different instruments gave 
results 42.22 ± 0.21, 43.40 ± .15 and 44.20 ± 0.18. 
What is the most probable value of the capacity ? 

For extended treatment of the subject illustrated in Exam- 
ples 143 to 145 see Johnson., "The Theory of Errors and 
Method of Least Squares," chap. vii. 

146. In an estimation of tenths what is the probable error 
of an observation ? What is the average deviation ? 0.025 

147. In dbtaining the angle of deflection of the needle of 
a, tangent galvanometer by the usual method what is the 
probable error of the result ? 
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148. If all the errors of a series of observations must fall 
between and a, and the frequency of any error is pro- 
portional to its magnitude, what is the Curve of Error? 

What are the values of r, a.d., and /i ? >• = 

V^ 

149. In Example 148 what is the probability that the error 
of a single observation will be as large as 0.5 a. 

150. If all values of x between and a are possible, 
and their probabilities are proportional to their squares, find 
the mean value of x and the probability that x will be as 
large as 0.5 a. Also draw the Curve of Error. 

151. What is the greatest probable error of a logarithm 
found by interpolation in a seven place table ? .000000015 

152. Given the following set of Normal Equations, together 
with [mm] = 1.3409, find the most probable values of 
the unknowns and their weights and probable errors. There 
were sixteen observations. 

3.1217 Zi + .5756 z^ - .1565 z^ - .0067 s^ 

- 1.5710 = 
.5756 2i + 2.9375 z^ + .1108 z, — .0015 s^ 

+ .9275 = 

- .1565 Zi -1- .1103 Za 4- 4.1273 z, -f .2051 Z4 

-(- .0652 = 

- .0067 Zi - .0015 Z2 + .2051 Zj -j- 4.1328 z^ 

4- .0178 = 
zi = 0.583 ± 0.018 
Zi = - 0.004 ± 0.015 
p,, = 4.12 

153. From ten observation equations, for which was found 
[mrn] z=z 2.6322, there resulted the normal equations 

5.2485 zi - 1.7472 z^ - 2.1954 z^ + 0.5399 = 

- 1.7472 zi + 1.8859 z^ + 0.8041 z^ - 1.4493 = 

- 2.1954 zi 4- 0.8041 z^ + 4.0440 S3 - 1.8681 = 



142 METHOD OF LEAST 8QUABE8. 

Find the most probable values of gj, Zj and Sj together with 
their probable errors. Si = 0.42 ±0.11 

154. Find the most probable values of the unknowns in 
the normal equations 

459 Si — 308 Zi — 389 Sg + 244 S4 — 507 = 

- 308 Si + 464 Sj + 408 83 — 269 z^ + 695 = 

- 389 Si -|- 408 Sj! -I- 676 z^ - 331 Z4 -f 653 = 
244 Zi — 269 Sjj - 381 z^ + 469 z^ — 283 = 

[mm] = 1129 
s, = - 0.488 ; A. = 281 

155. If thirteen observation equations give rise to the 
result [mm] ^ 100.34 and to the normal equations 

17.50 Zj - 6.50 Zj - 6.50 z^ - 2.14 = 

- 6.50 Zi + 17.50 Zj - 6.50 Zj — 13.96 = 

- 6.50 Zi - 6.50 Zj + 20.50 z, + 5.40 = 

show that the most probable values of the unknowns are 

zi = 0.67 ± 0.60, Zj = 1.17 ± 0.60, z, = 0.32 ± 0.55 
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BUSBIENTS OF THE THEORY OF FROBABIUTT. 

200. Definition. If an event can happen in a ways, and. 
fail in h ways, and all these ways are equally likely to occur, 

the probability of the happening of the event is — ^^ . , 

and the probability of its failure is T~r"- 

Since the event must either happen or fail, the sum of the 
above probabilities must represent a certainty. But 

"■ 4, \ 1 

That is, the probability of a certainty is expressed by unity. 
Also, if the probability, J^, of the happening of an event is 
known, the probability of its failure is given at once by 
1- P. 

201. Example A. A single throw is made with a pair of 
dice. What is the probability that the siun of the spots 
turned up will be 5 ? 

Number of ways of throwing the dice is 6 X 6 :^ 36 

Number of ways of throwing five is 4 

4 1 

.•. Probability of throwing five is ■ — := — 

Example B. A coin is tossed up six times. Find the 
chance that three heads and three tails will be the result. 
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Number 6f ways of throwing the coin is 2' = 64 

6x5x4 

Number of ways of throwing three heads is = 20 

•' 1X2X3 

20 5 

Probability of throwine three heads is — = — 

^ ^ 64 16 

202. Compotmd Events. A certain event can happen in 
a ways and fail in b ways : a second independent event can 
happen in a' ways, and fail in b' ways, all of these ways 
being equally likely to OQCur. 

To find the probability of the simultaneous occurrence of 
the two events. 

The total number of ways in which the events can take 
place together is (a -\t b) (a' -)- b') 

(1) Both events can happen in a a' ways. 

(2) Both events can fail in b b' ways. 

(3) First event can happen and second fail in a b' ways. 

(4) First event can fail and second happen in a' b ways. 



The probability of (1) is 

The probability of (2) is 

The probability of (8) is 

The probability of (4) is 



a a' 




(a + b) (a' 


+ *') 


bb' 




{a + b) (a' 


+ V) 


ab' 




(a + b) (a' 


+ V) 


a' b 





(a + b) (a' + b') 



But the probability of the happening of the first event is 
, and of the second event is — ; — j — ^„ etc. Hence it 



a + 6' a' -\- b 

will at once be seen that the probability of the simultaneous 

occurrence of two independent events is equal to the product 
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of the probabilities of the occurrence of the component events. 

Or, in general, if Pj, P^, . . . Pn are the probabilities of 

the occurrence of any number, n, of independent events, the 

probability of the simultaneous occurrence of all the events is 

Pi X A X ... A (A) 

By independent events is meant those such that the manner 
of occurrence of one has no influence upon the manner of 
occurrence of the others. 

203. Example C. The chance that A can solve a cer- 

2 
tain problem is — , and the chance that JB can solve it is 

—. Find, 
12 

(a) The probability that both will solve it. 

(b) The probability that the problem will be solved. 

For (a). This is a question as to the probability of the 
concurrence of two independent events. Therefore by an 
application of (A), the probability that both will solve the 
problem is 



2 A _ A 

3 ^ 12 ~ 18 



For (b). The problem will be solved unless both fail. 

17 7 

The probability that both will fail is — x — = 



The probability of getting a solution is 1 



3 12 36 

J7_ _ 29 
36 ~ 36 



Example D. A pack of cards is cut, and those taken off 
then replaced. In how many trials will it be an even wager 
that an ace will be cut ? 
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Let n be the number of trials. Then n is to be found 
from 

1 



/ 48 \» _ 1 
(52J ~ "2 



where the first member of the equation represents the proba- 
bility that we shall not fail n times in succession. 

Solving for n, 

log 2 

n = 2 

log 52 - log 48 

= 8.7 

In nine trials then there is a little more than an even chance 
of cutting an ace. 

204. Dependent Events. If we have a number of events 
whose modes of occurrence are dependent one upon another, 
the probability of their concurrence will be found by the same 
method as in paragraph 202 ; a' now denoting the number 
of ways in which after the first event has happened the 
second will follow, and b' the number of ways in which after 
the first has happened the second will not follow, etc. Accord- 
ingly, the general formula (A) of paragraph 202 applies to 
dependent events as well as to independent ones. 

Also, if an event can take place in a variety of ways, the 
total probability of its occurrence will be the sum of the 
probabilities of its occurrence in each of the different ways. 

205. Mxample E. Suppose two purses contain respect- 
ively five dimes and a copper, and six dimes. A coin is taken 
at random from the first purse and placed in the second, and 
then a coin is transferred from the second to the first. What 
is the probability that the copper will remain in the first purse ? 

The probability that the copper will be taken from the first 
purse and placed in the second, and then returned to the first 
purse is 

J_ 1 _ J_ 

6" ^ Y ~~ 42 
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and the probability that the copper will not be taken from 
the first purse at all is 

5 

Therefore the probability that the copper will finally remain 
in the first purse is 

15 _ 36 _ 6^ 

42 "^ ¥ ~ 42 ~" 7 



FUNCTIONS OF SEVERAL VARIABLES. 

206. For the application of Taylor's Theorem to the 
expansion of a function of several independent variables, see 
Osborne^ s " Differential and Integral Calculus," page 145. 
And for the conditions that lead to maxima and minima 
values of such functions, see page 155 of the same work. 
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Probability of Errors. Systematic Errors. The Method of Substitu- 
tion. 

Wright, " Treatise on the Adjustment of Observations." 

Special Metliods of Solution. Applications to Geodetic and Engi- 
neering Problems. 



vi METHOD OF LEAST SQUABES. 

Merriman, " Text-Book on the Method of Least Squares." 

Chauvenet, " Treatise on the Method of Least Squares." 

DevelqpTnent of the Theory. Applications to AstroTwmieal Observa- 
tions. 

Bobek, " Lehrbuch der Ausgleichsrechnung nach der Methode 
der Kleinsten Quadrate." 

General Synopsis of the Method, illustrated by ywmerous Examples. 
Koll, " Die Methode der Kleinsten Quadrate." 

Applications to Geodesy. 
Hansen, " Von der Methode der Kleinsten Quadrate." 

Applications to Geodesy. 
Helmert, " Die Ausgleichungsrechnung nach der Methode der Klein- 
sten Quadrate." 
lAagre, " Calcul des Probabilites." 

Holman, " Discussion of the Precision of Measurements." 

Problems in Physics and Electrical Engineering. 
Weinstein, " Handbuch der Physikalischen Maassbestimmungen." 

Applications to Physical Problems. 
Oppolzer, " Lehrbuch zur Bahnbestimmung der Kometen und Plane- 
ten." 
Jordan, " Handbuch der Vermessungskunde." 

For a complete list of works on the Method of Least 
Squares published up to 1876, see 

Merriman, " A List of Writings relating to the Method of Least 
Squares, with Historical and Critical Notes." 
Published in the Transactions of the Connecticut Academy, vol. iv, 
1877. 

Notice of works published since 1876 may be found in 
periodicals devoted to the progress of Mathematical Science. 
Such as 

" Jahrbuch fiber die Fortschritte der Mathematik." 
" Bulletin des Sciences Mathfimatiques." 



TABLES. 



TABLE I. 



Values of the Integral 



=f> 



i^dt for Argument 



t 



0.4769 



a 
r 


12 3 4 


6 6 7 8 9 


Diff. 


0.0 


0.0000 0.0064 0.0108 0.0161 0.0216 


0.0269 0.0323 0.0377 0.0430 0.0484 


54 


0.1 


0538 0591 0646 0699 0762 


0806 0869 0913 0966 1020 


64 


0.2 


1073 1126 1180 1233 1286 


1339 1392 1445 1498 1561 


63 


0.3 


1603 1666 1709 1761 1814 


1866 1918 1971 2023 2076 


52 


0.4 


2127 2179 2230 2282 2334 


2386 2436 2488 2639 2690 


61 


0.6 


0.2641 0.2691 0.2742 0.2793 0.2843 


0.2893 0.2944 0.2994 0.3043 0.3093 


50 


0.6 


3143 3192 3242 3291 3340 


3389 3438 S487 3536 3583 


49 


0.7 


3632 3680 3728 3776 3823 


3870 3918 3965 4012 4069 


46 


0.8 


4106 4162 4198 4244 4290 


4336 4381 4427 4472 4517 


45 


0.9 


4662 4606 4661 4696 4739 


4783 4827 4860 4914 4967 


43 


1.0 


O.BOOO 0.5043 0.5086 0.6128 0.6170 


0.5212 0.6264 0.5295 0.6337 0.5378 


41 


1.1 


6419 6460 6600 6540 6681 


6620 6660 6700 6739 6778 


39 


1.2 


5817 5856 5894 5932 5970 


6008 6046 6083 6120 6157 


37 


1.3 


6194 6231 6267 6303 6339 


6376 6410 6445 6480 6515 


35 


1.4 


6550 6684 6618 6652 6686 


6719 6753 6786 6818 6851 


32 


1.5 


0.6883 0.691B 0.6947 0.6979 0.7011 


0.7042 0.7073 0.7104 0.7134 0.7166 


30 


1.6 


7196 7225 7266 7284 7313 


7342 7371 7400 7428 7457 


28 


1.7 


7486 7612 7540 7667 7594 


7621 7648 7676 7701 7727 


26 


1.8 


7763 7778 7804 7829 7854 


7879 7904 7928 7952 7976 


24 


1.9 


8000 8023 8047 8070 8093 


8116 8138 8161 8183 8205 


22 


2.0 


0.8227 0.8248 0.8270 0.8291 0.8312 


0.8332 0.8353 0.8373 0.8394 0.8414 


19 


2.1 


8433 8453 8473 8492 8511 


8530 8549 8667 8686 8604 


18 


2.2 


8622 8639 8667 8674 8692 


8709 8726 8742 8759 8775 


17 


2.3 


8792 8808 8824 8840 8866 


8870 8886 8901 8916 8930 


15 


2.4 


8945 8960 8974 8988 9002 


9016 9029 9043 9066 9069 


13 


2.B 


0.9082 0.9095 0.9108 0.9121 0.9133 


0.9146 0.9168 0.9170 0.9182 0.9193 


12 


2.6 


9206 9217 9228 9239 9250 


9261 9272 9283 9293 9304 


10 


2.7 


9314 9324 9334 9344 9364 


9364 9373 9383 9392 9401 


9 


2.8 


9410 9419 9428 9437 9446 


9464 9463 9471 9479 9487 


8 


2.9 


9496 9603 9611 9519 9526 


9534 9541 9648 9666 9563 


7 


3.0 


0.9670 0.9677 0.9683 0.9590 0.9597 


0.9603 0.9610 0.9616 0.9622 0.9629 


6 


3.1 


9635 9641 9647 9662 9658 


9664 9669 9675 9680 9686 


6 


3.2 


9691 9696 9701 9706 9711 


9716 9721 9726 9731 9736 


5 


3.3 


9740 9744 9749 9763 9767 


9761 9766 9770 9774 9778 


4 


3.4 


9782 9786 9789 9793 9797 


9800 9804 9807 9811 9814 


4 


3. 


0.9570 0.9636 0.9691 0.9740 0.9782 


0.9818 0.9848 0.9874 0.9896 0.9916 




4. 


9930 9943 9954 9963 9970 


9976 9981 9985 9988 9990 




B. 


9993 9994 9996 9997 9997 


9998 9998 9999 9999 9999 




00 


1.0000 






a 
r 


12 3 4 


6 6 7 8 9 


Diff. 
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TABLE 


II.- 


— Common 


Logarithms. 






n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


DifE. 


10 


0000 


0043 


0086 


0128 


0170 


0212 


0253 


0294 


0334 


0374 


42 


11 


0414 


0453 


0492 


0631 


0669 


0607 


0645 


0682 


0719 


0756 


38 


12 


0792 


0828 


0864 


0899 


0934 


0969 


1004 


1038 


1072 


1106 


36 


13 


1139 


1173 


1206 


1239 


1271 


1303 


1336 


1367 


1399 


1430 


32 


14 


1461 


1492 


1623 


1663 


1584 


1614 


1644 


1673 


1703 


1732 


30 


16 


1761 


1790 


1818 


1847 


1876 


1903 


1931 


1969 


1987 


2014 


28 


16 


2041 


2068 


2096 


2122 


2148 


2175 


2201 


2227 


2253 


2279 


27 


17 


2304 


2330 


2355 


2380 


2406 


2430 


2465 


2480 


2504 


2529 


26 


18 


2553 


2577 


2601 


2626 


2648 


2672 


2696 


2718 


2742 


2766 


24 


19 


2788 


2810 


2833 


2866 


2878 


2900 


2923 


2945 


2967 


2989 


22 


20 


3010 


3032 


3054 


3075 


3096 


3118 


3139 


3160 


3181 


3201 


21 


21 


3222 


3243 


3263 


3284 


3304 


3324 


3346 


3365 


3386 


3404 


20 


22 


3424 


3444 


3464 


3483 


3602 


3622 


3641 


3660 


3679 


3698 


19 


23 


3617 


3636 


8655 


3674 


3692 


3711 


3729 


3747 


3766 


3784 


18 


24 


3802 


3820 


3838 


3856 


3874 


3892 


3909 


3927 


3946 


3962 


18 


25 


3979 


3997 


4014 


4031 


4048 


4065 


4082 


4099 


4116 


4133 


17 


26 


4150 


4166 


4183 


4200 


4216 


4232 


4249 


4265 


4281 


4298 


17 


27 


4314 


4330 


4346 


4362 


4378 


4393 


4409 


4426 


4440 


4466 


16 


28 


4472 


4487 


4502 


4518 


4633 


4648 


4564 


4579 


4594 


4609 


15 


29 


4624 


4639 


4664 


4669 


4683 


4698 


4713 


4728 


4742 


4767 


16 


30 


4771 


4786 


4800 


4814 


4829 


4843 


4857 


4871 


4886 


4900 


14 


31 


4914 


4928 


4942 


4955 


4969 


4983 


4997 


6011 


6024 


6038 


14 


32 


6061 


5065 


5079 


5092 


5105 


5119 


6132 


5145 


6169 


5172 


13 


33 


6185 


6198 


6211 


6224 


6237 


5250 


6263 


5276 


6289 


5302 


13 


34 


5316 


6328 


5340 


5363 


6366 


5378 


6391 


6403 


6416 


5428 


13 


35 


6441 


6453 


6466 


6478 


6490 


5502 


6614 


5627 


6639 


6561 


12 


36 


6663 


6575 


6587 


5599 


6611 


6623 


5636 


5647 


5668 


6670 


12 


37 


6682 


5694 


6706 


6717 


6729 


6740 


5752 


6763 


5776 


5786 


12 


38 


5798 


6809 


6821 


6832 


5843 


6866 


6866 


6877 


6888 


6899 


11 


39 


5911 


6922 


5933 


6944 


5965 


5966 


6977 


6988 


5999 


6010 


11 


40 


6021 


6031 


6042 


6063 


6064 


6075 


6085 


6096 


6107 


6117 


11 


41 


6128 


6138 


6149 


6160 


6170 


6180 


6191 


6201 


6212 


6222 


11 


42 


6232 


6243 


6263 


6263 


6274 


6284 


6294 


6304 


6314 


6325 


10 


43 


6336 


6345 


6365 


6366 


6375 


6385 


6395 


6406 


6415 


6425 


10 


44 


6435 


6444 


6454 


6464 


6474 


6484 


6493 


6503 


6613 


6622 


10 


45 


6532 


6642 


6561 


6561 


6571 


6580 


6590 


6699 


6609 


6618 


10 


46 


6628 


6637 


6646 


6656 


6665 


6675 


6684 


6693 


6702 


6712 


9 


47 


6721 


6730 


6739 


6749 


6758 


6767 


6776 


6785 


6794 


6803 


9 


48 


6812 


6821 


6830 


6839 


6848 


6857 


6866 


6875 


6884 


6893 


9 


49 


6902 


6911 


6920 


6928 


6937 


6946 


6965 


6964 


6972 


6981 


9 


50 


6990 


6998 


7007 


7016 


7024 


7033 


7042 


7050 


7059 


7067 


9 


61 


7076 


7084 


7093 


7101 


7110 


7118 


7126 


7135 


7143 


7162 


8 


62 


7160 


7168 


7177 


7186 


7193 


7202 


7210 


7218 


7226 


7236 


8 


53 


7243 


7251 


7269 


7267 


7275 


7284 


7292 


7300 


7308 


7316 


8 


54 


7324 


7332 


7340 


7348 


7356 


7364 


7372 


7380 


7388 


1396 


8 
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1 


2 


3 


4 


6 


6 


7 


8 


9 
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TABLE II.- 


— Common 


Logarithms. 






n 





1 


2 


3 


4 


6 


6 


7 


8 


9 


Difl. 


55 


7404 


7412 


7419 


7427 


7435 


7443 


7451 


7459 


7466 


7474 


8 


56 


7482 


7490 


7497 


7505 


7613 


7620 


7628 


7536 


7643 


7661 




57 


7659 


7566 


7574 


7682 


7689 


7697 


7604 


7612 


7619 


7627 




68 


7634 


7642 


7649 


7657 


7664 


7672 


7679 


7686 


7694 


7701 




69 


7709 


7716 


7723 


7731 


7738 


7746 


7762 


7760 


7767 


7774 




60 


7782 


7789 


7796 


7803 


7810 


7818 


7825 


7832 


7839 


7846 


7 


61 


7863 


7860 


7868 


7876 


7882 


7889 


7896 


7903 


7910 


7917 




62 


7924 


7931 


7938 


7946 


7962 


7069 


7966 


7973 


7980 


7987 




63 


7993 


8000 


8007 


8014 


8021 


8028 


8035 


8041 


8048 


8066 




64 


8062 


8069 


8075 


8082 


8089 


8096 


8102 


8109 


6116 


8122 




65 


8129 


8136 


8142 


8149 


8166 


8162 


8169 


8176 


8182 


8189 


7 


66 


8195 


8202 


8209 


8215 


8222 


8228 


8235 


8241 


8248 


8264 




67 


8261 


8267 


8274 


8280 


8287 


8293 


8299 


8306 


8312 


8319 




68 


8325 


8331 


8338 


8344 


8351 


8357 


8363 


8370 


8376 


8382 




69 


8388 


8395 


8401 


8407 


8414 


8420 


8426 


8432 


8439 


8445 




70 


8461 


8457 


8463 


8470 


8476 


8482 


8488 


8494 


8600 


8606 


6 


71 


8613 


8619 


8526 


8631 


8637 


8543 


8649 


8665 


8661 


8567 




72 


8673 


8679 


8685 


8691 


8697 


8603 


8609 


8615 


8621 


8627 




73 


8633 


8639 


8646 


8661 


8667 


8663 


8669 


8675 


8681 


8686 




74 


8692 


8698 


8704 


8710 


8716 


8722 


8727 


8733 


8739 


8745 




75 


8761 


8756 


8762 


8768 


8774 


8779 


8786 


8701 


8797 


8802 


6 


76 


8808 


8814 


8820 


8826 


8831 


8837 


8842 


8848 


8864 


8859 




77 


8865 


8871 


8876 


8882 


8887 


8893 


8899 


8904 


8910 


8915 




78 


8921 


8927 


8932 


8938 


8943 


8949 


8964 


8960 


8965 


8971 




79 


8076 


8982 


8987 


8993 


8998 


9004 


9009 


9016 


9020 


9026 




80 


9031 


9036 


9042 


9047 


9053 


9068 


9063 


9069 


9074 


9079 


5 


81 


9085 


9090 


9096 


9101 


9106 


9112 


9117 


9122 


9128 


0133 




82 


9138 


0143 


9149 


9164 


9169 


9166 


9170 


9176 


9180 


9186 




83 


9191 


9196 


9201 


9206 


9212 


9217 


9222 


9227 


9232 


9238 




84 


9243 


0248 


9263 


9268 


9263 


0269 


9274 


9279 


9284 


9289 




85 


9294 


9299 


9304 


9309 


9316 


0320 


9325 


9330 


9335 


9340 


S 


86 


9345 


0360 


9365 


9360 


9365 


9370 


9375 


9380 


9385 


9390 




87 


9395 


9400 


9406 


9410 


9415 


9420 


0426 


9430 


9436 


0440 




88 


9445 


0450 


9465 


9460 


9466 


9469 


9474 


9479 


9484 


9489 




89 


9494 


9499 


9504 


9609 


9613 


9618 


9623 


9628 


9633 


9638 




90 


9542 


9547 


9662 


9557 


9662 


9666 


9671 


9576 


9681 


9686 


5 


91 


9590 


9595 


9600 


9605 


9609 


9614 


9619 


9624 


9628 


9633 




92 


9638 


9643 


9617 


0652 


0657 


9661 


9666 


9671 


9675 


9680 




93 


9685 


9689 


9694 


9699 


9703 


9708 


9713 


9717 


9722 


9727 




94 


9731 


9736 


9741 


9745 


9760 


9764 


9769 


9763 


0768 


9773 




96 


9777 


9782 


9786 


9791 


9795 


9800 


9805 


9809 


0814 


9818 


4 


96 


9823 


9827 


9832 


0836 


9841 


9845 


0860 


9864 


9869 


9863 




97 


9868 


9872 


9877 


9881 


9886 


9890 


9804 


0809 


9903 


9908 




98 


9912 


9917 


9921 


9926 


9930 


9934 


9939 


9943 


9948 


9952 




99 


9966 


9961 


9966 


9969 


0074 


9978 


9983 


9987 


9991 


9996 




™ 





1 


2 


3 


4 


5 


6 


7 


8 


9 
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TABLE 


III.- 


— Squares 


of Numbers. 






« 





1 


2 


3 


4 


5 


6 


7 


8 


9 


DUE. 


1.0 


1.000 


1.020 


1.040 


1.061 


1.082 


1.103 


1.124 


1.146 


1.166 


1.188 


22 


1.1 


1.210 


1.232 


1.254 


1.277 


1.300 


1.323 


1.346 


1.369 


1.392 


1.416 


24 


1.2 


1.440 


1.464 


1.488 


1.513 


1.638 


1.563 


1.588 


1.613 


1.638 


1.664 


26 


1.3 


1.690 


1.716 


1.742 


1.769 


1.796 


1.823 


1.850 


1.877 


1.904 


1.932 


28 


1.4 


1.960 


1.988 


2.016 


2.045 


2.074 


2.103 


2.132 


2.161 


2.190 


2.220 


30 


1.5 


2.260 


2.280 


2.310 


2.341 


2.372 


2.403 


2.434 


2.465 


2.496 


2.628 


32 


1.6 


2.660 


2.592 


2.621 


2.657 


2.690 


2.723 


2.756 


2.789 


2.822 


2.856 


34 


1.7 


2.890 


2.924 


2.968 


2.993 


3.028 


3.063 


3.098 


3.133 


3.168 


3.204 


36 


1.8 


3.240 


3.276 


3.312 


3.349 


3.386 


3.423 


3.460 


3.497 


3.634 


3.572 


38 


1.9 


3.610 


3.648 


3.686 


3.725 


3.764 


3.803 


3.842 


3.881 


3.920 


3.960 


40 


2.0 


4.000 


4.040 


4.080 


4.121 


4.162 


4.203 


4.244 


4.286 


4.326 


4.368 


42 


2.1 


4.410 


4.462 


4.494 


4.537 


4.580 


4.623 


4.666 


4.709 


4.752 


4.796 


44 


2.2 


4.840 


4.884 


4.928 


4.973 


6.018 


5.063 


6.108 


5.153 


5.198 


5.244 


46 


2.3 


5.290 


6.336 


6.382 


6.429 


5.476 


6.523 


6.570 


6.617 


6.664 


6.712 


48 


2.4 


6.760 


5.808 


6.866 


5.906 


5.954 


6.003 


6.052 


6.101 


6.150 


6.200 


50 


2.5 


6.260 


6.300 


6..'550 


6.401 


6.462 


6.503 


6.564 


6.606 


6.6S6 


6.708 


52 


2.6 


6.760 


6.812 


6.864 


6.917 


6.970 


7.023 


7.076 


7.129 


7.182 


7.236 


64 


2.T 


7.290 


7.344 


7.398 


7.453 


7.508 


7.663 


7.618 


7.673 


7.728 


7.784 


66 


2.8 


7.840 


7.896 


7.962 


8.009 


8.066 


8.123 


8.180 


8.237 


8.294 


8.352 


58 


2.9 


8.410 


8.468 


8.526 


8.586 


8.644 


8.703 


8.762 


8.821 


8.880 


8.940 


60 


3.0 


9.000 


9.060 


9.120 


9.181 


9.242 


9.303 


9.364 


9.425 


9.486 


9.548 


62 


3.1 


9.610 


9.672 


9.734 


9.797 


9.860 


9.923 


9.986 


10.06 


10.11 


10.18 


6 


3.2 


10.24 


10.30 


10.37 


10.43 


10.50 


10.66 


10.63 


10.69 


10.76 


10.82 


7 


3.3 


10.89 


10.96 


11.02 


11.09 


11.16 


11.22 


11.29 


11.36 


11.42 


11.49 


7 


3.4 


11.56 


11.63 


11.70 


11.76 


11.83 


11.90 


11.97 


12.04 


12.11 


12.18 


7 


3.5 


12.25 


12.32 


12.39 


12.46 


12.53 


12.60 


12.67 


12.74 


12.82 


12.89 


7 


3.6 


12.96 


13.03 


13.10 


13.18 


13.25 


13.32 


13.40 


13.47 


13.64 


13.62 


7 


3.7 


13.69 


13.76 


13.84 


13.91 


13.99 


14.06 


14.14 


14.21 


14.29 


14.36 


8 


3.8 


14.44 


14.52 


14.59 


14.67 


14.75 


14.82 


14.90 


14.98 


15.05 


15.13 


8 


3.9 


16.21 


16.29 


15.37 


16.44 


15.52 


15.60 


16.68 


16.76 


15.84 


15.92 


8 


4.0 


16.00 


16.08 


16.16 


16.24 


16.32 


16.40 


16.48 


16.66 


16.65 


16.73 


8 


4.1 


16.81 


16.89 


16.97 


17.06 


17.14 


17.22 


17.31 


17.39 


17.47 


17.56 


8 


4.2 


17.64 


17.72 


17.81 


17.89 


17.98 


18.06 


18.16 


18.23 


18.32 


18.40 


9 


4.3 


18.49 


18.58 


18.66 


18.76 


18.84 


18.92 


19.01 


19.10 


19.18 


19.27 


9 


4.4 


19.36 


19.46 


19.64 


19.62 


19.71 


19.80 


19.89 


19.98 


20.07 


20.16 


9 


4.5 


20.25 


20.34 


20.43 


20.62 


20.61 


20.70 


20.79 


20.88 


20.98 


21.07 


9 


4.6 


21.16 


21.25 


21.34 


21.44 


21.53 


21.62 


21.72 


21.81 


21.90 


22.00 


9 


4.7 


22.09 


22.18 


22.28 


22.37 


22.47 


22.66 


22.66 


22.76 


22.86 


22.94 


10 


4.8 


23.04 


23.14 


23.23 


23.33 


23.43 


23.62 


23.62 


23.72 


23.81 


23.91 


10 


4.9 


24.01 


24.11 


24.21 


24.30 


24.40 


24.60 


24.60 


24.70 


24.80 


24.90 


10 


6.0 


26.00 


26.10 


25.20 


26.30 


26.40 


25.60 


26.60 


26.70 


25.81 


25.91 


10 


5.1 


26.01 


26.11 


26.21 


26.32 


26.42 


26.52 


26.63 


26.73 


26.83 


26.94 


10 


5.2 


27.04 


27.14 


27.25 


27.35 


27.46 


27.66 


27.67 


27.77 


27.88 


27.98 


11 


5.3 


28.09 


28.20 


28.30 


28.41 


28.62 


28.62 


28.73 


28.84 


28.94 


29.06 


11 


5.4 


29.16 


29.27 


29.38 


29.48 


29.69 


29.70 


29.81 


29.92 


30.03 


30.14 


11 


n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


DifE. 



TABLES. 



XI 







TABLE 


III.- 


— Squares 


of Numbers. 






n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Diff. 


6.5 


30.25 


30.36 


30.47 


30.68 


30.69 


30.80 


30.91 


31.02 


31.14 


31.26 


11 


5.6 


31.36 


31.47 


31.68 


31.70 


31.81 


31.92 


32.04 


32.16 


32.26 


32.38 


11 


6.7 


32.49 


32.60 


32.72 


32.83 


32.95 


33.06 


33.18 


33.29 


33.41 


33.62 


12 


5.8 


33.64 


33.76 


33.87 


33.99 


34.11 


34.22 


34.34 


34.46 


34.67 


34.69 


12 


5.9 


34.81 


34.93 


36.06 


36.16 


35.28 


36.40 


36.62 


SSM 


35.76 


35.88 


12 


6.0 


36.00 


36.12 


36.24 


36.36 


36.48 


36.60 


36.72 


36.84 


36.97 


37.09 


12 


6.1 


37.21 


37.33 


37.45 


37.68 


87.70 


37.82 


37.96 


38.07 


38.19 


38.32 


12 


6.2 


38.44 


38.66 


38.69 


38.81 


33.94 


39.06 


39.19 


39.31 


39.44 


39.66 


13 


6.3 


39.69 


39.82 


39.94 


40.07 


40.20 


40.32 


40.45 


40.68 


40.70 


40.83 


13 


6.4 


40.96 


41.09 


41.22 


41.34 


41.47 


41.60 


41.73 


41.86 


41.99 


42.12 


13 


6.5 


42.26 


42.38 


42.61 


42.64 


42.77 


42.90 


43.03 


43.16 


43.30 


43.43 


13 


6.6 


43.66 


43.69 


43.82 


43.96 


44.09 


44.22 


44.36 


44.49 


44.62 


44.76 


13 


6.7 


44.89 


45.02 


46.16 


45.29 


46.43 


46.66 


46.70 


46.83 


46.97 


46.10 


14 


6.8 


46.24 


46.38 


46.61 


46.65 


46.79 


46.92 


47.06 


47.20 


47.33 


47.47 


14 


6.9 


47.61 


47.75 


47.89 


48.02 


48.16 


48.30 


48.44 


48.68 


48.72 


48.86 


14 


7.0 


49.00 


49.14 


49.28 


49.42 


49.56 


49.70 


49.84 


49.98 


50.13 


60.27 


.14 


7.1 


60.41 


50.55 


60.69 


50.84 


50.98 


61.12 


51.27 


51.41 


61.65 


61.70 


14 


7.2 


51.84 


61.98 


62.13 


62.27 


62.42 


52.66 


62.71 


62.86 


63.00 


53.14 


16 


7.3 


63.29 


63.44 


63.68 


53.73 


53.88 


64.02 


64.17 


64.32 


64.46 


64.61 


15 


7.4 


64.76 


54.91 


66.06 


65.20 


55.35 


55.50 


56.65 


65.80 


66.95 


56.10 


15 


7.5 


66.25 


56.40 


66.66 


56.70 


66.86 


67.00 


57.15 


67.30 


57.46 


57.61 


15 


7.6 


67.76 


57.91 


68.06 


58.22 


68.37 


58.52 


58.68 


68.83 


68.98 


69.14 


15 


7.7 


69.29 


59.44 


69.60 


59.75 


59.91 


60.06 


60.22 


60.37 


60.63 


60.68 


16 


7.8 


60.84 


61.00 


61.15 


61.31 


61.47 


61.62 


61.78 


61.94 


62.09 


62.25 


16 


7.9 


62.41 


62.57 


62.73 


62.88 


63.04 


63.20 


63.36 


63.62 


63.68 


63.84 


16 


8.0 


64.00 


64.16 


64.32 


64.48 


64.64 


64.80 


64.96 


66.12 


65.29 


65.46 


16 


8.1 


65.61 


66.77 


65.93 


66.10 


66.26 


66.42 


66.59 


66 75 


66.91 


67.08 


16 


8.2 


67.24 


67.40 


67.57 


67.73 


67.90 


68.06 


68.23 


68.39 


68.66 


68.72 


17 


8.3 


68.89 


69.06 


69.22 


69.39 


69.56 


69.72 


69.89 


70.06 


70.22 


70.39 


17 


8.4 


70.66 


70.73 


70.90 


71.06 


71.23 


71.40 


71.67 


71.74 


71.91 


72.08 


17 


8.5 


72.26 


72.42 


72.59 


72.76 


72.93 


73.10 


73.27 


73.44 


73.62 


73.79 


17 


8.6 


73.96 


74.13 


74.30 


74.48 


74.65 


74.82 


75.00 


76.17 


76.34 


76.62 


17 


8.7 


76.69 


75.86 


76.04 


76.21 


76.39 


76.66 


76.74 


76.91 


77.09 


77.26 


18 


8.8 


77.44 


77.62 


7T.79 


77.97 


78.15 


78.32 


78.60 


78.68 


78.85 


79.03 


18 


8.9 


79.21 


79.39 


79.67 


79.74 


79.92 


80.10 


80.28 


80.46 


80.64 


80.82 


18 


9.0 


81.00 


81.18 


81.36 


81.64 


81.72 


81.90 


82.08 


82.26 


82.45 


82.63 


18 


9.1 


82.81 


82.99 


83.17 


83.36 


83.54 


83.72 


83.91 


84.09 


84.27 


84.46 


18 


9.2 


84.64 


84.82 


86.01. 


85.19 


85.38 


85.66 


85.76 


85.93 


86.12 


86.30 


19 


9.3 


86.49 


86.68 


86.86 


87.06 


87.24 


87.42 


87.61 


87.80 


87.98 


88.17 


19 


9.4 


88.36 


88.65 


88.74 


88.92 


89.11 


89.30 


89.49 


89.68 


89.87 


90.06 


19 


9.6 


90.25 


90.44 


90.63 


90.82 


91.01 


91.20 


91.39 


91.58 


91.78 


91.97 


19 


9.6 


92.16 


92.36 


92.64 


92.74 


92.93 


93.12 


93.32 


93.61 


93.70 


93.90 


19 


9.7 


94.09 


94.28 


94.48 


94.67 


94.87 


95.06 


95.26 


95.46 


95.65 


95.84 


20 


9.8 


96.04 


96.^ 


96.43 


96.63 


96.83 


97.02 


97.22 


97.42 


97.61 


97.81 


20 


9.9 


98.01 


98.21 


98.41 


98.60 


98.80 


99.00 


99.20 


99.40 


99.60 


99.80 


20 


n 





1 


2 


3 


4 


6 


6 


7 


8 


9 


Diff. 



